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We study a junction of three quantum wires enclosing a magnetic fiux. The wires are modeled 
as single-channel spinless Tomonaga-Luttinger liquids. This is the simplest problem of a quantum 
junction between Tomonaga-Luttinger liquids in which Fermi statistics enter in a non-trivial way. We 
study the problem using a mapping onto the dissipative Hofstadter model, describing a single particle 
moving on a plane in a magnetic field and a periodic potential coupled to a harmonic oscillator bath. 
Alternatively we study the problem by identifying boundary conditions corresponding to the low 
energy fixed points. We obtain a rich phase diagram including a chiral fixed point in which the 
asymmetric current flow is highly sensitive to the sign of the flux and a phase in which electron 
pair tunneling dominates. We also study the effects on the conductance tensor of the junction of 
contacting the three quantum wires to Fermi liquid reservoirs. 

PACS numbers: 



I. INTRODUCTION 

While electron-electron interactions can, in many 
cases, be adequately described by Fermi liquid theory in 
two- or three-dimensional systems, they lead to more ex- 
otic effects in one-dimensional Tomonaga-Luttinger liq- 
uids (TLL). One context in which striking TLL behavior 
occurs is in the response of a quantum wire ilj or a spin 
chain to a single constriction or weak link. The con- 
ductance through the constriction scales to zero at low 
temperatures in the case of repulsive interactions, corre- 
sponding to a TLL parameter g < 1 or to the ideal value 
of ge'^/h for attractive interactions, g > I. TLL behav- 
ior has recently been studied experimentally in carbon 
nano-tubes. 

In order to construct a useful circuit out of quantum 
wires it will be necessary to incorporate junctions of three 
or more wires. These and other closely related prob- 
lems exhibit rather rich TLL effects which have been the 
subject of a number of recent pap ers. p,p , 0, S S 
lIlEllIllIllIlllllillillfllll l3"^Much of the 
work in this field has used bosonization. As emphasized 
by Nayak et al. 0, when the number of wires meet- 
ing at a junction exceeds two, the Klein factors which 
give bosonized operators Fermi statistics play a crucial 
role. We recently introduced a new method to study 
this problem, mapping it into the dissipative Hofstadter 
model (DHM), which describes a single particle moving 
in a uniform magnetic field and a periodic potential in 
two dimensions and coupled to a bath of harmonic os- 
cillators. This mapping is useful because it allows us to 
take adva ntag e of earlier results of Callan and Freed on 
the DHM 1231 . When the three quantum wires enclose a 



magnetic flux, the mapping to the DHM also allows us 
to identify a new low energy chiral fixed point with an 
asymmetric flow of current that is highly sensitive to the 
sign of the flux. 

The purpose of this paper is to present a comprehensive 
study of the three- wire junction (or Y-junction) problem. 
We describe the mapping to the DHM and its implica- 
tions in more detail, and we further study the problem 
using a complimentary method which involves identify- 
ing boundary conditions on the boson fields which corre- 
spond to low energy fixed points. In combination, these 
methods provide a coherent picture of the possible stable 
fixed point conductances of the three-wire junction, and 
this picture shows a regime with chiral current flow and a 
regime in which electron pair tunneling dominates. This 
paper also addresses the problem of how the fixed point 
conductance tensor of the three- wire junction is modified 
if the TLL wires are connected to Fermi liquid leads far 
from the junction. 

The paper is organized as follows. In Sec.^we provide 
a brief summary of the results for the conductance ten- 
sor of the Y-junction and for the Renormalization Group 
(RG) flow diagrams for different ranges of the Luttinger 
parameter g. In Sec. lIIII we present our effective model for 
the three-wire junction. In particular, we show explicitly 
how the magnetic flux dependence comes into the prob- 
lem in which the three wires are connected to a circular 
ring by taking into account the discrete energy levels on 
the ring. In Sec. II VI we review the bosonization of the Y- 
jmiction model. In Sec. we discuss the strong coupling 
limit of the junction that is obtained by simply taking 
the hopping between the three wires to infinity. We show 
that such a simple argument does not provide informa- 



tion on the nature of the fixed points for a range of the 
parameter g, and that different methods that we describe 
in the paper are needed to shed light onto the nature of 
the strong coupUng fixed points. In Sec. I VII we estabhsh 
the mapping of the Y-junction model onto the DHM and 
deduce its consequences for the phase diagram of the Y- 
junction model. In Sec. IVIII we introduce the study of 
the system by looking at the boundary conditions satis- 
fied by the bosonic fields at the junction. In Sec lVIIll we 
discuss the twisted structure of the Hilbert space that is 
manifest in the compactification of the bosons so as to 
take proper care of the fermionic statistics of the elec- 
trons in the wires. In Sec. IIXI we introduce what we 
refer to as the method of delayed evaluation of boundary 
conditions (DEBC), which is less reliant on the appara- 
tus of BCFT. In Sec. |3 we apply all the results from 
the previous sections to analyze the RG fixed points in 
the junction problem. We also show in this section that 
asymmetries in the couplings between the three wires are 
irrelevant at low energies and temperatures, so the sim- 
ple Z3 symmetric model of the system may indeed be a 
good description of realistic junctions. In Sec. lXII we use 
simple arguments based on energy conservation to obtain 
constraints on the conductance tensor. We show that all 
the fixed points that we are able to understand in this pa- 
per correspond to physical situations where no energy is 
transfered to neutral (non-charge carrying) modes in the 
quantum wires. In Sec. IXIll we discuss how the conduc- 
tance tensors are modified if the Luttinger liquid wires 
are connected to Fermi liquid leads far from the junc- 
tion. Depending on details of the potential experiment, 
this may be a better model. In Sec. IXIIII we present 
and discuss open problems. In Appendix 1X1 we give re- 
sults on the non-interacting version of our model. In 
Appendix^ we review boundary conformal field theory. 
In Appendix lUl we review the boundary conformal field 
theory of standard free bosons. In Appendix^we review 
the boundary conformal field theory approach to calcu- 
lating the conductance. In Appendix El we discuss the 
related problem of a Y-junction of quantum Hall edge 
states. 



II. SUMMARY OF RESULTS 

We summarize in this section our results for the con- 
ductance tensor for the three wire junction, as well as our 
conjectured RG flow diagrams for several different ranges 
of the interaction parameter g. These results are based 
on our findings from three different methods discussed in 
the paper: mapping to the dissipative Hofstadter model, 
boundary conformal field theory, and delayed evaluation 
of boundary conditions. Our results rely on various as- 
sumptions and arguments which we believe to be reason- 
able in the light of several supporting evidences. The 
details of our analysis and assumptions are given in later 
sections. 

The device we consider is shown in Fig. ^ where the 



2 




FIG. 1: Junction of three quantum wires with a magnetic flux 
threading the ring. The Vi,2,3 are the voltages applied on each 
wire, and the /i,2,3 the currents arriving at the junction from 
each of the three wires. 

three quantum wires are connected to a ring which can be 
threaded by a magnetic flux. The electron transfer pro- 
cesses between the three wires meeting at the Y-junction 
is modeled by the effective Hamiltonian 

3 

Hb = -Y}^ e^^/^l (0)^,-1 (0) + h.c.] , (2.1) 

where ijjj [x = 0) is the electron operator at the endpoint 
of the j-th [j — 1, 2, 3, and we identify j = j + 3) quan- 
tum wire that connects to the junction (a; = 0). F is 
the tunneling amplitude between the leads, and (f) is a 
phase that depends on the flux that threads the ring. 
In addition to the boundary term Eq. H2.1|l . there is a 
bulk Hamiltonian -ffbuik that models the Luttinger liquid 
wires characterized by the interaction parameter g. The 
derivation of this effective Hamiltonian at the junction, 
and the effective Hamiltonian for the wire part, is given 
in section HTll 

The presence of a magnetic flux breaks generically 
time-reversal invariance, with the exception of = 0, tt, 
in which case time-reversal symmetry exists. In this pa- 
per, we assume the three wires to be identical. Never- 
theless, the junction can be asymmetric if the wires are 
attached to the ring in different ways. If all three wires 
are attached identically to the ring, the junction is sym- 
metric under cyclic permutations of the wires; this is the 
Z3 symmetric case. 

The important physical quantity for the three-wire 
junction problem is the conductance tensor. Within the 
linear response theory, the total current /-,- flowing into 
the junction from wire j is related to the voltage Vk ap- 
plied to wire k by 

k 

where j,k = 1,2,3 and Gjk is the 3x3 conductance 
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tensor. Note that current conservation implies that: 

^/,=0. (2.3) 

3 

Furthermore, a common voltage applied to all three wires 
results in zero current. Thus: 

5]G,fe=^G,,=0. (2.4) 

j k 

For a Z3 symmetric junction, the conductance tensor 
takes the form 

G,fe = ^(35,.-l) + ^e,fc, (2.5) 

where we separate the symmetric and anti-symmetric 
components of the tensor, and Gs and Ga are scalar 
conductances. (The e^j are defined as follows: £12 = 
£23 = €31 = 1, £21 = £32 = £13 = -1 and ejj = 0.) The 
anti-symmetric component, controlled by Ga, is present 
only when time-reversal symmetry is broken by a mag- 
netic flux. Namely, Ga should vanish in the time-reversal 
symmetric case (j) — 0,tt. Even when time-reversal sym- 
metry is broken in the microscopic model, Ga might van- 
ish in the low-energy fixed point if the RG flow restores 
the time-reversal symmetry. 

From the tensor conductance Eq. (|2.5() . one reads 

Gs — Gil = G22 = G33. (2.6) 

This represents the conductance of each wire when zero 
voltage is applied to the other two wires. For this reason, 
we shall sometimes refer to Gs in the text as the single 
terminal conductance. 

We also consider the same three-wire junction device 
when connected to Fermi liquid reservoirs, as depicted in 
Fig. 121 (the wide regions at the ends of the wires represent 
the Fermi liquid contacts). It is known that the contact 
to the reservoirs affect the two-terminal conductance of 
single quantum wires This renormalization of 

the conductance also takes place in the three-wire junc- 
tions, and we calculate these renormalization effects in 
this paper. In the presence of the Fermi liquid leads, 
we obtain that the "dressed" 3x3 conductance tensor 
G = G^/ loads of the system is related to the "bare" con- 
ductance tensor G = G^/o leads by 

, ^ G-] , , +G-H , (2.7) 

w/ leads w/o leads ' c ' \ / 

where Gc = "^g/ig — 1) e'^/h is effectively a contact con- 
ductance between the leads and the quantum wires. This 
relationship in Eq. H2.7|l can be interpreted physically 
as the addition in series of the lead/ wire interface resis- 
tances to the Y-junction resistance tensor. (Notice that, 
because of the J2k ^jk = condition, the inverse G~^ 
does not exist. Hence one must work directly with the 
conductance instead of the resistance tensor as discussed 
in section IXIII l 




FIG. 2: The same three-wire junction as in Fig. Q but now 
with the Fermi liquid reservoirs attached at the ends of the 
three wires. The ^^1,2,3 are now the voltages as measured at 
the three reservoirs. 

We are mostly interested in the stable RG fixed points 
of the problem. They would govern the physics in the 
low-energy limit ~ namely in the limit of low bias voltages 
and low temperature. Below we list our results for the 
low-energy fixed point conductances for the Y-junction 
devices. Which fixed point is the stable one is deter- 
mined by the Luttinger parameter g that depends on the 
strength of local electron-electron interactions. We find 
that the stable fixed point has the Z3 symmetry in most 
cases. Namely, even if the Z3 symmetry is broken at 
the junction, the Z3 symmetry would be restored in the 
low-energy limit, except for a few special cases. Thus we 
primarily focus on the RG flow in the presence of the Z3 
symmetry. The most important parameters of the junc- 
tion are then the hopping strength F and the enclosed 
magnetic flux (p. If we understand the RG flow most 
simply by the effective values of F and varying as func- 
tions of the energy scale, the RG flow diagram could be 
drawn on the T-(j) plane as shown in the following. 

As in the case of junctions of two wires, the interaction 
parameter g controls the RG flow and dictates the phase 
diagram. Now we present our results for several different 
ranges of g. 

A. g< 1 

When the interaction in the quantum wires is repulsive 
{g < 1), the hopping amplitude F decreases along the RG 
flow. The stable fixed point corresponds to F = 0, namely 
to completely decoupled wires, as pictorially illustrated 
in Fig. O We will call this fixed point the N fixed point, 
where N stands for Neumann boundary condition. We 
note that the value of the effective flux (f> does not matter 
in the decoupled limit, so in the N fixed point there is no 
breaking of time reversal. Moreover, even if the junction 
is Z3 asymmetric, the system is generically renormalized 
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FIG. 3: Pictorial representation of the N fixed point. It just 
corresponds to three decoupled wires. 
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FIG. 4: RG flow diagram for g < 1. The junction is renor- 
malized toward the decoupled N fixed point. 

into the N fixed point; the Z3 asymmetry turns out to 
be irrelevant. The leading irrelevant perturbation to the 
N fixed point is the electron hopping between two wires, 
which has scaling dimension > 1. 

The conductance tensor at the N fixed point is obvi- 
ously 

Gjk = Gjk = 0. (2.8) 

Of course there is no difference between the "bare" con- 
ductance G and the "dressed" (by the reservoirs) con- 
ductance G at the N fixed point. 

The corresponding fiow diagram in the T-(f) plane is 
shown in Fig. 01 

B. g = l 

For the noninteracting case g = 1, one can solve the 
original electron model as a single particle problem. The 
junction is then characterized by the 3x3 scattering ma- 
trix of a single free electron. Thus there is a continuous 



FIG. 5: Pictorial representation of the chiral fixed point. The 
incoming electrons from one wire is diverted to one of the 
other wires. 



manifold of RG fixed points. In other words, electron 
hopping is exactly marginal. 

Several special cases on the "free electron" manifold 
are worth mentioning. Let us consider the Z3 symmetric 
case first. The N fixed point, as introduced above for 
g < 1, corresponds to a complete refiection for each wire 
and thus is a special point on the RG fixed manifold. If 
we further impose the time reversal symmetry, we always 
have Ga — while the maximal single-terminal conduc- 
tance on the "free electron" manifold is Gs — 8/9 e'^/h. 
The constraint on the maximum of Gs in the time rever- 
sal symmetric case is due to the unitarity of the single 
electron scattering matrix as discussed in Ref. 0. On 
the other hand, if the time reversal symmetry is broken 
by the flux (f>, the maximal single-terminal conductance 
Gs ~ e'^/h can be realized by a complete transmission 
of electrons from wire j to wire j + 1 (or j — 1), where 
again we identify j = j + 3. We call them chiral fixed 
points x±; they are pictorially represented in Fig.|Sl and 
they are also on the RG fixed manifold. At the x± points 
Ga = ±65. The conductance tensor at x± reads 

Gjk = ^ [{5Sjk - 1) ± ejTc] = -J^i^jk - Sj.k±i)- (2.9) 

If we allow breaking of the Z3 symmetry, a simple fixed 
point is given by connecting perfectly two of the wires, 
while the other is left decoupled completely, as depicted 
in Fig. El This fixed point, which we call Da, is also 
described by a single electron scattering matrix and thus 
a special point on the "free electron" manifold. It is 
clear that there are other fixed points on the manifold 
with varying degrees of Z3 asymmetry. 

In fact, even for noninteracting wires or 17 = 1, other 
RG fixed points, which cannot be characterized by the 
single electron scattering matrix, exist in the presence 
of the interaction at the junction. However, these fixed 
points are more unstable, and the RG flows (shown in 
Fig. [TJ are towards the "free electron" manifold. 
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FIG. 6: Pictorial representation of the asymmetric fixed point 
Da, in which two of the wires (for example wires 1 and 2 in 
the figure) are perfectly connected while the other (wire 3) is 
left decoupled. This fixed point is unstable for all values of g 
with the exception of gr = 1, 3, in which case Da may belong 
to the continuous manifold of fixed points. 
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FIG. 7: RG fiow diagram for g = 1. There is a stable RG 
fixed manifold which includes N and x± ^ special points. 



tively. The chiral fixed points exhibit the asymmetric 
conductance tensor 

G%^^ [{?.5,u-l)±9^,k] (2.10) 

with the single terminal conductance or symmetric com- 
ponent 



and anti-symmetric component Ga — i<? Gg. We note 
that, in the limit 5 ^ 1, the chiral x± fixed points reduce 
to the complete transmission of single electron from wire 
j to wire j ± 1 discussed previously. 

An intriguing aspect of the conductance (|2.1()|l is that, 
when a voltage is applied to one of the wires while the 
other two are kept at zero voltage (as illustrated in 
Fig. IHt), the current is rather "sucked in" from one of 
the wires with the zero voltage. However, if we include 
the non-interacting leads to the system (as shown now 
in Fig. 13)), the effective conductance tensor for x± ^-re 
given by the conductance tensor Eq. (|2.9|) obtained for 
X± at (7 = 1 (in which Ga — ±^5 with Gs — e^/h). 
Namely, the "sucking" effect disappears in the presence 
of non-interacting leads, while the asymmetry holds. 

The RG flow diagram (shown in Fig. ^ implies that 
even a tiny flux (p brings the junction to the completely 
asymmetric conductance in the low-energy limit. The 
leading irrelevant perturbation to the x± fixed points 
has scaling dimension 4(7/(3 -I- g^) > 1. A remarkable 
aspect of the fixed points x± is that the conductance, as 
well as the scaling dimension, exhibit non-monotonic de- 
pendence on the interaction parameter g, unlike in other 
known applications of TLL. Making the electron interac- 
tion more attractive can decrease the conductance! 

We note that the Da fixed point, which breaks the 
Z3 symmetry, is unstable for this range of g. Hence, 
small asymmetries in the physical device do not affect 
the conductance of the x± fixed points at low voltages 
and temperatures for this range of the g parameter. 

D. g = 3 



C. 1< < 3 

The N fixed point is now unstable against the inclusion 
of electron hoppings, which have scaling dimension 1/g < 
1. For the time-reversal symmetric case (j) — 0,tt, the 
system is renormalized into the M fixed point. The M 
fixed point should have the single-terminal conductance 
Gs in the range < Gs < (4(?/3) (e^/h), as shown in 
Section lXll The M fixed point is unstable against turning 
on the fiux (j). Unfortunately, we do not have an exact 
solution corresponding to the M fixed point, and other 
properties are not known. 

For < < TT and — tt < (j> < 0, the RG flows are to- 
wards the stable, chiral fixed points x+ and x~ , respec- 



The N fixed point is again unstable. On the other 
hand, there is an RG fixed manifold (shown in Fig. I10|) . 
as we have seen for 17 = 1. For g = the manifold 
is again characterized by the scattering matrix of a free 
fermion, which is not the original electron. The chiral 
fixed points x± become the special points on the manifold 
as g — > 3. The largest single-terminal conductance is 
realized on another special point Dp on the manifold, 
which is explained in the next subsection for 3 < 17 < 9. 

The system with (j) ^ 0, tt is renormalized into the fixed 
manifold. On the other hand, for the time-reversal sym- 
metric case = 0, TT, the infrared M fixed point is not on 
the manifold, and it is unstable against the inclusion of 
the flux (f). 
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FIG. 8: Response of the Y-junction system at the x+ fixed 
point when a voltage V is applied to wire 1, while wires 2 
and 3 are grounded. Case (a) is when wire 1 is directly at a 
voltage V , and case (b) is when wire 1 is in contact with a 
reservoir at voltage V . Notice that in case (a) some current is 
sucked into the junction region from lead 3, while in case (b) 
this effect goes away and the junction behaves as a perfect 
"circulator" , transmitting all incoming current from lead 1 to 
lead 2. The response for the X- fixed point is similar, with 
the current flowing from 1 to 3 instead. 
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FIG. 9; RG flow diagram for 1 < g < 3. For <?!> / 0, tt the 
system is renormalized toward the chiral fixed point x±- Foi^ 
the time-reversal symmetric case (j) — 0,tv, the infrared fixed 
point is the nontrivial M fixed point. It is unstable against 
the addition of flux (f). 



derstood as a consequence of pair tunneling, and this is 
pictorially illustrated in Fig. ^2 

The RG flow diagram for this range of g's is shown in 
Fig. El For </) ^ 0, TT, the system is generically renor- 
malized into the Dp fixed point. Although the Dp fixed 
point itself does not break the time-reversal symmetry, 
the system does not reach there with (p — 0,tt. In the 
time-reversal symmetric case (/) = 0, tt, the property of 
the infrared M fixed point is not known precisely. It 
should be unstable against inclusion of the flux. 



E. 3 < 3 < 9 

There is a stable fixed point Dp , which becomes a spe- 
cial point on the fixed manifold in the limit g —^ 3 as 
discussed above. The conductance tensor at Dp takes 
the symmetric form {Ga = 0) with Gs = (4g/3) (e^//i). 
Including the non-interacting leads, the effective conduc- 
tance tensor is again symmetric {Ga — 0), but the effec- 
tive single-terminal conductance is Gs = (4/3) {e^ /h). 

Notice that the maximum values Gs — (45/3) (e^//i) 
or Gs = (4/3) (e^//i) is larger than the maximum con- 
ductance Gs = (8/9) /h for unitary scattering of single 
particle states; this shows that multiparticle scattering at 
the junction takes place for interacting quantum wires. 
Indeed, the enhanced conductance at Dp may be un- 



F. 9 <5 

The time-reversal invariant system with (p — 0, tt is 
now renormalized into the stable Dn fixed point. For 
(f) ^ 0, TT, the infrared stable fixed point is still the Dp 
fixed point as in the case 3 < g < 9. The fiow diagram 
for g > 9 is shown in Fig. El 

As far as the conductance is concerned, the Dn fixed 
point is identical to the Dp fixed point. The difference 
is in the scaling dimension of the leading irrelevant op- 
erator; it is (7/3 for Dp but is (?/9 for Dn. As both Dn 
and Dp are stable, there must be another critical fixed 
point between these two; its properties however are not 
known. 
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FIG. 10: RG flow diagram for g = 3. For <j!> 7^ 0, tt the system 
is renormalized toward the RG fixed manifold, which includes 
X± as special points. For the time-reversal symmetric case 
(j) = 0, TT, the infrared fixed point is the nontrivial M fixed 
point. It is unstable against the addition of fiux 0. 




FIG. 11: Pictorial representation of the Dp fixed point. The 
conductance is enhanced by the Andreev reflection process. 
We only illustrate the process where the electron pair exits 
through lead 3; but notice that, since the Dp flxed point is 
time-reversal symmetric, the pair has exactly the same am- 
plitude for exiting through lead 2. 




FIG. 12: RG flow diagram for 3 < 3 < 9. For / 0, tt 
the system is renormalized to the Dp fixed point exhibiting 
Andreev refiection. For the time-reversal symmetric case <j) = 
0, TT, the infrared fixed point is the nontrivial M fixed point. 
It is unstable against the addition of fiux 0. 



71/2 



-71/2 



9<g 




III. MODEL 

In this Section, we define the model Hamiltonian pre- 
cisely, for the junction introduced in Fig. ^ 

We study only the simplest model of a Y-junction 
which includes TLL effects. Thus we consider only a 
single channel of electrons in each wire, ignore electron 
spin, assume that the interactions are short-ranged, and 
ignore phonons and impurities. We also assume that all 
three wires are equivalent in the bulk, and consider the 
Z3 symmetric junction. (The generalization to the Z3 



FIG. 13: RG flow diagram for 9 < g. For / 0,7r and 
(j) — 0,n, the system is renormalized to Dp and Dn flxed 
points respectively. Both exhibits the enhanced conductance 
due to the Andreev reflection. 



asymmetric junctions is straightforward, and we find in 
Sec. IX Dl that the Z3 asymmetry turns out to be irrele- 
vant.) 

It is convenient to start by defining a tight-binding 
version of our model. Let ^l)n,i annihilate an electron on 
site n on wire i. Here n = 0, 1, 2, 3, ... 00 and j = 1, 2 or 
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3. The Hamiltonian is: 
where: 

oo 3 



(3.1) 



Ho = -tJ2 T.^^h'^n+i.j + h.c), (3.2) 

n=0 j = l 



,^0/3 ,/,t 



0,i-l 



h.c. 



and 



(3.3) 



(3.4) 



Here n„j = jV'nj is the number density and the wire 
index j = is identified with j = 3. Note that we have 
allowed for a non-zero phase (j) in the tunneling matrix 
element Te^'^/'^ . By making redefinitions of the phases 
of the fields V'n,i, independently on the different wires, 
we can make the phase of each tunneling matrix element 
equal. However, we cannot transform away a uniform 
phase, (/i/3. We may think of the phase as being the 
fiux through the junction. It must vanish if the system 
is time-reversal invariant but will generally be non-zero 
otherwise, for example if a magnetic field is applied to 
the junction. 

While it seems intuitively reasonable that (j) should be- 
come non-zero in the presence of a magnetic field, it is not 
so obvious how (j) actually depends on the field. To study 
this issue it is convenient to consider a slightly more real- 
istic model in which each lead is connected to a circular 
conducting ring, which may contain a flux. In particular, 
if the ring is very small and the Fermi energy of the leads 
is in resonance with a single unoccupied energy level in 
the ring, so that we may ignore all other energy levels of 
the ring, then (p actually vanishes. This follows since the 
tunneling Hamiltonian then takes the form: 



ip\ -d -I- h.c] 



(3.5) 



i=i 



where d annihilates an electron at the single level in the 
ring. We are now free to redefine the phases of ipn,] 
independently on each lead, making all the a^'s real. This 
is nothing but the junction constructed with a "quantum 
dot" as introduced in Ref. 0. 

On the other hand, if two or more energy levels in the 
ring are considered, (j) will generally be non-zero. More 
generally, if the operators di annihilate electrons in the 
ring in cigcnstates with energies i?;, then we may write 
the tunneling Hamiltonian as: 



E E ['^j'' + 



(3.6) 



The tunneling amplitudes, ajj from the lead to the 
^th j.jj^g g^g^^g niay be assumed to be j-independent con- 
stants multiplied by the wave-functions of the ring at 
the locations of the wire, at angle 9 = 27rj/3. We 
may take the ring wave- functions to be fi{9) = e*'^ for 
integer I, with energies Ei — iP /2. A non-zero dimen- 
sionless fiux, (f), changes the energies to /(/ — 0/27r)^/2 so 
that the states with angular quantum numbers ±/ do not 
have the same energy. We may assume that the phases 
of the tunneling amplitudes aj^i are given by the phases 
of the ring states at the locations 9 = 27rj'/3 where the 
jth -yyjj-g joins the ring: 



a,j^i ^ ai e 



(3.7) 



To get a tunneling Hamiltonian of the form of Eq. H3.3() 
we assume that \ai\ << AE, the level spacing in the ring. 
Then we consider processes of second order in the oj's: 1) 
in which an electron tunnels from a lead to an unoccupied 
ring state I and then tunnels off to a different lead, and 
2) in which an electron in an occupied ring state I first 
tunnels off to a lead and then an electron from a different 
lead tunnels into the ring to fill the hole. This gives a 
low-energy effective tunneling Hamiltonian of the form of 
Eq. (EiSl) with: 



Ep — El 



(3.8) 



We see that if only one (occupied or unoccupied) level I 
dominates the sum, then: = 27r/, equivalent to zero. 
If 2 levels dominate the sum, say h and I2, with the 

contribution from the li^^ level much bigger than the 

contribution from the ^2^^, then 



\ai2 



E,, 



E, 



\ai, |2 El, - Ef 



sin[2^(?2-?i)/3], 



(3.9) 



with 101 << 1. So, (/> becomes small near a resonance. 
Of course, when there is zero flux {(j) — 0) in the ring, 
(f) also vanishes due to the cancellation of the ring levels 
with quantum numbers ±L In general, <j) will be some 
periodic odd function of the flux (j) in the ring. 

We now consider the continuum limit of the tight- 
binding model of Eq. (|3.1|l . We introduce left and right 
movers as usual, keeping only narrow bands of wave- 
vectors near the Fermi "surface" at ±fci?: 



ipRj{na) + e 



— ikpna 



^Lj(na). (3.10) 



(a is the lattice spacing.) Linearizing the dispersion re- 
lation. Ho then becomes: 



Ho 



IVp 



J 



dx 



(3.11) 



where vp = 2tsa\kp. The interaction Hamiltonian iJj, 
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takes the form 

POO 

poo 

+ VrlY, JL,jix)JR,j{x) , (3.12) 

where 

Jl,]{x) ^ ipl ji^L,] and Jrj{x) ^ ^pl^ .ipB.,j (3.13) 

and Vrl — 2V [1 — cos{2k pa)]. The first term in 
Eq. (|3.12(l can be absorbed into a renormalization of the 
Fermi velocity. It is the second term, the one propor- 
tional to Vrl, that controls the Luttinger parameter g. 
For repulsive interactions (Vrl > 0), 5 < 1, while for 
attractive interactions {Vrl < 0), g > 1. 

The open boundary conditions in the tight binding 
model are equivalent to a vanishing boundary condition 
at X = —a, i.e. 



e-^fc^-^^(_a) + e'''^>L(-a) = 0. 



(3.14) 



Since ipRj is a function only of (x — vpt) and ipLj of {x + 
vpt), Eq. (|3.14() implies that we may regard the right- 
movers as the analytic continuation of the left-movers to 
the negative x-axis: 



^Rjix) 



2ikf 



"i^Lji-x), (x>0). (3.15) 



Thus we can apparently write Hq entirely in terms of 
left-movers: 

dx i^lj—^Lj- (3.16) 



However, treating the boundary condition of Eq. (|3.14l) 
more carefully, we see that H actually contains an ex- 
tra term when kp 7^ 7r/2a. To obtain this boundary 



condition on the left-movers we must add a term to the 
Hamiltonian: 



(3.17) 



The corresponding Schrodinger equation (which does not 
mix the leads at this point) is: 



iwi^— - r(5(x)*j = (3.18) 



The energy eigenvalues are E ~ vpk, and the wave- 
functions are: 

*j(x) = A„,je~''", (x>0) (3.19) 
= Autj e-*'=^ (x < 0). (3.20) 

Interpreting ^'j(O) in the second term of Eq. H3.18|l as 
[■^j{0+) + *j(0-)]/2, this implies: 



_ iVF - r/2 



(3.21) 



tVF + r/2 

From Eq. (|3.14l) we see that the phase of the wave- 
function advances by e~^'''^° upon scattering so that: 



„2ikf 



and hence: 



ivF — r/2 
ivp + r/2' 



-2vp cot{kpa) 



(3.22) 



(3.23) 



In the continuum limit, the total boundary term in H 



is: 



Hb « ^ {[re*^/3 ^pl^mL.j-m + h.c] + r V1,,(0)Vl,,(0)} , (3.24) 

I 



where the relation between the coupling F in the contin- 
uum limit and the F in the tight-binding model is ob- 
tained in Appendix ^ by matching the scattering matri- 
ces for non-interacting electrons calculated in both cases: 



zero temperature limit), so we henceforth ignore it. It is 
important to note that the continuum form of Hq and Hq 
in Eqs. (|3.16|) and H3.24|) is very general and not restricted 
to the underlying tight-binding model. 



r^\r-2ivp\t/2, (3.25) 

with r given by Eq. H3.23|l . For interacting systems, r 
has no effect on the conductance (in the zero frequency. 
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IV. BOSONIZATION 

The effects of electron-electron interactions on the 
three ID quantum wires is treated using bosonization 
methods. (For a review, see e.g. [s^)- Here we shall fo- 
cus on the simpler case of spinless electrons. In the low- 
energy limit, the effective theory for each quantum wire 
is given by the free boson field theory, which is nothing 
but the Tomonaga-Luttinger liquid. In formulating the 
junction of three wires, we can generalize the "folding" 
trick used in Rcf. 36] and define the three wires on the 
half-line x > 0, so that the boundary at a; = represents 
the junction. The (imaginary-time) low-energy effective 
action of the wires is thus given as 




where j = 1, 2, 3 labels each wire. The electron-electron 
interaction strength is essentially contained in a single 
parameter g, which determines various physical quanti- 
ties and exponents, as we will see later. The cases g < I 
and g > I correspond, respectively, to repulsive and at- 
tractive interactions, while g = 1 is the Fermi liquid. 

The same theory can be also described in terms of the 
dual field 9j for each wire, with the similar action 

3 

J drdxp^^^id^e.f. (4.2) 

The precise definition of the dual field 9j and the deriva- 
tion of the action (|4.2|) is discussed in Sections IVIII and 
IVIIII Here we just note that, while (pj represents the 
quantum mechanical phase of the electron, 6j represents 
the phase of the charge-density wave. In the bulk, both 
phases tpj and 9j are not fixed to a particular value, re- 
flecting the strong fluctuation in one dimension. 

Of course, the effective action (|4.1f) by itself does not 
give an answer to the problem of our interest. In fact, the 
theory defined on the half-line is not completely well de- 
fined without specifying the boundary condition at x — 0. 
In the present formalism, it is the boundary condition 
that represents the physics of the junction. In particu- 
lar, the renormalization-group (RG) fixed points of the 
problem correspond to conformally invariant boundary 



conditions. It turns out that there is a rich variety of RG 
fixed points (conformally invariant boundary conditions) 
for the present problem, as we presented in Sec. ^ and 
will discuss in detail in Sec. |3 

In this Section, we begin our analysis starting from the 
simplest RG fixed point of the problem: the limit of three 
disconnected wires without any transfer of electrons be- 
tween different wires. It is clear that this limit should 
exist as a RG fixed point. The corresponding boundary 
condition is known from the previous studies on the open 
end of a single quantum wire: each boson field 4>j obeys 
the Neumann boundary condition, d4>j/dx = at the 
boundary a; = 0. In terms of the dual field 6j, it is equiv- 
alent to the Dirichlet boundary condition Oj — const, at 
a; = 0. Thus, the disconnected limit can be described by 
the effective action H4.1(l with the Neumann boundary 
condition on ipj imposed. 

We defer the technical discussion of the boundary con- 
dition to Sections IVIII and IVIIII Instead, now we in- 
troduce the electron hopping term H3.3|l as a perturba- 
tion. In the low-energy effective theory, the perturbation 
should also be bosonized. 

The electron annihilation operators at the open ends 
{x — 0) of the three ID wires can be written in terms of 
bosonic fields ipj, j = 1,2,3 [lISllll: 

V'j ~ Vj e'f'l^, (4.3) 

where the r\j are the so-called Klein factors satisfying 
{'7j,?yfc} = 2(5jfe, which are necessary to ensure the an- 
ticommutation relations of the fermion operators in dif- 
ferent wires. Notice that for electron operators on the 
same wire j, the correct anticommutation relations fol- 
low, via the bosonization scheme, from the commutation 
relations of the boson field <pj, and the Klein factor r\j 
plays no role. However, for i ^ j, [(/^i, v'j] = 0; hence, for 
electron operators on different leads, it is only the Klein 
factors that are responsible for the correct anticommuta- 
tion relations. The Klein factors may be represented by 
the Pauli matrices. An alternative representation using 
the boson zero modes is described in Appendix IeI when 
we discuss tunneling in junctions of fractional quantum 
Hall liquids. 

The boundary action Eq. H3.3|l . when rewritten in 
terms of the boson fields, is thus equivalent to 



(4.4) 



It is useful to define the rotated basis 



$2 = -^(i,5l + (/32 - 2(^3) 



and similarly 



(4,5) 



1 



61 

62 = ^(01+02-20, 



2 + ^3) 
2) 



(4.6) 
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The a; = tunneling operators between the three wires 
are more easily expressible in terms of these rotated 
fields. Also, the product of two Klein factors is simplified 
by using Tjjrjj^i — —irjj^i (recall that lead = 3), we can 
write 



where 



and 



iK„<i 



$ = ($1,$2), 

Ki =(-l/2,V3/2) 
i?2 = (-l/2,-V3/2) 
Ks -(1,0). 



h.c. ) + rVS d^^o 



(4.7) 



(4.8) 



(4.9) 



The boundary interactions lead to renormalization to 
an infrared fixed point with a different boundary con- 
dition from Neumann on the two fields $1,2 • However, 
current conservation at the junction requires the "cen- 
ter of mass" field $0 to always obey Neumann boundary 
conditions. Thus the remaining degrees of freedom at 
the boundary comprise the two component boson field 
<f> ~ ($i,$2)- Notice that the term proportional to r in 
Eq. H4.7|l is a constant due to the boundary condition on 
$0, and hence it can be dropped. Hence, the boundary 
action can finally be expressed as 



3 



h.c. 



(4.10) 



The scaling dimension of the hopping term in the dis- 
connected limit (r = 0) is calculated by standard meth- 
ods as 1/g, which is not affected by the Klein factors 
nor the magnetic fiux. The hopping term is thus irrele- 
vant for g < 1, and the N fixed point corresponding to 
disconnected wires is stable. 

When g > 1, the hopping term is relevant and the N 
fixed point is unstable. The system is renormalized into 
an infrared fixed point different from the N fixed point. 
The infrared fixed point should correspond to some con- 
formally invariant boundary condition. Identification of 
such a fixed point in general is a rather difficult prob- 
lem, even in the pure bosonicproblem without the Klein 
factors as discussed in Refs. [l|,[la,li6, 33]. 

In the present case, moreover, additional complicacy 
arises from the Klein factors, namely due to the Fermi 
statistics of the electrons. In the following Sections, us- 
ing various approaches, we study possible infrared fixed 
points that are reached when the effect of the Klein fac- 
tors is included. 



V. STRONG HOPPING LIMIT 

As we have discussed in the previous section, for g > 1 
the hopping term is a relevant perturbation to the N 



(disconnected) fixed point. We would like to know where 
the system flows into, once the relevant hopping term 
is added to the N fixed point. Tracing the RG fiow ex- 
actly, starting from the N fixed point down to the infrared 
fixed point is generally difficult, except in fortunate cases 
where the flow is integrable. Even when we do not have 
the exact solution of the integrable RG flow, it often hap- 
pens that some of the flxed points can be given exactly. 
Lacking the solution of the RG flow itself, we generally 
resort to conjecturing the RG flow, based on the known 
flxed points. 

Thus the flrst step to construct the phase diagram is 
to construct the candidates for the infrared flxed points. 
Since the small hopping is relevant and grows under the 
RG transformation, a naive guess is that the hopping 
strength "goes to inflnity" in the infrared limit. This 
simple guess turns out to be rather useful in many cases. 
However, as it will become clear later in this paper, the 
limit of the "inflnite hopping strength" actually depends 
on the representation of the given problem. 

A natural idea may be to take the hopping ampli- 
tude r to inflnity in the original electron representation 
Eq. (|3.1|) . However, in the presence of interaction, it is 
difflcult to know what happens in the T — > 00. Instead, 
we could discuss the strong tunneling limit by taking the 
hopping strength F to inflnity, in the bosonized represen- 
tation Eq. H4.1L)|I . This is generally different from taking 
the hopping amplitude to inflnity in the original electron 
models. 

This approach was taken in Ref. ^ for the resonant dot 
model of the junction. Here we apply a similar argument 
to our model, namely the junction with direct hoppings 
between the wires, or equivalently the "ring" junction. 
The hopping term in the bosonized representation corre- 
sponds to the boundary potential for $. Thus, the fleld $ 
would be pinned to the potential minimum at the bound- 
ary in the limit F — )■ 00. Namely, the strong hopping 
limit corresponds to the Dirichlet boundary condition on 
$. However, in fact, there are different types of Dirichlet 
boundary conditions (or more precisely, boundary states) 
with different stability. All the Dirichlet boundary con- 
ditions, on the other hand, give rise to the same conduc- 
tance tensor. As was discussed by Nayak et al. HI, the 
Dirichlet boundary conditions exhibits an enhanced con- 
ductivity which may be understood as a result of Andreev 
reflection (see Ref. [25|). 

In order to analyze the stability of the Dirichlet bound- 
ary condition, we have to identify the potential minima 
of the hopping action Eq. H4.10|l at the boundary. As 
was emphasized by Nayak et al. Q , the stability of the 
Dirichlet boundary condition is indeed affected by the 
presence of the Klein factors. 

To demonstrate the importance of the Klein factors, 
let us first consider the fictitious problem without them. 
If there were no Klein factors in Eq. (I4.4|l to begin with. 
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the boundary action would read 
f 3 

S'B = -2r / dT^cos(i?a-$ + 0/3). (5.1) 

a=l 

This is equivalent to the model studied by Yi and 
Kane O (see also Ref. [s^) in the context of the quan- 
tum Brownian motion on a triangular/honeycomb lat- 
tice. The minima of this boundary potential is given by 
a honeycomb lattice with lattice constant (distance be- 
tween nearest neighbors) 47r/3 if 0/3 = 7r(mod27r), and a 
triangular lattice with lattice constant 47r/\/3 otherwise. 

The stability of the F — s- cx) fixed point is conveniently 
studied after integrating out the fields <i>(T, x) for a; > 0, 
thus deriving an effective boundary action. [A more de- 
tailed discussion in a related context is given in Subsub- 
section (|X C 2|l .] Consider a general boundary field ^(r), 
in the imaginary time formalism at temperature 1 = 1//3, 
with Fourier expansion: 

^b{r)^-Y.H^n)e'^--. (5.2) 

n 

(Here LOn = 2?™//?.) The solution of the classical equa- 
tions of motion, 9^$(r, x) = with boundary condition 
$(t, 0) = 3'b(T), (defined for a; > only) is 

$''\t, x) = ^J2 $„e*""^-l'^"l^. (5.3) 

n 

Substituting into Eq. (|4.1|) gives the non-interaction term 
in the boundary action (at /3 —> oo): 

/d 
^|w|$*(w) •$(.;). (5.4) 

The full boundary action also includes the potential term 
of Eq. (|5.1|l . which can be written entirely in terms of 
the boundary field <&b(T), the Fourier transform of ^{oj). 
When r — > oo the path integral is dominated by con- 
figurations where the fields $h(T) stay close to one of 
the minima of the potential. Rare tunneling events (in- 
stantons) occur in which ^^(t) goes from one minimum 
to a neighboring one. In general, multi-instanton con- 
figurations must be considered in the dilute gas approx- 
imation. This is reviewed in some detail in a related 
context in Subsec. IjVI C|l . The conclusion is that, for a 
multi-instanton configuration in which $ tunnels between 
minima of V separated by Mi at time t^, the instanton 
interaction term in the action, for large time separations, 
|ri — Tjl, takes the form: 

This implies that a single tunneling process corresponds 
to an operator in the effective Hamiltonian, at the infinite 
F fixed point, with a RG scaling dimension: 



This is 4(7/9 and 4g/3 respectively for the honeycomb and 
triangular cases. Thus these Dirichlet boundary condi- 
tions are stable for g > 9/4 and 5 > 3/4 respectively. 
This would be troublesome especially in the noninteract- 
ing case g — 1, a,s the "triangular" Dirichlet boundary 
condition is stable. Because the Dirichlet boundary con- 
dition corresponds to the enhanced conductance due to 
the Andreev reflection, it should not be realized if there 
is no interaction. 

However, the correct bosonized representation for the 
junction of electron systems should include the Klein fac- 
tors (as Pauli matrices) as in Eq. (|4.1(J|) . to reproduce 
the effect of the Fermi statistics of electrons. Following 
Ref. 01, the classical value for the boundary field $ is 
obtained by diagonalizing the boundary term in the aux- 
iliary two-dimensional Hilbert space on which the Pauli 
matrices act. The resulting potential is 




with two branches corresponding to ± for a given bound- 
ary value of $. The potential minima of $, in the 
—branch, are those which make the argument of the 
square-root maximum. As in the case without the Klein 
factors, the structure of the potential minima depends 
on the flux 0. For the time-reversal symmetric case 
(/) = 0, TT, the potential minima form a honeycomb lat- 
tice with the nearest neighbor distance 27r/3. For other 
cases (0 ^ 0, tt) the potential minima is instead a trian- 
gular lattice with the nearest neighbor distance 27r/-\/3. 
The scaling dimension of the leading perturbation, which 
is proportional to the nearest neighbor distance squared, 
is (7/9 and g/3 respectively for the honeycomb and trian- 
gular cases. 

As was mentioned above, the "strong hopping" limit 
in the bosonized representation corresponds to Dirichlet 
boundary condition for the boson field $ in either case. 
However, the operator content (and thus the stability) 
differs in these cases. Thus we distinguish them by re- 
ferring to the honeycomb (</> = 0, tt) and the triangular 
(0 7^ 0, tt) cases as Dn and Dp fixed points respectively. 
N of Dn stands for Nayak et al., as it is identical to the 
one discussed in Ref. [3 for the junction with the quan- 
tum dot. P of Dp stands for pair-tunneling; let us now 
show that the Dp fixed point may be regarded as the 
limit of strong tunneling of electron pairs. 

We emphasize again that the strong hopping limit 
r — *■ cxD in the bosonized representation is different from 
the F 00 limit in the original electron model. This 
is evident in the non- interacting case g = I, where the 
original electron model should be always described by a 
free electron scattering matrix. In this case, the strong 
hopping limit in the bosonized representation gives the 
Dp or Dn fixed point, which exhibits the Andreev re- 
fiection and thus cannot be described by free electron 
scattering. For 5 = 1 both Dp and Dn fixed points are 
unstable. Thus, as far as the infrared stable fixed point is 
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concerned, the original electron model and the bosonized 
representation give consistent conclusions. 

If we introduce the hopping of electron pairs instead 
of electrons, the boundary action should be given by 

S'ts = 2Tp drY^ cos (2Ka ■ $ + 20/3^ + h.c. (5.8) 

a=l 

The sign of the pair-tunneling parameter Fp, may be 
deduced from second order perturbation theory in the 
single-electron tunneling parameter T. From the operator 
product expansion: 

for a positive constant, C, and taking into account that 
T]f = —1, we conclude that Fp (x F^ > 0. We note that 
the Klein factors are absent from the action. This is nat- 
ural, considering that the electron pair is a bosonic ob- 
ject. For (/) ^ 0, TT, the potential minima form a triangular 
lattice with the nearest neighbor distance 2tt/^/3. Thus 
the boundary condition in the strong pair-tunneling limit 
Fp — !■ oo is identified with Dp . For i/i = 0, tt, the poten- 
tial minima form a honeycomb lattice with the nearest 
neighbor distance 27r/3. Thus the strong pair-tunneling 
limit is again identified with the more unstable Dn , in 
the time-reversal invariant cases 4> — 0,tt. 

Therefore, the strong pair-tunneling limit Fp ^ oo 
gives the same boundary conditions as the strong hop- 
ping limit T —^ oo. As the pair-tunneling term would 
be generated through the RG transformation even if the 
microscopic model only contains the single electron hop- 
ping, this serves as a consistency check. 

We note that even though the junction with cj) ^ 0,71 
is not time reversal symmetric, the strong hopping limit 
for ^ 0, TT is the Dp fixed point which is time-reversal 
symmetric in itself. This is of course not a contradic- 
tion because the symmetry may be higher at the in- 
frared fixed point than in the original microscopic model. 
This symmetry restoration may be understood naturally 
by considering the maximally time-reversal breaking flux 
(j) ~ ±7r/2 for the electron hopping model. At this value 
of the flux, the pair-tunneling model Eq. (|5.8(l is actually 
time-reversal invariant, as the Aharonov-Bohm phase for 
the electron pair is tt instead of ±7r/2. 

What is more surprising is that the model with (/) = 0, tt 
apparently cannot reach the Dp fixed point, although 
there is no symmetry forbidding it. This was deduced 
from the structure of the potential minima in the Fp — > 
oo limit, as discussed above. We note that this result re- 
lies on the assumption that the sign of the coefficient F p 
of the potential does not change under the RG flow. A 
similar assumption was used in Ref. ,16j for the related 
model of the quantum Brownian motion. Of course, if Fp 
is the only coupling constant in the theory, it should not 
change sign because the RG flow cannot cross the fixed 
point corresponding to the Neumann boundary condi- 
tion. However, in the present problem, the validity of 



the assumption is subtle. This is because one should, 
in principle, consider both the single-electron hopping F 
and the pair hopping. Fp simultaneously in discussing the 
RG flow. 

While we have not yet resolved this question, the above 
arguments lead to the following result on the phase dia- 
gram of the system. 

A. time-reversal symmetric (p = 0,n case 

The N fixed point (decoupled wire) is stable for g < 1 
and unstable for g > 1. The Dn fixed point is unstable 
for g < 9 and stable for g > 9. Therefore, the system 
should flow to some "intermediate" stable fixed point for 
I < g < 9. The intermediate fixed point corresponds to 
neither Dirichlet nor Neumann boundary conditions. 

This conclusion is similar to that in Ref. |^ for a dif- 
ferent setup of the junction. The only difference is that 
the N fixed point in their model is unstable for g > 1/2 
because of the resonant transmission through the dot. 

B. time-reversal breaking 7^ 0,7r case 

In our model, we have the flux </) as a new parameter 
which breaks the time-reversal symmetry, and can lead 
to new physics. 

In fact, as it is seen above, the new Dirichlet fixed 
point Dp can be realized for (p Q,7t. It is stable for 
g > 3, so we expect a generic junction without time- 
reversal symmetry to fiow to Dp for g > 3. We note 
that although the Dp fixed point itself is time-reversal 
invariant, it seems that the time-reversal symmetry in 
the original junction model must be broken to reach Dp 

For 1 < g < 3, as both the N and Dp fixed points 
arc unstable, there must be an intermediate stable fixed 
point to which the generic system flows. However, the 
present argument does not shed light on the nature of 
the stable fixed point. In the following sections, we will 
argue that the stable fixed point for 1 < g < 3 is solvable 
with different techniques that we present in this paper. 

VI. CONNECTION WITH THE DISSIPATIVE 
HOFSTADTER MODEL 

As we have seen in the last section, there are infrared 
stable fixed points that cannot be identified with either 
weak- or strong- hopping limit in the bosonized represen- 
tation Eq. H4.10|l with the Klein factors. 

The presence of the phase factors due to the fermionic 
statistics, encoded in the Klein factors, makes the anal- 
ysis more difficult than that for the standard problems 
with only bosonic terms in the boundary action. The 
extra phases due to the fermionic statistics appear, for 
example, if one attempts a perturbative expansion in the 
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hopping amplitude F. The coefficients appearing in each 
order in perturbation theory are different than those ob- 
tained without the phase factors, and hence they alter 
the nature of the strong hopping limit. 

In this section, we develop an alternative representa- 
tion, in which the fermionic statistics of electrons are en- 
coded without using the Klein factors. Namely, we will 
show that the perturbative series in the hopping ampli- 
tude r for the Y-junction problem is identical to the series 
obtained for the dissipative Hofstadter model (quantum 
motion of a single particle under a magnetic field and 
a periodic potential, subject to dissipation) studied by 
Callan and Freed in Ref. |2^. In the dissipative Hof- 
stadter action, phase factors arise from the motion in 



a magnetic field, but the action contains solely bosonic 
terms, and hence it is amenable to standard methods for 
determining the strong coupling behavior, such as instan- 
ton expansions. 



A. Generalized Coulomb gas from the Y-junction 

We consider the perturbation theory in the hopping 
parameter F, starting from the bosonized representation 
Eq. (|4.10l) with the Klein factor. The partition function 
Zy of the system can be expanded as 



Zy 



E 



dTn \ dTn- 



dri 



giL2-*(T2)gi-Li-*(ri)\ 



(6.1) 



A V 



(a) 



(b) 



FIG. 14: A typical closed path occurring in perturbative ex- 
pansion of partition function 



FIG. 15: (a) a counter-clockwise loop around an "up trian- 
gle" . (b) A counter-clockwise loop around a "down triangle" . 



where Lj = QjKaj is one of the six vectors ±Ki 2.3 {qj = 
±1 and aj = 1,2,3), and ( )o denotes the expectation 
value with the Neumann boundary condition on The 
= e"i^'^/^r]a^ in the [. . .] term in Eq. (EU keep 

track of the phase factors. 

Now, each term in the expansion Eq. (|5.1|l corre- 
sponds to a closed directed path formed by the vectors 
Li, L2, . . . , Ln, on the triangular lattice spanned by the 
primitive vectors Ki and K2. An example of a closed 
path is shown in Fig. A charge vector Lj is asso- 
ciated to each vertex operator e''^^ *, and the condition 
that the sequence of vectors must form a closed path is 
a consequence of the charge neutrality condition of the 
Coulomb gas expansion. 

Let us now argue that the total phase factor for a given 
path C = {Li, L2, ■ ■ ■ , Ln), 



(6.2) 



is the sum of the phase contributions from all the ele- 
mentary triangles enclosed by the loop. The phase con- 
tributions from "up" and "down" triangles, however, are 
different. The phase factor for a counterclockwise loop 



on the up triangle [see part (a) of Fig. I15j is determined 
to be 

e**^^ = (_j,^3)e*^/3(-iry2)e*'^/^(-z77i)e"^/3 = e'^ (6.3) 

while that for the counterclockwise loop on the down tri- 
angle [see part (b) of Fig. ^] is 

e'^v ^ (i,73)e-^'^/3(i772)e-''^/3(^^^)g-#/3 ^ ^^i^-<^>) ^ 

(6.4) 

Namely, loops on the triangular lattice pick up a phase 
as if there is a staggered magnetic flux of <j) and tt — </) in 
each elementary triangle. The total phase accumulated 
due to the loop C is thus 



(6.5) 



where N/^iC) and Nx^{C) are the net numbers of up 
and down triangles enclosed by the path (a triangle con- 
tributes ±1 to JVa and depending on whether it is 
traversed clockwise or counter-clockwise). 

The argument why the phase for a given path C is 
the sum of the phases due to the enclosed elementary 
triangles is constructed inductively. First consider the 
special case where the closed path consists of a single non- 
intersecting closed loop. This implies that each vertex on 
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the loop is visited a single time except for the first vertex 
which is visited twice, at the beginning and end of the 
loop. The phase from this loop is equal to the sum of 
phases of elementary triangular loops into which it can 
be decomposed. Say it is so for all paths enclosing up 
to m triangles, and then consider a loop enclosing m + 
1 triangles; we must argue that the phase accumulated 
around the loop enclosing m + 1 triangles is equal to the 
sum of the phases enclosed by two smaller area paths 
with nil a-iid m2 triangles {mi +m2 = m+l, mi 2 l£ so that 



Consider then an arbitrary non-intersecting loop C = 
{Li, L2, ■ ■ ■ , Ln), which we break into two loops by 
adding an internal segment that splits the loop into 
two loops. Let the segment be the sequence of vectors 



L'i,L'2, . 



,L' . It is easy to show that 



(6.6) 



e = 



Cf!' ■ ' ' Cf' Cf' Cr, ■ ■ ■ Cl^ Cf* 



(6.7) 





+ 

A 



FIG. 17: Decomposition of a path into non-crossing loops 



FIG. 16: Decomposition of a non-crossing loop into elemen- 
tary triangles 



for two loops £i_2 that encircle areas smaller than the 
original loop £ , i.e., li^2 < ™ and with li-\- 12 = m -\- \, 
as required by the induction argument. An example of 
a decomposition of a non-crossing loop into elementary 
triangles is given in Fig. [TCI 

Finally we must argue that an arbitrary closed path 
can be decomposed into a set of non-intersecting loops 
such that the total phase of the path is the sum of 
phases for each non-intersecting loop. This can be seen 
as follows. Consider following an arbitrary oriented path. 
Consider the first time that any vertex on the path is 
visited for a second time. Label this vertex Vi. The 
sequence of edges between the first and second visit to 
Vi defines a non-intersecting closed loop, Li. (This fol- 
lows since the first intersection on the path occurs at Yi 
at the instant of the second visit.) The phase due to 
this closed non-intersecting loop may be calculated un- 
ambiguously and makes an additive contribution to the 
total phase of the path. This follows since the product of 
the Klein factors around a closed loop is always propor- 
tional to the identity matrix. Now excise the closed loop, 
Li, from the path. The remaining path remains closed. 
Consider the first time that any vertex on this excised 
path is visited twice, at some other vertex, V2. This de- 



fines another closed loop L2 which we again excise. We 
continue in this way until we finally revisit the first point 
on the path. This final revisitation defines a final closed 
loop L„. An example of the decomposition of a path into 
non-crossing loops is given in Fig. El Note that when a 
path (or an excised path) makes a U-turn and retraces a 
link this counts as a closed loop of length 2 and zero area. 
The total phase associated with the path is the sum of 
phases of each non- intersecting closed loop. These phases 
can each be expressed as a sum of phases of elementary 
triangles into which they can be decomposed. Adding up 
the net number of up triangles and down triangles in all 
the non-intersecting loops gives us Eq. H6.5|l for the total 
phase. 

Let us now turn to the correlation function entering in 
Eq. EU: 



„»L„-<E.(r„) 



^ giL2-*(r2)giLi-*(Ti)\ 



^KiJ^-jLj) exp 



i^Lj • Lk ln|T,- - Tfel 



j>k 



(6.8) 



where Sk{L) = 1 if L = 0, and Sk{L) = otherwise. The 
factor SK{J2j ^j) enforcing the charge vectors Lj must 
form a closed path, as we mentioned above, explicitly 
follows from of charge conservation in the Coulomb gas 
expansion. 
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Using this expression, in conjunction with the phase 
factor we discussed above, we can rewrite the partition 



function as: 



J 



(6.9) 



r 



Thus, the partition function Zy is written in a simi- 
lar form to the so-caUed Coulomb gas (particles inter- 
acting with the logarithmic potential) in one dimension. 
The difference from the standard Coulomb gas is in the 
phase factor e*^^, which is indeed the consequence of 
the Fermi statistics of the electron and of the Aharonov- 
Bohm phase. We may call this a "generalized Coulomb 
gas." In fact, below we show that this is indeed identical 
to the generalized Coulomb gas which appeared in a quite 
different problem - the dissipative Hofstadter model on 
a triangular lattice. 



Generalized Coulomb gas from the dissipative 
Hofstadter model 



Dissipative quantum mechanics (also known as quan- 
tum Brownian motion) is the problem of the quantum 
motion of a single particle subject to dissipation (see 
Ref. for a review). In the presence of any periodic po- 
tential, single-particle quantum mechanics implies that 
the particle never localizes due to the tunneling effect. 
However, if the dissipation is sufRciently strong, the tun- 
neling can be suppressed and the particle may be local- 
ized. When the dissipation is Ohmic, there is a surprising 
mapping of the problem to the boundary problem of the 
free boson field theory. Here, the bulk of the free boson 
is related to the heat bath in the dissipative quantum 
mechanics, and only the boundary field is related to the 
original particle. By integrating over the bulk degrees of 
freedom, the boundary problem can be reduced to the 



one-dimensional Coulomb gas. 

Dissipative quantum mechanics in the presence of the 
periodic potential and the magnetic field is known as the 
dissipative Hofstadter problem. It is again related to 
the boundary problem of the free boson, or the following 
one-dimensional generalized Coulomb gas. The "free" 
action for the dissipative Hofstadter model, namely the 
effective action in the absence of the potential (but with 
the magnetic field and the dissipation), is 



(6.10) 

where — 1,2 and £12 = —£21 = 1, en = £22 = 0. 
a and P are related to the dissipation and the magnetic 
field, respectively. This determines the propagator: 



D^.ir) - (X^(t)X.(O)) = 
Inr Ouv + nT- 



sgnr e^^. (6.11) 



We now introduce a potential term with the periodicity 
of the triangular lattice (notice the difference from the 
"rectangular" one in Ref. 0), as 



Sv[x] = -ve^i^ dTj2 



e^^-^ + c.c. 



(6.12) 



a=l 



where V and S are chosen to be real. Expanding the 
partition function in powers of V, one obtains 



J 



dTn / dT„_i ••• / dri ^ '^Li 



iL„-X(r„) „iL2-X{T2)„iLi-X{Ti)^ 



0, (6.13) 



where again Lj — qjKa^ is one of the six vectors ±Ki^2,3 correlation function is given by 
(q, = ±1 and a,- = 1,2,3), and the Kf = e"'^^/^ in ^ ^ .... 

r 1 /„ii„-X(T„) iL2-X{r2) iLrX(Ti)\ _ 

the [. . .] term keep track of phase factors due to 6. The \^ • ■ ■ ^ ^ /o ~ 

^liiEj Lj) exp [^^2^^ Lj ■ Lk In \tj - Tk\^ 



(3 

'a2 + /32 



X Lk sgn(T, 



(6.14) 
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(Because of the time ordering, notice that sgn{Tj—Tk) = 1 
for all terms in the sum.) 

Notice that the K£ in the [. . .] of Eq. H6.13|l commute 
with one another, in contrast to the Klein factors in 
the [. . .] of Eq. (EUJ. However, the ( )o term in Eq. 
contains phase factors not present in Eq. H6.8|l . While the 
phases entering the perturbative expansions for Zy and 
Zbhm are arising from different sources, their net effect 
turns out to be identical, as we show below. 

The phases in the expansion of the dissipative Hofs- 
tadter model can also be associated with the closed paths 
C = (Li, L2, . . . , Ln). The contribution from the Lj x 
terms amounts to a flux through the oriented area Ac en- 
closed by £. To see this we again decompose any closed 
path into a set of closed non-intersecting loops. The ori- 
ented area of such a non-intersecting loop is 



1 V- - 



X Lk 



(6.15) 



for this is carried out by induction, identically to what 
we have used for the products of the (r above. 

The phase contributions from "up" and "down" tri- 
angles are again different. The phase factor for a coun- 
terclockwise loop on the up triangle is determined deter- 
mined to be 



gii5/3 ^i5/3 ^i6/3 _ ^iS 



(6.19) 



while that for the counterclockwise loop on the down tri- 
angle is 



g-i(5/3 _ ^-iS 



(6.20) 



j>k 



Namely, loops on the triangular lattice pick up a phase 
as if there is a staggered magnetic flux of 5 and —6 in 
each elementary triangle. The total phase accumulated 
due to the loop C is thus 



Thus, the accumulated phase A£ due to the magnetic 
flux is 



Ar = 27r 



2-K 



(6.16) 



A[N^{C) + N^{£)] , (6.17) 



where A = ^ is the area of an elementary triangle, and 
7Va(>C), N^{£) are the net numbers of enclosed elemen- 
tary up and down triangles, as before. 

Now we turn to the phase due to the [. . . ] in Eq. 16.13|l . 
The total phase factor for a given closed path £ = 

{Ll,L2, . . . , Ln), 



(6.18) 



Ec^N^i£)6-N^{£)S , (6.21) 
where Ni^{£) and N^{£) are the net numbers of up and 
down triangles enclosed by the loop. 

The total phase factor due to the loop £ is the sum 



where 



Tc = {£) + (-C) , (6.22) 



■Da = 27r 



tA + 6 



9 ^ 
^ + P^ 



A-5 



(6.23) 
(6.24) 



is the sum of the phase contributions from all the ele- 
mentary triangles enclosed by the path. The argument 



We can now rewrite the partition function for the dis- 
sipative Hofstadter model on the triangular lattice as: 



C={Li 



Following Callan and Freed [23j , we may regard the DHM 
partition function as the partition function of a general- 
ized classical 1-dimensional Coulomb gas. Tj is the posi- 
tion of the j*'' particle and Lj is its generalized "charge" , 
a vector quantity taking on six possible values. 



C. Mapping of the Y-junction to the dissipative 
Hofstadter model 

The expression for Zdhm in Eq- 16.25|1 is identical to 
that for Zy in Eq. H6.9|l , provided the following conditions 
are satisfied. First we need the relation 
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for the scaling dimension of the leading perturbation to 
be matched between the two theories. Furthermore, in 
order to match the phase factor, we require Tc = Tc 
mod 2tt, or equivalently 



mod 2n, 
mod 27r 



(6.27) 
(6.28) 



number of different choices of a and (3 labeled by the 
integer n. Eqs. H6.26|l and l|6.3()(l give: 



a 



3.g 



3 + (2n - l)2g2 
V3(2n - l).g2 



(2n 



1)V 



(6.31) 



The above two equations imply: 



73- 



mod 27r, 



mod 27r. (6.29) 



5 



(2n- 1) 

+ 7r(l/2-n). 



(6.30) 



for an arbitrary integer, n. Because the phases, and S, 
are only defined modulo 2tt, there is actually an infinite 



To understand the ambiguity in the integer n, let us 
define the Coulomb gas charge density 



JV 

p(r)^^L,<5(r- 



(6.32) 



The generalized "energy" , E, of an arbitrary classical gas 
configuration is determined by the iV positions, Tj and 
charges, Lj. This generalized "energy" is now a com- 
plex quantity and exp{—E/T) (where T is the classical 
temperature) is given by Eq. H6.13|l . i.e. 



J 



-| = ^[lny + z5g,/3] 



T 



a^+f3' 



j>fc 



j>k 



r 



(6.33) 



where is defined by Lj = qjKa^ as before. Since x 
Lk is always or ±-\/3/2, exp[— iJ/T] is invariant under 
a shift: 

7rV3 /3 tta/S 



/3 



27r. 



(6.34) 



Eqs. 
under 



/32 2 a2 4-/32 
imply that there is also a symmetry 



7r\/3 /3 



2 a2 

(5 — TT, 



which corresponds to shifting the integer n in Eq. (|6.31|l 
by 1. Since the energy of any configuration is inde- 
pendent of 71, it follows that the Coulomb gas charge 
density correlation function, (p^(r)pi^(r')), is also n- 
independent. 

On the other hand, the correlation function 
{X^^{t)X''{t')), does depend on n. These two correla- 
tion functions can be simply related to each other by 
adding an external source, a(T), to the action, coupling 
to the Coulomb gas charge density. That is, we modify 
the energy of the Coulomb gas by: 



~E/T -E/T -i J dr a(r) • p{t). (6.36) 
The Coulomb gas correlation function is given by the 



functional derivative of the partition function with re- 
spect to the source term: 



92 InZ 



1 



;{p^,{-UJn)p,A-^n))- (6-37) 



Note that the Coulomb gas energy arises from the func- 
tional integral over X: 



-E/T 



1 



(6.35) where: 



y"e*(''/3)i:.9. [dX{T)]e-^'''^^\ (6.38) 



(6.39) 



5, = So{X)^i drXir) ■ p{t 



Here Zq is the functional integral with Sq replaced by Sq, 
which is given by Eq. H6.10|l . Thus, adding a source term 
to the Coulomb gas corresponds to a modification of the 
potential term in the DHM action: 



Sv[X] 



a=l 



(6.40) 

Alternatively, we may shift X{t) X{t) — d{T), remov- 
ing the source term from the potential term in the DHM 
action but inserting it into the non-interaction term, Sq: 
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^0-^1 - a^]*K) {D-')^^ K) [X^'K) - an(c.„). (6.41) 

I 

Here we have introduced the inverse of the DHM propa- Taking a second derivative with respect to we see that: 
gator, defined in Eq. (|6.11l) . which is simply the matrix 
appearing in S'o, in Eq. (|6.1(J|) : 

{D-^)^^ = a5,,M + /?e^,w. (6.42) 
I 

{p,{-c.)pA^c.')) = Siu; + J) {D)-l i-u;) - ^ {D-') {-u:){X^ {-u:)X\~u')) (D)-'^ (u;'), (6.43) 
or, equivalently: 

(X^(-c.„)X-(-^;)) - (2^)25(c.„ +^;)i^^-(-c.„) - {27rfD^^^~LU^){p„{-co,,)p^{-u;'„))D^''{Lu'J (6.44) 

I 



Due to the comphcated dependence of a and f3 on n, 
given by Eq. (|6.31|l . and hence the comphcated depen- 
dence of D'^'^ on n, we see that {X'^X''') is n-dependent. 

However, what we are interested in here, for obtain- 
ing properties of the Y-junction, is only the Coulomb gas 
correlation functions, which are n-independent. Hence 
we may choose any value of n that we find convenient. 
It will sometimes be useful to deduce the long time be- 
havior of {X'^X'^) by RG and physical arguments and 
then use Eq. H6.43|l to deduce the Coulomb gas correla- 
tion function. In so doing we may choose any value of n. 
In general, the behavior of (X'^X^) will depend strongly 
on n but this n-dependence will cancel when we compute 
{p^iP^} using Eq. H6.43|l . 

Let us now use this mapping between the Y-junction 
problem and the dissipative Hofstadter model and con- 
sider the case of g > 1. Then, as we have already dis- 
cussed, the electron hopping is a relevant perturbation 
in the disconnected limit. In the DHM, for any choice 
of n, y is a relevant perturbation for g > I. We would 
like to find the infrared stable fixed point reached in the 
low energy limit. A simple guess is that it occurs at 
V oo. The stability of the ^ ^ oo fixed point can 
be determined using the instanton method p3 | . Namely, 
in the strong potential limit, the X field is pinned at 
one of the minima of the potential (|6.12|) . The leading 
perturbation to this limit is given by a tunneling pro- 
cess between the neighboring minima, represented by an 
"instanton." This is a classical solution, Xci{t) which in- 
terpolates between two adjacent minima of the potential 
as T goes from — oo to oo. The nature of these solutions 
is discussed in Callan and Fried and, in more detail, in 
Ref. "2^ . To determine the scaling dimension of the tun- 
neling operator we must consider the classical action of 



a dilute multi-instanton configuration which tunnels suc- 
cessively from minimum to minimum with each tunneling 
event well separated in time from the rest. This instan- 
ton gas is dilute because, in the limit of large V, the 
action of an instanton goes to infinity, so the density of 
instantons goes to zero. In the dilute limit, we can ap- 
proximate such a multi-instanton by a sum of the form: 



X(r)«Xo + ^X,(r-T,), (6.45) 

where 

X,{t) ^ 0, (r^-oo) 

^ Mj, (r^oo). (6.46) 

Here the j^^ tunneling event takes place at time tj and 
involves the particle tunneling between two adjacent min- 
ima of V{X), displaced by Mj. In the dilute limit each 
function Xj{T) is a single instanton solution. In the long 
time low energy limit of interest to us, we may regard 
the single instanton solutions as having the form: 

XaAr)^Mjf{r), (6.47) 

for a real function, /(r) which interpolates between and 
1. The Fourier transform of the multi-instanton solution 
can then be approximated as: 

X(a>)«^A/,e''^^^/(c.). (6.48) 
3 

The low frequency behavior of /(w), which follows from 
the long time asymptotic behavior of the interpolating 
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function /(r), is given by: 



/H - - 

lUJ 



(6.49) 



The important interaction term in the instanton action 



Er dm 



comes simply from Sq and is obtained by substituting 
Eq. into Sq in Eq. This gives the instanton 

interaction term in the action: 



(6.50) 



In the dilute limit, where the time separations are large, 
we may use the asymptotic form of f{uj) in Eq. I|6.49(l . 



In this we recover the Coulomb gas energy: 



--^^M..A/,ln|r. 

^ ' i>3 



r 



ITT/? X - 
(2^^ 



Mi X Mj sgn(Ti - Tj). 



(6.51) 



We see that this expansion is equivalent to the weak V 
expansion and therefore that the tunneling events are 
relevant only if 



\M,\'a'/{2T,Y < 1. 



(6.52) 



Thus we see that the separation of the minima of the 
potential, \Mi\ is a crucial quantity in determining the 
stability of the localized fixed point. This separation is 
determined by a simple classical analysis of the potential, 
V{X) in Eq. H6.12|l . The stationary points of V{X) occur 
at the points: 



^ z=l 



(6.53) 



independent of 5 and forming a triangular lattice. How- 
ever, the energies (values of V) at these stationary points 
depend on 5. For general values of 5, these energies are 
different at the 3 sub-lattices into which the triangular 
lattice can be partitioned: 



V{Q) 



-6Fcos(5/3) 



v[A-KKi/i^ = -6V^cos[(27r- 5)/3] 
F(-47ri?2/3) = -6Fcos[(27r + ^)/3]. (6.54) 

The sub-lattice containing 0, is a triangular lattice with 
basis vectors, 2'KRi {i — 1,2) with 



i?, = (2/V3) X z. 



(6.55) 



Of course, the other sub-lattices are simply displaced 
from this one by the fixed vectors, {ATT/'i)Ki and 
(47r/3)i?2- Therefore, the vectors Mi occurring in the in- 
stanton solutions are ±27r/ii, ±27ri?2 and ±27r(J?i -I- R2) 
with length 47r/'\/3. As we vary 5 the relative energies 
on the three sub-lattices change. For 5 — Q, the potential 
minima occur on the X = sublattice. As we increase 
5 from zero, the minima remains on the X — Q sublat- 
tice until we reach 5 — tt. For t: < 5 < 27r, the minima 
are rather on the X = (47r/3)^i sub-lattice. Right at 
5 = 7:, there are degenerate minima on both X = Q and 
X = (47r/3)^i sub-lattices. These two sub-lattices define 
a honeycomb lattice with lattice spacing 47r/3, smaller by 
a factor of 1/ a/S compared to that of the triangular lat- 
tices. 

Thus, from the dilute instanton analysis, we conclude 
that, for general (5, the localized phase is stable when 



A = 



a\2TiRi 



4a 

T 



45 



3 + (2n 



1)V 



> 1. 



(6.56) 



While this does not prove rigorously that the DHM is 
in a localized phase for small V when this inequality is 
satisfied, it is reasonable to believe that the system is in 
a localized phase for large enough V , in this region of a. 

For n 7^ or 1, the localized phase is never stable, for 
any value of g. This follows since the maximum value of 
A, with respect to g, for fixed n is 2/[x/3|2n- 1|], which 
is < 1 for n ^ or 1. On the other hand, for n = 
or 1, A < 1, for 1 < 5 < 3. Thus is it reasonable to 
assume that the behavior of the Y-junction, at least for 
large F, corresponds to the localized phase of the DHM 
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in the n = 0, or 1 representation. This result appears at 
first pecuhar because we might just as well have chosen 
a different value of n, in which case we would conclude 
that the DHM is not in the localized phase. For g > 1, 
the DHM is not in the freely diffusing phase cither, since 

V is relevant. It therefore must be in some intermedi- 
ate phase which is neither freely diffusing nor localized. 
Remarkably, Eq. Ht).44|l allows us to determine {XX) 
in these intermediate phases by determining (pp) using 
71 = or 1 H. 

However, an ambiguity still remains because both n = 
and 1 give stable localized fixed points. However, the 
localized fixed points for these two values of n predict 
different (pp). Consistency requires that, while both of 
these fixed points is stable at large V, the RG flow goes 
all the way to the localized fixed point, starting from a 
small V, for only one of these two values of n. For the 
other value of n the RG flow should go to an intermediate 

V flxed point, which must give the same (pp). 

The particular cases 4> = ±7r/2, which maximally 
break time-reversal, are instructive in understanding the 
resolution of this ambiguity. From Eq. Ht).30|l . for 
(j) = 7r/2, the choice n ~ 1 {f3 positive) gives 5 — 0, 
and the potential minima forms a triangular lattice as 
usual. However the other choice n — {P negative) gives 
6 — TT, for which the potential minima form a honeycomb 
lattice with spacing 47r/3, making the strong potential 
limit unstable for all g. Similarly, for (j) = n — 

{(3 negative) is the unique choice that gives a stable flxed 
point. This suggests that these localized flxed points x± 
reflect the breaking of time reversal symmetry due to the 
magnetic flux (j). Indeed, the conductance at these flxed 
points exhibits a chiral behavior breaking the time rever- 
sal invariance, as we show below. 

Let us next consider the case when (jj is very close to, 
but not exactly equal to, 7r/2. If one chooses n = 1, S will 
be close to zero, and the potential minima forms a tri- 
angular lattice with a large separation between the true 
minima and secondary minima. However, if one chooses 
n — instead of n = 1, the minima of the potential 
V{X) will form a triangular lattice with spacing An/y/S, 
but there are local minima with only slightly larger en- 
ergy such that the distance between a true minimum and 
the nearest local minima is 47r/3. We could then envision 
another type of approximate instanton solution where the 
particle tunnels from a true minimum to a local minimum 
where it lingers for a long but flnite time but must even- 
tually tunnel to a true minimum. The typical time that 
the particle spends at a minimum is determined by the 
instanton fugacity, and grows exponentially as ^ oo. 
Thus this other type of instanton might never be impor- 
tant for sufficiently large V. However, for intermediate 
values of V, these instantons could play an important 
role, effectively preventing the RG flow to go to inflnite 

V for the n = choice. This argument suggests that the 
choice of n = 1 should represent the stable flxed point 
around a vicinity of = 7r/2. Of course, similar state- 
ments hold regarding the case (j) very close to — 7r/2, for 



which n = should represent the stable flxed point. 

We are thus led to conjecture that the localized flxed 
point for n — 1 {n — 0) corresponds to the stable flxed 
point of the Y-junction, for 1 < g < 3, even for small 
bare V, for a range of around 7r/2 (— 7r/2). While 
none of our arguments have proven rigorously that the 
RG flow really goes all the way from F = to ^ = oo 
(i.e. to the localized phase) for any values of g and (j), 
this seems like the most "economical" assumption since 
otherwise there would necessarily be intermediate V flxed 
points for all values of n. On the other hand, ior cj) — 
or TT, the n — and n — 1 flxed points appear equally 
stable. Since, as we show in the next sub-section, the 
conductance exhibits broken time-reversal symmetry for 
these localized flxed points, we expect that neither of 
them is stable at </> = 0. There must be some other stable 
fixed point, corresponding to intermediate V, in this case. 
We refer to this fixed point as M (for mysterious). 

It is convenient to now think about the RG flow in 
the space of the physical parameters of the Y-junction, 
r and (j). We have proposed (unique) stable flxed points 
for (j) near tt/2 (x+) or —tt/2 (x-) and another flxed 
point, M, stable for (/> = 0. While the x± flxed points 
are stable against a small change in the flux, we do not 
know whether the M flxed point has such a stability. The 
simplest assumption is that it is does not. In this case 
the RG flows might go to the x± fixed points for any non- 
zero (j), starting from arbitrarily small V. Alternatively, 
it is possible that the M fixed point is stable against 
adding a small flux. In that case, there would have to be 
additional unstable fixed points defining the boundaries 
between the basins of attraction of the M, x+ and X- 
fixed points. So again, an "economy" principle, suggests 
the simple picture with only 3 stable fixed points, for 
1 < g < 3, as sketched in Fig. El 



D. Calculation of the conductance tensor 

Let us calculate the conductance tensor Gjk as defined 
in Eq. H2.2|) at the chiral flxed points x±- There are 
two convenient ways of introducing the voltages drops. 
They can either occur right at the junction, or else can 
be spread over flnite regions of the wires which are long 
compared to microscopic length scales like kp ^ but could 
still be short compared to the length of the wires and 
other long length scales such as an inelastic scattering 
length. Both lead to the same results. The former ap- 
proach has the advantage that it allows the conductance 
tensor to be expressed in terms of the Coulomb gas cor- 
relation function. The latter approach is useful when the 
BCFT formalism is used and is the only method at our 
disposal in the case of the Dirichlct flxed point. We dis- 
cuss each in turn, one in this section and the other in the 
next. For convenience, we set the electron charge, e, to 
— 1. (In our convention the actual electron charge e < 0.) 
We also set ?i = 1. 

We flrst consider the case where the voltage drop oc- 
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curs right at the junction. We then modify the tunnehng that: 
Hamihonian by time-dependent phases, Aj{t): 

Ht -re"^/^^ VlV-j-ie'^^^"-^^-^^ + h.c. (6.57) 



dAj/dt = -Vj. (6.59) 



It is convenient to define Ij as the current on wire i, 
headed towards the junction: 

Note that, in the case where Aj{t) — —Vjt, that this 

corresponds to the replacement: Ij ~ —dNj/dt = i[Nj, H]. (6.60) 

^^.(^t) — > e^^^^V^l^)- (6.58) Here Nj is the total number of particles on wire j. Since 

only the tunneling terms destroy conservation of particle 
This corresponds to a potential Vj on wire j, so we see number on each wire separately, we readily obtain: 



+ h.c. (6.61) 



Expanding to first order in the Aj^s gives: 

- {1°) ' [{A, - A,_i) + {A, - A,+i)] Et/3. (6.62) 

Here is the tunneling energy: 

Et = -re"^/3^(V'jV'j-i) + c.c. = -3re''^/3(V'jV'j-i) + c.c. (6.63) 



J 

^th 



where the symmetry was used in the last step. is the current operator to zero*'' order in the Ak 



+ h.c. (6.64) 



(Ij) is calculated to first order in Aj{t) using: 

H^Ho + H' = Ho+^I^ Aj . (6.65) 



(/°(^))=-^ / di'([/°(t),i7'(0])=-*E / dt' {[I^^(t)jUmAu{t')^Y. / dt' Kf^{t-t')Au{f). (6.66) 

I 

Here K^j* is the retarded Green's function of the /j* op- and the w-independent symmetric matrix, Cjk is given 
erators. Fourier transforming, by: 



Aj{uj) = / dte"^' Aj(t), (6.67) 



Cjj = 2Ej^/3, 

Cjk = -St/3 (jy^fc). (6.70) 
Note that 

and including the "diamagnetic term" in Eq. H6.62|l . we 

obtain the Kubo formula in the form: Cjk — Cjk — , (6-71) 



(^X^) - '^ok{^) ^fc(^)- (6-6^) which combined with Eqs. and ^^oW^ implies 

Y^Kf^ = Y.^Tk-^- (6-72) 



fc 

where the conductance matrix is given by 



1 The retarded Green's function may be obtained by ana- 

Cjk{^) — -T^ \J^jk il-^) — Cjk\ , (6.69) lytic continuation from the Matsubara Green's function: 
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dre^-"^ {TI"{t)IUO)). (6.73) 



Here ujn = 2-KnT for integer n, at temperature T, r is 
imaginary time and T represents time ordering. The 
retarded Green's fmiction is obtained, as usual, by: 



(6.74) 



The conductivity can be conveniently derived from the 
imaginary time path integral, with a vector potential, 
Aj{T), defined as in Eq. (|6.57|l which depends on 
imaginary time. The imaginary time partition func- 
tion, Zi[Aj{T)]^ is obtained from the path integral of 
exp[— 5*0 — J dr Ht{t)], where 5*0 is the imaginary time 
action for the three independent leads. We see that: 



SlnZi 



(6.75) 



and, at Aj 0, 

(52 In Zj 
SAj{t) Mfc(r') 

Fourier transforming: 

6^ InZi 
SAj{ojn)SAk{uj'J 2tt 

(6.77) 

The conductivity can be obtained from this by analytic 
continuation. Explicitly, writing: 



KfUr-T') + 5{T-T')C,,. (6.76) 



In Zj 



1 



the dc conductance is given by: 



Gjk = lim 



1 



Zjk{^n)- 



(6.79) 



Here the limit w„ — > must be taken by first continuing 
LOn into the complex plane and then taking a;„ — > along 
the imaginary axis with a small positive real part. This 
version of the Kubo formula is convenient because this 
second derivative of In Zj is proportional to the Coulomb 



gas density correlation function in the DHM. To make 
this connection it is convenient to rewrite the three vector 
potentials, Aj{T) in terms of a pair of vector potentials, 
a^(T) such that: 



Aj - Afc = ^ tjki Ki 



(6.80) 



Here ejki is the antisymmetric unit tensor (ei23 = 1 etc.) 
and the Greek indices run over 1 and 2 only. Then the 
tunneling term is modified to: 



(6.81) 



The corresponding modification of the DHM action is the 
shift: 



X ^ X + d. 



(6.82) 



This is precisely the source term coupled to the Coulomb 
gas density discussed above, which appeared in Eq. 

(lonii . 

We now see that the Y-junction conductivity is directly 
related to the Coulomb gas correlation function. To make 
the connection we use: 



to invert l|6.8()(l : 



a = ^ Kjtju{Ak - A/)/3. 



r5(t^„+c^:jZ,fe(u;„), (6.78) This implies: 



da^' 2 



dAk 3 



3^ -.-J 



e,fc. 



(6.83) 



(6.84) 



(6.85) 



Here we use an anti-symmetric 3x3 tensor, ajfc, defined 
by: 



E 



(6.86) 



with ei2 = £23 = £31 = 1. We may thus write: 



The derivative on the right hand side of Eq. H6.87|l is proportional to the Coulomb gas correlation function. 
{PtJ.{^n)pv{t^'n)) as we observed in Eq. (|6.37|) . 

There are two simple limits of the DHM. One is the case where the potential, V, is irrelevant. Then the correlators 
of the X-fields are determined by So- 

{X^~Lu^)X\-u;'„)) = (27r)2<5(,.„ + c.^) (i?-i)-^(-c.„), (6.88) 
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and so the Coulomb gas correlation function vanishes. 
This corresponds to the limit of disconnected wires and 
zero conductance. The other simple limit the corresponds 
to the potential, V, renormalizing to infinity and pinning 
the X co-ordinates at one of its minima. The Coulomb 
gas correlation function is then given by the first term 
only in Eq. Thus, 



Zjk{uJn) 



ej^efenif^if;: D^,i-Un). (6.89) 



Using this gives: 

2 2 
Zjk{u!n) = -^{3Sjk - l)a\u!n\ + —^Pujnejk- (6.90) 

Thus, from Eq. (|6.79|) . the dc conductivity is given by: 



Gjk - ^ 



^(3<5j7c 



(6.91) 



We have restored the factors of and h which we had 
previously set equal to one. The values of a, (3 for the 
chiral fixed points are obtained from Eq. with n = 

1 and n — 0, and lead to the conductance tensor 



2g 



h (3 + 5^) 



[{35jk - 1) ±gejk], 



(6.92) 



with ± corresponding to the fixed points x±- It is in- 
teresting to check that the same result can be obtained 
by applying electric fields uniformly to the entire wire j, 
which we do in the next section, where we address the 
Y-j unction problem using BCFT. 



VII. APPROACHES BASED ON BOUNDARY 
CONDITIONS 

After studying the junctions of three quantum wires 
using the mapping into the dissipative Hofstadter model 
in the previous section, we now turn to a study based on 
boundary conditions satisfied by the bosonic fields de- 
scribing the degrees of freedom of the wires. The bound- 
ary condition of the boson fields represents the physics at 
the junction. In particular, conformally invariant bound- 
ary conditions corresponds to RG fixed points of the junc- 
tion problem. Among them, stable fixed points are of 
special importance, as they govern the physical proper- 
ties of the junction at low bias voltages and low tem- 
peratures. From these boundary conditions, one can de- 
termine the conductance tensor of the junctions as well 
as the scaling of perturbations to the fixed point, i.e, 
the perturbing operators and the physical processes that 
they correspond to. The study of conformally invariant 
boundary conditions will allow us to identify other fixed 
points in addition to the ones discussed in section IVTl 

We carry out the analysis based on boundary condi- 
tions in two ways. The first approach is based on the 



formalism of boundary conformal field theory (BCFT), 
which is a systematic treatment of conformally invariant 
boundary conditions for a given (bulk) CFT. The sec- 
ond way is an approach that we introduce in this paper, 
which we refer to as the method of delayed evaluation of 
boundary conditions (DEBC). 

A theory with free c-component bosons is a conformal 
field theory with central charge c. The classification of 
all possible conformally invariant boundary conditions, 
even for free bosons, is a rather difficult and rich problem 
especially for c > 2. For c > 2, it is possible to have 
a "magnetic flux" at the boundary, as was discussed in 
Refs. /3] in the context of Dissipative Hofstadter Model. 
Nontriviai boundary conditions, which have no simple 
description in terms of the original boson field, have been 
discussed in Refs. 0, Is^l for c = 2 free boson. As the 
junction of 3 wires is mapped to the boundary problem 
of c = 2 free bosons, we can expect these features to 
be relevant in the present problem. App. [CI and IdI 
review aspects of BCFT, BCFT for free bosons and the 
calculation of the conductance using BCFT, respectively. 

In the bosonization approach, the Fermi statistics of 
the electrons may be taken care of by the Klein factors. 
On the other hand, in the BCFT, conformally invariant 
boundary conditions, in particular the nontriviai ones, 
are formulated in terms of boundary states. Boundary 
states belong to the Hilbert space of the theory with 
the periodic boundary condition in the "space" direction. 
The implication of the Fermi statistics of electrons on the 
boundary state formulation is not obvious. 

Considering the Fermi statistics of the electrons, we 
show in section IVIIII that the structure of the Hilbert 
space is "twisted" compared to that of the standard free 
boson. The twisted structure affects the possible bound- 
ary states, the scaling dimensions of the boundary oper- 
ators, and stability. We show that the derived boundary 
states are indeed consistent with electrons with Fermi 
statistics. Thus the twisted structure encodes the effect 
of the Fermi statistics. 

The formulation of the problem of tunneling between 
multiple quantum wires within the DEBC approach, pre- 
sented in section IIXI provides another angle of attack, 
which gives results consistent with the BCFT approach. 
The DEBC method allows us to give a simple physical 
interpretation of the operators that become relevant and 
enforce the correct boundary conditions, as well as the 
operators that correspond to perturbations around the 
fixed point boundary condition. 

The DEBC method consists of doubling the number of 
degrees of freedom in the problem by keeping both the 
bosonic fields ipi, i = 1, 2, 3 in each lead as well as its con- 
jugate momentum 6i. Notice that such construction, in 
terms of both sets of bosonic fields (pi and 9i , is the usual 
one when describing systems on the infinite line, but is 
unusual if we want to describe the system on the half line 
a; > 0, appropriate for describing the quantum wires that 
meet at a junction. There must be a boundary condition 
imposed at a; = which cuts in half the number of degrees 
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of freedom, when describing the system on the half-hne. 
It is only with this imposed boundary condition that the 
operators above describe the electronic degrees of free- 
dom on the the half-line. The essence of the method is 
that we delay the choice of boundary condition (hence 
the name "delayed evaluation" ) until after we write the 
operators corresponding to the physical tunneling of elec- 
trons (and pairs or multiplets), working with twice the 
physical degrees of freedom. With the tunneling oper- 
ators in hand, we can then study the correct boundary 



condition that must be imposed on ipi and 9i. 

In Sec. 13 we use both methods to discuss various fixed 
points of the Y-junction. 



VIII. FERMI STATISTICS AND THE TWISTED 
STRUCTURE IN BOSONIZATION 



The BCFT of free bosons can be applied to the system of interacting electrons via bosonization. However, owing 
to the Fermi statistics of the electron, there are various differences compared to the standard compactified boson field 
theory, reviewed in App.jCl In this section we formulate the effect of the Fermi statistics in terms of a twisted structure 
in the Hilbert space of the free boson theory. This approach was first developed in Ref. [s^. In this paper, we further 
propose a simplified treatment of the twisted structure for the junction problem, and apply it to the junction of three 
wires. 



A. Bosonization of a single wire 



Let us start from the bosonization of a single wire, carefully taking the effect of the Fermi statistics into account. 



1. Derivation of the twisted structure 



Since the Matsubara formalism implies anti-periodic boundary conditions in imaginary time, anti-periodic boundary 
conditions in the space direction will also occur after a modular transformation. These must be imposed independently 
on left and right- movers: 

i^L/R{t.x)=-il;L/R{t,x + (3). (8.1) 

It will be important to determine the corresponding b.c.'s on the bosons, using the bosonization formulas reviewed in 
Sec. IV. Let us first consider the case, g = 1, corresponding to free fermions. The b.c.'s are then: 

eM^^VL/R{t.P)] = -exp[z\/2(^i/j^(t,0)]. (8.2) 

Because ['-Pl/RtVl/r] 7^ it is not so straightforward to deduce the boundary conditions on the boson fields from 
Eq. (|8.2|l . One way of proceeding is to determine the b.c.'s on the bosons which reproduce the finite-size spectrum 
implied by anti-periodic b.c.'s on the fermions. The allowed wave- vectors for single-particle fermion states are: 

27r(n-H/2) 
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(neZ). (8.3) 



Recognizing that k is measured from ±kp for right and left-movers respectively, we see that n may be positive or 
negative in Eq. (|8.3() for both left and right movers. The negative n states correspond to hole states for right-movers 
and electron states for left-movers. It is straightforward to calculate the energy of an arbitrary multi-particle fermionic 
state with this linearized dispersion relation. Consider the energy of an arbitrary state of the right-movers. The lowest 
energy state of charge Qr S Z, is: 

EiQR)='i E (-+1/2) = ^%. (8.4) 

Here Qr may be positive or negative, corresponding to creating particles or holes. An arbitrary multi particle-hole 
state can be constructed, on top of the charge Qr groundstate, by first raising the lowest rim fermions by m levels, 
then raising the next rim-i by m — 1 levels, etc., finally raising ni fermions by 1 level. Thus the energy of an arbitrary 
excited state of the right-movers may be written: 
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where Qn is an arbitrary integer and the arbitrary non-negative integers. 

Precisely this spectrum arises from right-moving bosons. We assume that the right and left moving boson field, 
4>R/L are periodic variables: 

'PR/L{t,x + (3) (PR/L{t,x) + V2TTn, {neZ). (8.6) 

This is reasonable since all physical operators, local in fermion fields, are single-valued with this identification. We 
then find that the anti-periodic b.c.'s on fermions arise from periodic b.c.'s on bosons. To see this, consider the mode 
expansion for periodic right-movers: 



TT 1 °° 1 

^ n—1 ^ 



R r - ^'^1 1 

exp[— J H- n.c. 



(8.7) 



Here Pr, the momentum variable conjugate to (pQ, has eigenvalues times an integer, consistent with the assumed 
periodic b.c.'s. Substituting this mode expansion into the right-moving Hamiltonian: 



1 

Hr = ^J^ {d^^Rf dx, 



gives precisely the spectrum of Eq. H8.5|l . upon the identification, Pr <-> ^/2Qr. Thus we see that anti-periodic b.c.'s 
on fermions turn into periodic b.c.'s on the chiral boson field. 

Alternatively, it is also possible to derive the periodic b.c.'s on right-moving bosons from anti-periodic b.c.'s on right- 
moving fermions from Eq. (|8.2|l . taking into account the boson commutation relations. The equal-time commutator 
follows from Eq. H8.7|l : 

[ipRix),^Rm = iin/2)ep{x), (8.9) 

where ep{x) is the function, anti-symmetric in x, which is constant everywhere except at a; = nf3, for n ^ Z, where it 
jumps by 2: 

ep(x) = (2n+ 1), [nP < x < {n + l)/3]. (8.10) 

We multiply both sides of (the "right version") of Eq. H8.2|l on the right by exp[—i^/2ipR{S)] where 5 is a small 
positive quantity, giving: 

eM^V2ipR{p)] exp[-z\/2^fl(<5)] = - exp[i\/2(^j^(0)] exp[~iV2ipR{S)]. (8.11) 

This can be rewritten: 

eM^V2[^R{l3) - ipR{S)]}exp{[ipR{P),^R{S)]} - - exp{^V2[(pfl,(0) - ^Mi^R.iO), miS)]} (8.12) 

Since < /3 — (5 < /? and — /3<0-(5<0, using Eq. H8.9|l this becomes: 

i exp{iV2[(pR{f3) - 'fR{6)]} = i exp{iV2[(pR{0) - ^r{S)]} (8.13) 
We may finally let 5 — + 0+, giving: 

exp{iV2[^RiP) - ipRiO)]} = 1, (8.14) 

and hence periodic b.c.'s: 

ipj^{f3) = ^pJ^{0) (mod V2Trn). (8.15) 

Note that if we had instead chosen the infinitesimal S < 0, so that (3 < (3 — 5 < 2(3 and < — (5 < /3, the two 
commutators in Eq. (|8.12|) would still differ by tt from Eq. H8.10|l so that we would still obtain Eq. H8.14|l . 
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2. Finite-size spectrum and bulk operators 



The left-moving fermions with anti-periodic b.c.'s have an identical spectrum to the right-movers which again is 
given by periodic b.c.'s on the left-moving bosons. Adding left and right contributions, we obtain the mode expansions: 



(fa 



[{QL+QR)t+iQL~QR)x] 



1 °° 1 

— Y — 



R I ^TT T 27r 

a„ exp \—inx- — ) + a„ exp {—inxj^ — ) 



0{t,x) 



V2t: 



[iQL-QR)t + iQL+QR)x] 



a„ exp [-mx^ — ) - a„ exp [-inx+ — ) 



(8.16) 



Here 



While it is convenient to define a new set of integer quantum numbers: 

Q = Ql+Qr 

Q = Ql — Qr, 



(8.17) 



(8.18) 



it is important to realize that they do not run independently over the integers, but rather obey the "gluing condition" 



Q^Q (mod 2) 



(8.19) 



Thus, the periodic identifications of the fields ip(t,x) and 9{t,x) do not simply correspond to them being angular 
variables. They are of the form: 



If ^ ip ->r Qypi'K, 



.20) 



but rather with the gluing condition (|8.19() . We refer to this as a "twisted structure" . It simply refiects the fact that 
the bosons arise from bosonizing fermions. It can be verified that the finite size energies (with anti-periodic b.c.'s 
on the fermions) are in one-to-one correspondance with the (bulk) operator content. For example, the right-moving 
oscillator vacuum with quantum number < corresponds to the lowest dimension (right-moving) operator with 
charge Q^. Taking into account Fermi statistics, this operator is: 



Qr-1 

n 

3=0 



fR — J- 

n [(^^ 



with the (right-moving) conformal dimension 



Qr-1 -i 



Qr^ 



.21) 



.22) 



Comparing with eq. 1)8.41) . we see that the usual relationship between the finite size energy and scaling dimension. A, 
is obeyed: E — (27r//3)A. The corresponding bosonized operator is exp[z(5fl-\/2</'_R]- 

We have discussed extensively the case g — I, corresponding to free fermions. However, these results can be seen 
to carry over to the case of interacting fermions, with general values of g. In general, the mode expansion for if is: 



ip{t, x) 



^0 



V27r 
IT 
1 1 



-Qt + Qx 
9 



R I ^TT r 27r 

a„exp(-^n..-)-,a„exp(-.n.,-)+h.c. 



(8.23) 
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and the corresponding spectrum is: 



again with the gluing condition l|8.19|l . The corresponding bulk primary operators are: 

exp 



with dimensions: 



4(7 4 ■ 



(8.24) 



(8.25) 



(8.26) 



When we form boundary states from the Hilbert space with the anti-periodic b.c.'s on the fermions, the condition 
\Hr — T-LlWA) = wiU simpUfy this twisted structure, as we will see. 



3. Boundary states and boundary operators 



Let us now consider the boundary operators, finite size spectrum and boundary states. We consider BC's on the 
fermions either of the free type: 



V'L(i,0) = 7Afl(i,0), 

corresponding to 9{Q) — 9q (a constant), i.e. N b.c.'s on Lp or else of Andreev type: 
corresponding to D BC's on ip. The N EC turns the operators of Eq. (|8.25|) into: 



of dimension: 



exp [-!-{Q^{t, 0) + Qe{t, 0))] ^ exp [-!-{Qip{t, 0) + Q6q)] cx exp[iQ\/2v3i(t, 0)], 



(8.27) 
(8.28) 

(8.29) 
(8.30) 



We find that the value of Q is now immaterial to the scaling dimension. Consequently, we can allow Q to take any 
integer value, even or odd, despite the gluing condition (|8.19|) . Also note that, in the free fermion case, g = 1, we 
obtain A^r = (5^/2, corresponding to the operator spectrum of a single left-mover, which we can obtain by regarding 
the right movers as the continuation of the left-movers to the negative x-axis. 

As reviewed in IVIII the spectrum of the boundary operators is related to the FSS on the finite- width strip with 
the same boundary condition imposed on the both ends. To find the corresponding FSS, let us consider the strip 
< a; < Z, and impose the N b.c. on (f at the both ends a; = 0, L We notice that the N b.c. (on ip) corresponds to the 
D b.c. on the dual field 9. The mode expansion consistent with the D b.c.'s on 9 



9{{)) = (mod 

9{l) = (mod V27r) 



reads: 



This implies the Hamiltonian: 



^ OO 

y 2 sir 



H 



2<7 



—j— [a„ exp {—int—) + h.c. 



.31) 



.32) 



(8.33) 



29 



This is indeed consistent with the scahng dimension of the boundary primary operator in Eq. (|8.30() . 

On the other hand, as discussed in Sec. IV ill the boundary condition generaUy can be described in terms of the 
corresponding boundary state upon modular transformation. Naturally, the N boundary state for the present case is 
similar to the N boundary state of a standard free boson. However, here the gluing condition H8.19|l is important in 
determining the boundary state. Since the N boundary state (on ip) is Dirichlet boundary state on 9, the winding of 6 
along the boundary is zero. Consequently, we have the quantum number Q = Q. Because of the gluing condition H8.19|l . 
the quantum number Q is now restricted to even integers. The N boundary state is: 



Here the ' on the summation indicates that Q is restricted to even integers. Upon performing the modular transfor- 
mation we obtain: 



Here m is summed over all integers. This corresponds to the finite size spectrum of Eq. ()8.33|l . demonstrating the 
consistency. 

We can repeat the similar analysis for the D b.c. on ip. The primary boundary operators have dimensions: 




(8.36) 



where Q can be even or odd, and the mode expansion becomes: 



oo 




-I Q 



(8.34) 




(8.35) 




(8.37) 



where Q can be even or odd. This implies the Hamiltonian: 




(8.38) 



The D boundary state is: 



oo 




1/4 



exp 




(8.39) 



n=l 



J Q 



where the sum is restricted to even Q since Q is now set to zero. 



We now obtain: 




(8.40) 



consistent with the FSS of Eq. (|0^ . 



B. Junction of two wires 



Now let us consider the implications of the twisted structure for the junction of two wires. 



1. Structure of the boundary states 



We introduce 2 fermion fields tpj^L/R (j — 1; 2), corresponding to each wire. In order to apply the boundary state 
formalism, we define them on a circle of circumference p with the antiperiodic boundary condition. Bosonizing each 
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fermion fields, we introduce 2 boson fields ipj and their duals, 9j. They independently obey the twisted periodic 
identification of Eq. 



It is convenient to define linear combinations: 



QjV27r 
(mod 2) 



..41) 



$0 = 

eo = 



V2 

h+02 

V2 



and 
and 



e = 



9i — 62 



(8.42) 
(8.43) 



As already discussed in Sec. IVIII as long as the charge is conserved at the junction, <i>o always obeys the Neumann 
boundary condition. The Hamiltonian will thus decouple into separate terms for $ and <i>Oj with $0 decoupling from 
the interactions. Namely, only $ and O will appear in the boundary interactions. 

However, the periodic identification conditions (referred to as "gluing conditions" in [s^) now link these 2 fields, 
as well as linking the <!> and 8 fields. Explicitly, defining: 



Qo = Q1+Q2 
Q = Qi — Q2, 



(8.44) 



corresponding to <I>o and $, and similarly defining Qq and Q, we see that these obey the restrictions, which are 
equivalent to Qj = Qj (mod 2): 



Qo = Qa = Q = Q (mod 2) 
Qo + Q + Qo + Q = (mod 4). 

the mode expansions for the fields $0 and $ are: 



..45) 



$o(t, x) 



-Qot + Qox 
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-Qt + Qx 
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RO 



1 ^ 1 



exp \—%nx^ — ) 



27r 



,L0 



exp {—mxj^ — ) + h.c 



27r, 



a„ exp {—inx-——) + a„ exp (— ma;_|_ — ) + h.c 



(8.46) 



and similarly for Oq, 6. The Q, Qqj Q and Qo quantum numbers obey the gluing conditions of Eq. H8.45f) . 
The ESS is given by two copies of Eq. (|8.24|) . This can be written as: 



E = 



27r 
1 



LO 



(8.47) 



where, again the gluing conditions of Eq. H8.45() are obeyed. 

Analogous to what was found in Sec. IVIII in the 3- wire case, $0 always obeys N BC's at a; = 0. This is the "bare" 
BC corresponding to decoupled leads, and it is not changed by the interactions since $0 does not appear in them. 
On the other hand, <i> can obey N or D BC's corresponding to perfect refiection or perfect transmission through the 
junction, respectively. 

It is straightforward to determine directly the full boundary operator spectrum with either set of b.c.'s. (Henceforth 
we will label these two sets of b.c.'s as simply N or D.) In the N case we simply get 2 copies of the operator spectrum 
of a single wire with N b.c.'s, since (pj {j = 1, 2) independently obey N b.c.'s in this case. As usual, this spectrum of 
operator dimensions is in one-to-one correspondance with the ESS with N b.c.'s at both ends of a strip of length /. 
This is, as in Eq. 1H!S3: 



E 



N 



Qi , Q2 



(8.48) 



31 



where Qi and Q2 range independently over all integers. Equivalently, we can write this in terms of the charge quantum 
numbers Qq and Q defined in Eq. (|8.44l) and the new oscillator operators, which are associated with the fields $0 
and $: 



En ~ T 



where now: 



Qo,Q^, ST^ r 



Q^Qo (mod 2) 



(8.49) 



.50) 



The N boundary state is just a product of the single fcrmion N boundary states. Taking the case where both phases 
Qq and 8q are set to zero, this can be written as: 



\N) 



exp 



^(aft„^ot+,Lt„«t) 

.n=l 



E l(Qo,o))®|(g,o)). 



.51) 



Here, as follows from Eq. H8.45|l . with Q = Qo = 0, Qo and Q are even and Q + Qo = (mod 4). 

The FSS with the D b.c.'s at both ends of the strip is equivalent to the FSS of a single wire with periodic EC's on 
an interval of length 21, as follows from an unfolding transformation [3^. This is given by Eq. H8.24|l with f3 replaced 



by 21. We may also replace the n. 



L 
m ' 



quantum numbers by a new basis, n^, and rim, and identify the total charge 



Q with Qo and Q with the same quantity related to the linear combination (p in Eq. 



Ed = t 



with 



17 + — + L "^("™ 

^ 7n=l 



Qo = Q (mod 2). 



+ rim) 



The corresponding boundary state is: 



\D) = exp 



E l(Qo,0))®|(0,Q)), 
Qo,Q 



.52) 



.53) 



(8.54) 



where the sum is restricted by the gluing conditions of Eq. H8.45|l . in the case Qo = Q = 0, to Qo and Q even with 
Qo + Q = [mod 4]. 

In fact, we are generally only interested in the Qo = boundary operators since these arc the only ones which 
can appear in the Hamiltonian. Operators with Qo 7^ corresponds to a vertex operator involving Qq and thus are 
forbidden by charge conservation. Setting Qo = in Eq. (|8.49|) or (|8.52|l we see that Q or Q respectively must be 
even. The dimensions of the primary operators can thus be written as: 



Am = — 



gm 



(8.55) 

with n and m arbitrary integers. These are precisely the same as the boundary scaling dimensions of primary fields 
for a single periodic boson with N or D EC's [ Eq. (|C44ll or (|C23|I .] These reproduce the known results on the junction 
of two wires. QHHH 

We note that, however, in the presence of a resonant level at the junction, the charge-nonconserving boundary 
operators with Qo 7^ should be also considered. Suppose for example that there is a single resonant level, exactly 
at the Fermi energy, at the junction. One may imagine the two wires connected to a quantum dot. Let d and d^ be 
the annihilation/creation operator of an electron in that level. Then the hopping of an electron between a wire and 
the dot can be written as ipjd^ + h.c. As d and d^ have zero scaling dimension, the scaling dimension of the hopping 
operator is determined by the electron annihilation operator ipj at the end of the wire. Therefore, in the presence of 
the resonant level, even if the charge is overall conserved, we have to consider the "charge-nonconserving" boundary 
operators. Indeed, the effect of the resonant tunneling at the junction has been studied in Refs. QjE^ for the case of 
two wires and in Ref. for three wires. 

Nevertheless, in this paper, we focus on the problem without the resonant levels. We discuss only the charge- 
conserving boundary operators, for which we develop a simplified treatment in the following. It should be kept in 
mind, however, that it does not cover the resonant case. 
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2. Projection to an effective single component boson 



While we have succeeded in correctly reproducing known results on the junction of two wires, the gluing conditions 
obviously complicates the calculation. On the other hand, if we restrict the analysis to the charge-conserving boundary 
operators, there is a simpler method to find the dimensions of such boundary operators. This will be generalized in 
the next sub-section to the 3-wire case. 

Note that in any boundary state which is N for $0, only Qo = states occur. The restrictions of Eq. H8.45|l reduce 
to: 



Qo = Q = Q = (mod 2) 
Q + Q + Qo = (mod 4). 



(8.56) 
(8.57) 



Defining: 



Qo 
Q 
Q 



2no 

2n 

2n, 



these integers obey the condition: 



n + h ^ fiQ (mod 2). 
A general boundary state of this type can be written: 

4m+2,4m+2/- 



(8.58) 

(8.59) 
(8.60) 



Here j-^g g) are states in the Hilbert space of <f>, proportional to \{Q, Q)). Actually only |^2fi,o) is non-zero in the N 
state and only |V'o,2n) is non-zero in the D state but we since we wish to consider both cases at once we use this more 
general notation. IA^q) and \Nq) are N states in the Hilbert space of $0, given explicitly by: 



I AO) = (2g)i/4exp 
\Nl) = (2g)i/4exp 



00 



E l(4'?io,0)) 



E |(4mo + 2,0)). 



(8.61) 



Here we have let fio ~ 2too in |A^o) and no = 2mo + 1 in |A^^). The diagonal partition function is given by: 

ZbB — Zj^o E] ^4m,4m + Zj^l E] ^4m-|-2,4m-|-2 • (8.62) 



m.rn£Z 



Here Z^^i is a diagonal partition function in the $o Hilbert space and the second factors are diagonal partition 
functions in the $ Hilbert space: 



^2n,2n = (V'2n,2n | exp[-ZiJ^] |V'2fi,2n) 



.63) 



The former are given explicitly by: 



ZatO 



Qoez 
1 



Qo„QoV(49) 



.64) 



Qo&Z 



The complete spectrum of boundary operators can be read off from the partition function Zbb- However, the 
scaling dimension coming from Zjyo and Zj^i corresponds to a boundary operator involving the center-of-mass field 
$0- In particular, those includes a nontrivial dimensions Qq'^ /(Ag) comes from the vertex operator e*'5o*o/2 -^^j^j^j-^ 
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changes the total charge. Thus they do not appear if we consider only the perturbations which respects the charge 
conservation at the junction. The other operators with integer dimensions corresponding to the expansion of l/r]{q) 
comes from derivatives of $o- They are all irrelevant or marginal. 

Thus, in order to identify potentially relevant operators, we can consider the boundary operators which involves 
the dynamical field $ only. For this purpose, we may replace both Z^o and Zj^i by 

Zno,Z,mi 1 (8.65) 

The resultant effective partition function is 

ZbB ^^2ri,2n, (8.66) 

n,n 

from which we can read off the spectrum of boundary operators involving $ only. 

This "reduced partition function" can be obtained by simply ignoring the $0 sector and dropping the constraint of 
Eq. H8.57|l . thus allowing n (= Q/2) and fi {= Q/2) to be arbitrary integers. Namely, Eq. (|8.66|) is identical to the 
diagonal partition function arising from a "reduced" boundary state entirely in the $ sector of the form: 

|B)flCX^|l^2S,2n), (8.67) 
n,n 

in which the second constraint of Eq. H8.57|l has been simply ignored. 

Thus, to obtain the dimensions of boundary operators in $ sector, we may simply forget about the second constraint 
of Eq. (|8.57|) and hence replace the entangled boundary state of Eq. H8.60|) by the simpler reduced one, written only 
in the $ sector, of Eq. H8.67|l . This corresponds to a boundary state for a single boson obeying standard periodic 
b.c.'s. In particular it leads to the primary operators spectra of Eq. (|8.55|) . Therefore, as long as we restrict ourselves 
to Qo = operators, we may drop the 4>o boson entirely and ignore the twisted structure that results from Fermi 
statistics. This simplicity appears to be related to the fact that the Klein factors are trivial for a 2 wire junction. 

In contrast to the above, we will see that a similar reduction leaves a nontrivial effect for the junction of three wires. 



C. Three wires 



Let us finally turn to the model of interest here: a 3- wire junction. As in Subsection IVIII BI on the 2- wire junction, 
we define 3 fermion fields on a circle of circumference (3. Upon bosonizing, we introduce 3 fields and their duals, ipj 
and 9j, which obey the fermionic type (twisted) periodic identification of Eq. (|8.41l) where now j = 1, 2 or 3. We then 
introduce the 3 linear combinations of fields $j as in Eq. H4.5|l . We then combine <I>i and $2 into a 2-component vector 
field, $ and similarly for Q. We now consider the mode expansions for the fields <&, $0, © and Qq. We introduce 
charge and dual charge quantum numbers Qj and Qj for each of the original boson, (pj, as in Eq. (|8.23|) . H8.18|) . We 
define the total charge and dual charge: 

Qo = Q1 + Q2 + Q3, 

Qo = Qi + Q2 + Q3- (8.68) 

We are only interested in b.c.'s which are N for ^q, so that the corresponding b. states are constructed from Qo = 
states. It is convenient to parameterize the vectors of integers, Qi satisfying this condition by: 

(Qi, Q2, Q3) = "^l(0, 1, -1) + m2(-l, 0, 1). (8.69) 

Clearly mi — Q2, m2 = —Qi and letting the run over all integers uniquely parameterizes all Qi satisfying Qo = 0- 
On the other hand, Qq need not be zero and it is convenient to parameterize the Qi in a different way: 

(Qi, Q2, Q3) = «o(l, 1, 1) - m(0, 1, 1) - n2(l, 0, 1). (8.70) 

This can be inverted to give: 

(no, m, "2) = (Qi + Q2-Q3, -Q3 + Qi, -Q3 + Q2). (8.71) 

Thus any set of integer dual charges, Qj, can be uniquely expressed in terms of the integers, nj. The conditions 
Qj ~ Qj (mod 2) imply: 

no = 2po (i-e. no is even) 

Uj = n-ij (mod 2). (8.72) 
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In this parameterization, the total dual charge is: 

Qq = 2(3po - ni - 712). 

This implies that Qq is even and: 

Qa/2 = —ni — n2 (mod 3). 
Thus we see that a general boundary state which is N with respect to $0 will be of the form: 



(8.73) 
(8.74) 



\B)= J2 K 

fc=-l,0,l 



^ |i/)(ni,n2,TOi,m2)). 

ni+n2 — — :? {mod 3) 



(8.75) 



Here the ' over the summation sign indicates that, in addition to the explicitly noted restriction, the conditions: 

Hk = TUk (mod 2), j = 1,2 (8.76) 

must also be imposed. The states 1-0(71,1, 712, mi, 7712)) are in the sector of the 2-component field $, with charges and 
dual charges labelled by the integers 77ij and rij. The states I^Vq), for k = —1,0 or 1, are in the $0 sector and are 
given explicitly by: 



IK') - ( y 



1/4 



exp 



i?,Ot 



El(2(3fo + fc),0)). 



.77) 



The corresponding "partial Neumann" partition functions in the $0 sector may be readily calculated. 



Qo&z 



As in Subsection IVIII Bl we are only interested in the boundary operators in the $ sector, and replace Zj^k by unity. 

We may equivalently get the dimensions of all such primary boundary operators by using a reduced boundary state, 
which exists in the reduced Hilbert space of the 2 component boson field, $. That is: 



\B) E |V'(^^l,?^2,"^l,7^^2)). 



(8.79) 



The restrictions of Eq. H8.76|) remain, but the other conditions, that entangle the $ sector with the $0 sector disappear 
upon restricting to boundary states with Qq = and boundary operators with Qq — 0. 

The windings A $ = $(/3)— $(0) and A 8 = 8(/3)— 8(0) along the boundary follow from the integers (771,772,7711,7^72). 
A$ and A8 can be obtained from the windings of the three pairs of boson fields, Ayj^ — Qi\/2TT and A9i = QiV27r, 
using the rotations in Eas. l|4.5l4.6|) and the vectors Kj and Rj defined, respectively, in Eas. H4.9|l and H().55|l . One 
obtains 



$(/?)- $(0) = 27r(77i:^+ 772 ^) 



(8.80) 

e(/3)-e(0) = 2Tr{miKi+m2K2). (8.81) 
Thus the "reduced boundary state" H8.79|l would arise from a 2-component boson field with the mode expansion: 
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1 ^ 1 ^ 



1 °° 1 

— Y — 



27r, 



27r, 



a;^ exp[— 77ia;+— J + a^^ exp[— ma;_ — J + h.c. 



J = l,2 i=l,2 

where the integers rij and 777^- are restricted by Eq. (|8.7t)|) . The corresponding energies are 



E : 



27r 

T 




(8.82) 



(8.83) 
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where the . . . represents the oscillator mode terms. Of course the restriction of Eq. (|8.76() must be applied. 

Let us compare the mode expansion (|8.82|l to the corresponding Eq. (|C53|I for the standard compactified boson, with 
A being the triangular lattice spanned by Ri and R2. While they are similar, we find that they are in fact different 
and cannot be identical with any choice of the compactification lattice A for the standard free boson. Namely, the 
fermionic nature of the problem has modified the mode expansion, via the factor of 1/2 multiplying the second term 
inside the first square bracket in Eq. H8.82() and the restriction of Eq. (|8.76|) on the quantum numbers rii and m^. 
In contrast to the two wire case, the fermionic nature of the problem survives and makes important changes to the 
boundary states and the corresponding "reduced mode expansion" of Eq. (|8.82|l . This may be related to the fact that 
the Klein factors cannot be eliminated for the 3-wire juncti on whereas they can be eliminated in the 2-wire case. 



IX. THE METHOD OF DELAYED 
EVALUATION OF BOUNDARY CONDITION 
(DEBC) 

In this section we introduce the DEBC method. The 
method consists of doubling the number of degrees of 
freedom in the problem by keeping both the bosonic fields 

and its conjugate momentum 9i, for each semi- infinite 
lead i. The redundancy is eliminated by imposing the 
boundary condition at a; = 0, which cuts in half the 
number of degrees of freedom, down to the correct num- 
ber for the semi-infinite wires. The basic idea behind the 
method is that we delay the choice of boundary condition 
until after we write the operators corresponding to the 
(multi)particle tunneling processes at the junction. The 
correct (or stable) boundary condition that must be im- 
posed on and 6i is determined a posteriori, taking into 
account the scaling dimensions of the tunneling operators 
given the different choices of boundary conditions. 

Explicitly, we write the electron operators, both right 
and left movers (or incoming to and outgoing from the 
junction at x = 0), for each of the semi- infinite leads, 
labelled by i, defined in the half line x > 0: 

Wj cx e v5^^' " . (9.1) 

Klein factors can be absorbed into the zero mode terms of 
the mode expansion of the ipi and 9i bosons. For example, 
the fermionic commutations between different leads can 
be taken care of by imposing that ['pf'\'pf'^] = i27ray , 
with an anti-symmetric a^j = —ctji and unit entries 
aij = ±1. (A different represention of Klein factors using 
zero modes is introduced in Appendix IeI when we discuss 
fractional quantum Hall junctions.) The 0^ fields can be 
taken to commute for i ^ j. 

Notice that the above construction, in terms of both 
sets of bosonic fields ipi and 9i, is the usual one when 
describing systems on the infinite line, and thus we are 
working with twice the number of degrees of freedom for 
the semi-infinite line. What is missing is the boundary 
condition at the origin x — (the choice of which we 
are delaying) which contraints ipi and 0i and thus cuts in 
half the number of degrees of freedom for systems on the 
half-line. 

The subsequent step is to construct all tunneling op- 
erators between the leads using the electron operators in 
Eq. (|9.1|) . keeping both ip and 6 fields, and imposing later 
the appropriate boundary condition. 



For concreteness, let us present below the DEBC 
method through an example, using the simpler problem 
of the two- wire junction for this purpose. In the next 
section, we apply the method to the problem of interest 
in this paper, the three-wire junction. 

For the two- wire problem, the single particle tunneling 
terms at a; = can be expressed in terms of the ^'^ ' as: 

Ti"^ = vl/fl^^^ cx e'^- e'— (9.2) 

r r r t r I i vi-V2 ; ^1-^2 

Tii^ =v|/^T v|/[|^^^(xe* (9.3) 

The fields and 9i above are all evaluated at the bound- 
ary X — Q. Other terms can be obtained by simply taking 
the Hermitian conjugate or by exchanging RL in these 
three combinations above. Physically, these terms cor- 
respond to tunneling between the leads [terms (|9.2() and 
(|9.3|) ]. and backscattering [term H9.4|l ]. 

Written in terms of the rotated basis Eqs. (|8.42I8.43|1 . 
we have 

T^i^ oc e'* e*®° oc e** (9.5) 
T^^^ oc e^* e'® (9.6) 
T/J^ oc e'(e+eo) ^ ^ (9 j) 

where we have used the fact that charge conservation re- 
quires that $0 must obey the Neumann boundary condi- 
tion, or equivalently, Oq must obey the Dirichlet bound- 
ary condition (Og is pinned at the boundary). 

The boundary conditions on the other pair of bosons, 
$ and 8, can be expressed in terms of chiral fields as 

0fl - n^L + C, (9.8) 

with $ = -^(0;j -I- 0i) and 6 = ^ {(j)R - <t>L), where Tl 
is a rotation and C is a constant. Since we just have a 
single pair of RL fields, the only possibilities are 7?. = ±1. 
We will soon see that in the problem of three leads there 
will be two pairs of RL boson fields, and more general 
rotations TZs are allowed. 

The scaling dimension of a general vertex operator 

is given by 

Ay = i|7^(a + 5) + (a-6)|2 . (9.10) 
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For N boundary conditions, Q is pinned, TZ = 1 and 
Ay — a?. For D boundary conditions, $ is pinned in- 
stead, TZ — —I and Ay = 5^. These results for Ay could 
have been easily guessed, but going through the exercise 
with the rotation TZ is useful as a warm up for the three 
lead case. 

We now analyze the scaling dimensions of the three sets 
of operators H9.5I) - (|9.7|) for different boundary conditions 
on $. 

• Neumann boundary condition - In this case 

At«^ = 5"' , ^T,\^ = ff"' , ^T^^ - 0. (9.11) 

The operator is proportional to the identity (dimen- 
sion zero), which is consistent with the fact that it is this 
backscattering operator that is responsible for pinning 
the O field for N boundary condition (on 4>) . This is the 
perfect backscattering fixed point, stable as long as the 
other two operators and are irrelevant, which 
is the case for g < 1. 

• Dirichlet boundary condition - In this case 

Aj.Ri, = , ^T^^'^ — 9 ^ '^T/l^ = 9- (9-12) 

Now it is the operator T^'" which is proportional to the 
identity (dimension zero), and this transmission opera- 
tor is responsible for pinning the (p field for D boundary 
condition. This is the perfect transmission fixed point, 
stable as long as the other two operators Tg^^ and 
are irrelevant, which is the case for 5 > 1. 

It is clear from this approach that the two-wire junc- 
tion can be understood in terms of the boson fields $ and 
Q, which are dual to each other. They can be regarded as 
subject to the standard compactification <i> ~ $-t-27r and 
9 ~ 6 -I- 27r, as can be seen from Eqs. (|9.2I9.3I9.4|I . Also, 
the DEBC approach reproduces rather simply the results 
for a junction between two wires as in the Kane-Fisher 
model P|. We now apply it to the three- wire junction. 



X. FIXED POINTS FOR THE JUNCTION OF 
THREE WIRES 
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• backscattering: 
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The LL and RR combinations lead to another group: 
• LL- RR: 
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To these sets of operators we must add their Hermitian 
conjugates. 

In the following, we consider several conformally in- 
variant boundary conditions which can be expressed as 

$R = TZc$L+C, (10.16) 

where 

7^, = (''^°^^ -%"^V (10.17) 

* \ sm^ cos 4 y ' ^ ' 



Now let us discuss the RG fixed points in the junction 
of three wires. As we have discussed, each fixed point cor- 
responds to a conformally invariant boundary condition 
of the boson field theory. 

In the DEBC approach, we need the x = tunneling 
operators between the three wires. They are easily ex- 
pressible in terms of the vectors Kj defined in Eq. H4.9|l . 
with the two-component field $ defined in Eqs. H4.5I4.6|I . 

The tunneling operators involving the RL combina- 
tions can be split into three groups: 



I = — {$R + $l) (10.18) 
e = ^{$l-$r) (10.19) 

and C represents a constant vector. Qq is pinned by 
charge conservation, and can be ignored in all operators 
above. 

Consider a general operator 

y ^ g"V9 e ^ ^i(S-b)-$R+Ka+b)-$L _ (10.20) 
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For the boundary condition (|10.16() , its boundary scaling 
dimension is: 



1 



-|7^-l(a-6) + (a + 6)p 
\ |7^e(a + fe) + (a-fe)p 
i [\a\^ + \b\^ + {n^a)-a~{n^h)-h 

-{n^a) ■b + [n^b)- d] 



(10.21) 
(10.22) 



With the scahng dimension of the tunneling operators 
in hand, we can now analyze the stability of the fixed 
points that are characterized by a choice of boundary 
condition. We will carry out this analysis in the follow- 
ing section, where we present, discuss, and parallel the 
results obtained using both BCFT taking into account 
the twisted structure of the Hilbert space and the DEBC 
approach. 



A. Disconnected fixed point ~ Neumann boundary 
condition on $ 



The choice ^ = in Eq. (|10.1()|l pins 8 = C (back- 
reflection condition). This amounts to the Dirichlet 
boundary condition on 0, equivalently the Neumann 
boundary condition on $. 

Let us calculate the scaling dimensions for the ± cycle 
operators. Since TZq = 1, Av = |ap. Take as 
an example, so a 

Eq. mm -. 



h=Kz and b ~ x enter in 



A^RL = - |i?3|2 



1 

9 



(10.23) 



The result is the same for all 6 ± cycle operators, as 
well as for the 6 LL — RR operators. Upon fixing the 
boundary condition, many of these operators are redun- 
dant; the number of independent operators can be easily 
accounted for, and we find that there are only 3 indepen- 
dent operators plus their 3 Hermitian conjugates. All 
these 6 independent operators are irrelevant for g < I, 
and so the backscattering fixed point is stable for this 
range of g. 

The case just discussed above is the simplest boundary 
conditions for the problem of the junction of three wires. 
It corresponds to disconnected wires. Indeed, Neumann 
boundary condition on $ (or Dirichlet on O) is equiv- 
alent to Neumann boundary condition on all </Ji,2,3 (or 
Dirichlet on 0i,2,3)- 

Let us now connect this to the standard BCFT ap- 
proach, and construct the boundary state and the cylin- 
der amplitude explicitly. As it is reviewed in Add. lUlfor 
the standard free boson, the Neumann boundary state 
is given by a superposition of Ishibashi states. Since the 
boundary condition for O is Dirichlet, the boundary state 



must have zero winding A0 = 0, namely r/ii = TO2 = 0. 
Let we define the bosonic Ishibashi state 



|(ni,n2,0,0))) = exp 



7«t 



1(711,712,0,0)), 



(10.24) 

where |(7ii, 772, 0, 0)) is the oscillator vacuum with the 
quantum number (771,772,7771 = 0,7772 = 0). 

The most stable Neumann-type boundary state would 
be given by the linear superposition over all the possi- 
ble Ishibashi states. The important consequence of the 
twisted structure is that, because of the condition 18.76|l . 
77i and 772 cannot run over all integers when 7771 = 7772 = 
0. Instead, they can only take even integers values. Thus 
the most stable N boundary is given by 



1^) - 9N 



J2 |(277'i,277^,0,0))). 



(10.25) 



where the prefactor gi^ represents the ground-state de- 
generacy |41| of the N boundary state. (It is not 
the Tomonaga-Luttinger parameter g introduced in 
Eq. g3J.) 

The calculation of the cylinder amplitude with the 
Neumann boundary state at the both ends then follows, 
exactly as in the case of the standard boson. The diago- 
nal cylinder amplitude is given by 



(9) 



\N) 



{N\e 
9n \ - 

7?((7)2 ^ 



(10.26) 



This can be identified with the corresponding partition 
function l|C63|l of the standard compactified free boson, 
with A being the triangular lattice of lattice constant 
2/V3, spanned by i?i and i?2. Thus we find 



gN 



1/4 



(10.27) 



because Vo(A) — The modular transform of 

Eq. HC65|I should also be given by Eq. HC66|I . namely 



(9) 



1 



E 

h.hez 



(10.28) 



From Eq. H10.28|) we can immediately read off the scaling 
dimensions of the boundary operators. The leading ones 
have the smallest scaling dimension 1/g, and correspond 
to the ± cycle operators discussed above based on DEBC. 
The BCFT analysis confirms that they are indeed the 
leading operators and the N fixed point is stable for g < 1 
and is unstable for g > 1. 

For the Neumann boundary condition on $, the con- 
ductance tensor vanishes: Gjk = 0. This is most easily 
seen by observing the Neumann boundary condition $ 
is nothing but the Neumann boundary condition on the 
original fields (^1,2,3 • Alternatively, we find that Gjk van- 
ishes by an explicit calculation of Eq. (jD23|l . 
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B. Andreev reflection fixed point - 
boundary condition on $ 



Dirichlet 



The other simple boundary condition for the free boson 
is the Dirichlet boundary condition on $, which is dual 
to Neumann. In fact, the Dirichlet boundary condition 
on $ is the Neumann boundary condition on O. The 
choice of angle ^ = tt pins $ = C. Plugging TZtt = — 1 
into Eq. 11U.21II . we obtain Ay — |6p. Using this, we 
find, for example, that 



^ 01 



45 I- r> |2 45 



9 = 



9 ■ 



(10.29) 
(10.30) 
(10.31) 



There are 3 independent operators with dimension 4g/3, 
3 with dimension g/3, and 3 with dimension g, plus their 
respective Hermitian conjugates. The leading irrelevant 
terms are the ± cycle tunneling operators. This fixed 
point is stable for g > 3. 

Within the BCFT approach, in the standard compact- 
ified boson case, the scaling dimensions of the bound- 
ary operators in the Neumann and Dirichlet boundary 
states are given by the reciprocal lattices as discussed 
in App. [ni If the present problem of the junction were 
mapped to the standard compactified free bosons, the 
Dirichlet boundary state would have the leading bound- 
ary operator with the scaling dimension 4(7/3, which is 
irrelevant for g > 3/4. Thus, the Dirichlet boundary 
condition would be stable for 3/4 < g. In particular, 
both the Neumann and Dirichlet boundary conditions 
are then stable for 3/4 < g < 1. This is discussed by 
Yi and Kane 16j in the context of the quantum Brow- 
nian motion on a triangular lattice. There is a duality 
g <-> 4/(35) which exchanges the Neumann and Dirichlet 
boundary conditions. 

These considerations for the standard compactified bo- 
son exactly matches the potential argument presented in 
Sec. As discussed there, it would not apply to the 
present problem in which electrons hop between differ- 
ent wires. In Sec. the fermi statistics of electrons 
were taken of care by the Klein factors, to derive sensi- 
ble results. In the present context, the difference from 
the standard boundary problem of free boson field the- 
ory is given by the twisted structure of the Hilbert space 
discussed in Section IVIIII This affects the stability of 
the fixed point in an important manner. As it will be 
shown below, the most stable Dirichlet-type boundary 
state turns out to be the strong pair-hopping limit. Dp, 
discussed in Sec. 

The most stable Dirichlet-type boundary state is given 
by a superposition of all the Ishibashi states allowed un- 
der AQ = 0, namely ni = 7^2 = 0. These Ishibashi states 



are given as 



|(0,0,TOi,TO2))) = exp 



E 



1(0,0, mi, 7712)), 



(10.32) 

As in the construction of the N boundary state, be- 
cause of the parity condition (|8.7t)|) . 777,1 and 777,2 take 
only even integer values. Thus we have 



9Dr 



E 

,m!^ £Z 



1(0, 0,2777;, 2777^))). (10.33) 



The corresponding diagonal cylinder amplitude is 



ZDpDp{q){Dp\e 
2 

E 

1^ ,m2tzZ 



-IH 



9Dr 



f \Dp) 

nT^k'^ra'lK^^lra'^Ki.f ^ (10.34) 



Modular transforming using Eq. (|C59|) . now with A 
being the triangular lattice with the lattice constant 2, 
we obtain 



Z 



DpDp 



1 



5(i/li?l + ii2fl2)^ 



(10.35) 



Now the dimension of the leading boundary operator is 
given by 



(10.36) 



Furthermore, the dimensions of the second-leading and 
third-leading operators are related to the distance be- 
tween next-nearest and next-next-nearest neighbors on 
the triangular lattice and given by g and 4(7/3, respec- 
tively. These exactly agree with the scaling dimensions of 
tunnehng operators H10.30() . H10.31|I and H10.29|l obtained 
from the DEBC method. Furthermore, these scaling di- 
mensions coincide with those that can be obtained from 
the instanton argument for the Dp fixed point, in which 
$ is pinned to the triangular lattice with the lattice con- 
stant 

The present analysis confirms that g/3 is the scaling 
dimension of the leading boundary operator, and thus the 
Dp fixed point is stable for 3 > 3 and unstable for g < 3. 
This is quite different from what would obtain for the 
stability of the Dirichlet boundary state if the problem 
were mapped to the standard compactified boson. The 
difference is due to the twisted structure, which originates 
from the fermionic nature of the electron. 

The present result can also be understood from the 
duality. Eq. H8.82|l shows that the twisted boson theory 
has the duality 



e 
3 

9 



(10.37) 
(10.38) 



instead of g ^ 4/(3(7) for the standard boson. The dual- 
ity (|10.38|) was found in Ref. by a different argument. 
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Our analysis shows that the duahty exists as a structure 
of the effective Hilbert space, beyond a particular choice 
of the perturbation. 

Let us now discuss the conductance. For the Dirich- 
let boundary condition, a similar calculation to that in 
Add. IdI gives 



lim 



1 



dx 



uj^o+ ttujL Jo 
We thus find, from Eq. (|D2l 



Gki = 2.g— - ^ Vkf_iVifj, 

^"=1,2 



(10.39) 



(10.40) 



The summation X];i=i 2 interpreted as an inner 

product of the two-dimensional vectors Vk and vi. From 
eq. HD19|) , we find 



Vk • Vl = {Ski - 



Therefore we recover the results in Ref. Q as 



Gki — 



3 h 



ik = l) 



(10.41) 



(10.42) 



This derivation and the result on the conductance equally 
apply to Dp and Dn fixed points; they are indistinguish- 
able in terms of the conductance. As pointed out in 
Ref. 0, the single-terminal conductance in this case is 
larger than the single ideal wire. This exhibits an en- 
hanced conductance due to the Andreev reflection. 



For example, pinning the phases in the 
sponds to the boundary condition 



X+: K,-<i>+\ -{zxK,)-e^C. 



cycle corre- 



(10.43) 



And pinning this combination corresponds to the choice 
of angle ^ such that 



tan - = — . 

2 9 



(10.44) 



With this choice of boundary condition or angle ^, the 
scaling dimension, for example, of a backscattering term 

Oand6= -2y^zxK3 



such as T^'" is obtained with a 



in Eq. H10.21|l . Using that 



(z X v) — (TZ^v) ■ V — \v\ cos^ for any vector v: 



and [TZ^{zv)] 



2^ 



|zxi?3|'-[7^^(^xi?3)]•(5xX3) 

2 r. .1 45 



2||i^3|^ [l-cose] = 3 



9' 



(10.45) 



after some simple algebra. 

Similarly, the scaling dimensions of the — cycle terms 
(such as T^z") as well as those of the LL — RR terms 
(such as Ti^) can be shown to be also 3^^- Again, 
these operators are not independent after the choice of 
boundary condition; in total, there are 3 independent 
operators and their 3 Hermitian conjugates. This agrees 
with the dimension of the leading operator for the fixed 
point x+, from the instanton analysis in Sec. I VII 

Now, had we chosen to pin the phases in the — cycle 
set of operators, or select the boundary condition 



C. Chiral fixed points x± 

1. Chiral fixed points from DEBC 

In the DHM representation discussed in Section lVTI the 
quantal phases due to the Fermi statistics and the mag- 
netic flux piercing through the junction was accounted 
by the "boundary magnetic field" . In the representation 
with n = 0, 1 the chiral fixed points x± correspond to 
the localized phase of the DHM, which amounts to the 
Dirichlet boundary condition on the field X. We empha- 
size that the field X is not identical to the fields Q and 
$ we have been using. The x± fixed point corresponds 
to neither Dirichlet nor Neumann boundary conditions 
in terms of and $. As we will see in the following, the 
field X in fact corresponds to a chiral rotation of 8 or $. 

Indeed, by inspecting the operators in the ± cycles, we 
would find that if these coupling constants dominate the 
low-energy physics (flow to strong coupling), then they 
set boundary conditions that are mixed in the and 
fields. 



X-: i^. 



- (z X if,) . e 



C, 



(10.46) 



we would have obtained an analogous case, in which 

tan I = — . The dimensions of the backscattering and 
2 g o 

-f- cycle terms, as well as of the LL — RR terms, are again 
g^^. This fixed point is stable for 1 < g < 3. Which of 
the x± is selected depends on the time-reversal symmetry 
breaking flux. 



2. Relation to the DHM representation 

Let us connect the present analysis to the ones based 
on the dissipative Hofstadter model, more explicitly. 

First we formulate the DHM (for y = 0) as a boundary 
problem of a c = 2 free boson field theory. This can be 
done in two different ways. One of them is to begin with 
the (1+1) dimensional action: 



S 



dT 



dx 



— {d^Xy + ILe^^-eabd^X^d.Xb 
47r 47r 



(10.47) 
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We do not impose any b.c. on X at a: = 0, instead 
letting the field be dynamical on the boundary. By in- 
tegrating out the bulk part of the free boson fields, X, 
this is reduced to the (non-interacting part of the) one- 
dimensional action of Eq. (|6.10(l introduced in Ref. [23| • 
To see this explicitly, we may decompose the fields 
X(t, x) inside the path integral into a sum of classical 
and quantum terms: 

X{t,x) = X''\t,x) + X{t,xY. (10.48) 

Here X'^^ obeys the classical equations of motion: 

d^X''^ = 0. (10.49) 

[Note that the term proportional to jS in the action of 
Eq. 1)10. 47|) makes no contribution to the Euler-Lagrange 
equations because it is a total derivative.] X'^^ is chosen 
to obey the boundary condition: 



X'\t,0) = Xbir). 



(10.50) 



On the other hand, the quantum part, X'{t, x) obeys a D 
b.c, X'{t,0) = 0. With this choice of b.c. on X' , using 
the fact that X"^^ obeys the Euler-Lagrange equations, we 
see that: 



S{X''^ + X') = S{X''^) + S{X'), 



(10.51) 



i.e. the cross-term vanishes. Furthermore, in calculating 
Green's functions of X{t,Q), X' doesn't appear due to 
its D b.c. Therefore, the calculation of these boundary 
Green's functions reduces to a functional integral over 
the boundary fields, Xi,{t) with a classical action: 



Sb[Xb] = S{X''') 



(10.52) 



where S is the two-dimensional action of Eq. (|10.47|l and 
the dependence on the boundary fields, Xtir) is induced 
through the b.c. on X^' of Eq. p0.50)l . Fourier expand- 
ing X^: 



gives the classical solution: 



X''\T,x)=J2^n 



-\uJn\x 



(10.54) 



and the boundary action H6.10|l . 

The alternative approach is to take the bulk action 
with no magnetic field term: 



dr / dx 



and impose a chiral boundary condition: 



Xr — TZ^'Xl 
for some rotation angle . 



(10.55) 



(10.56) 



We wish to show that correlation functions at x = for 
Eq. 1)10.55(1 with 1)10.56)1 are the same, for an appropriate 
value of as those obtained from the boundary action 
of Eq. ()6.10)) . This can be shown using a limiting pro- 
cedure where we first calculate correlation functions for 
fields X(Ti,Xi) and then take the limit where Xi . 
Since we are considering the non-interacting theory at 
this point, it is sufficient to show this for the 2-point 
function. Using: 



{X^{t + ix)Xl{T' + ix')) = ^ ln[(T - t') + i{x - x')] 

(10.57) 

and the fact that the boundary condition of Eq. ()10.56l) 
implies: 



Xr{t,x) ^TZ£_'Xl{t,~x), 



(10.58) 



X{r) 



Xr, 



(10.53) for a: > 0, we obtain: 



{x''{t,x)x\q, x')) - {[xUt + tx) + ni?xi{T - ix)][xl[zx') -f nfxti-ix')]) 



2a 



{6''^\n[T + i{x~x')]+6''^\n[T-i{x-x')]+nf\n[T + i{x + x')]+n\'f\iY[T-i{x + x')]}. (10.59) 



Using the explicit form of TZ^i given in Eq. 1)10.17)) this can be written: 



1 



(X''(t,x)X\Q,x')) = - — {[l + coseV''ln[r^ + (a;-a;') 
2a 

Now letting x, a;' — > 0+, we obtain: 

(X'^(T,+0)X^(0,-fO)) - 



e sm 



iin (In [t - i{x + x')]-\n[T + i{x + x')])} . (10.60) 



I + COSC' ,a. 



2a 



5"" luT^ + in 



sine' a 



2a 



e SI 



5n(r) - 1]. 



(10.61) 
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Choosing: 

tan^ = -, (10.62) 
2 a 

this reduces to the correlation function of the DHM 
model, Eq. H6.11|l . apart from a constant term: 

6{X^T, +0)X''{0, +0)) = -tn . ^ ^„ £°^ (10.63) 

This constant term has no effect on the perturbative ex- 
pansion discussed in Sec. IVII and used to connect the 
DHM with the Y-junction model. This follows since the 
extra contribution to Eq. (|6.14l) is proportional to 

eS:i=i 1:^=1 ^ 1. (10.64) 

Thus it is shown '2^ that imposing the chiral b.c. of Eq. 
(Il().56|l on 2-component bosons with the standard action 
of Eq. I10.55|l with the rotation angle, C of Eq. H10.62|) 
gives the DHM model. 

In fact, it is possible to derive a general connection be- 
tween the fields X{t, x) introduced here and the original 
fields $(r, x) introduced in Sec. IIVI Naively, by compar- 
ing the actions, we might expect that: 

X = y/gja^ ? (10.65) 

However, this does not take into account the b.c. cor- 
rectly. Recall that the disconnected wire fixed point cor- 
responded to N b.c. on 

$L = (10.66) 

However, the above discussion establishes that this N. 
V = 0, fixed point corresponds to the chiral b.c. of Eq. 
(|10.56|l on X. Thus we tentatively identify: 

Xr = ^/gJaU^'^R 

Xl = \fgl^^L- (10.67) 

To check the consistency of this identification, consider 
the chiral b.c. on $, given by Eq. 1)10. 16|l with the ro- 
tation angle of Eq. (|10.44|) . As we argued in Sec. IVII 
this fixed point corresponds to the localized phase of the 
DHM which in turn corresponds to a simple D b.c. on 
X: 

Xl = -Xr. (10.68) 

Consistency of Eq. (|10.67|l [with f ' given bv Eg. I|10.62|l ]. 
which maps X into <I>, with the b.c. on $ of Eq. (|l().16|l . 
(|10.44|l and the b.c. on X of Eq. H10.68|) requires: 

7^e = -7^^^ (10.69) 

or 

C' = 7r-^. (10.70) 

This implies tan(^'/2) = cot(^/2) and hence is consistent 
with Eqs. H10.62I) and H10.44|l for ^' and i together with 
the value of (3 /a determined in Eq. H6.3()|l for the case 
n = 1, corresponding to the stable chiral fixed point. 



3. Boundary states and the generalized chiral fixed points 

We have now established that the chiral fixed points x± 
corresponds to the chiral boundary condition (|10.16() on 
the field Let us construct the corresponding boundary 
state explicitly, taking the twisted structure discussed in 
Sec. IVIlil into account. 

The winding number along the boundary can be writ- 
ten as 




(10.71) 
(10.72) 



For the chiral boundary condition (|10.16(l . we have 

^A$_^Ae = 7^^(V5A$ + ^Ae) (10.73) 

V 5 \ 9 

For ^ 7^ 0, TT, it is clear that we need both A$ and A8 
nonvanishing. Moreover, TZ^ is a rotation matrix and 
thus preserves the length of the vector. This requires 
A$ and A8 to be mutually orthogonal. 

Thus the only oscillator vacua which can appear in the 
chiral boundary state is 

\{rh,rl2,sh,sl2)). (10.74) 

for arbitrary integers li and l2- Here (r, s) is a fixed set 
of integers. In order to satisfy the parity condition (|8.76(l 
for arbitrary integers Zi^2, 

r = s (mod 2), (10.75) 

is required. The chiral rotation angle is fixed by (r, s) as 

tan| = -^. (10.76) 

That is, the chiral rotation angle is "quantized" . 

For each vacuum H10.74() . we can construct a chiral 
bosonic Ishibashi state 

\{rli,rl2,sli,sl2))) = exp(^a^''' ■ TZ^d^'')\{rli,rl2, sh, si 

n=l 

(10.77) 

which is conformally invariant. 

A single Ishibashi state does not satisfy Cardy's con- 
sistency condition. Similarly to the case of the standard 
Neumann/Dirichlet boundary state, the summation over 
all integers li^2 could give a consistent boundary state. 
Namely, 

lx>=5x E \{rlurl2, slush))), (10.78) 

where is again the ground-state degeneracy |^ of the 
chiral boundary states. 
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The diagonal cylinder amplitude is given by 

2 



~i ['leys'- 4^ + 75 Jfi)+'2(y9riil + ^K2)J 



This can be written as 



E 



(10.79) 
(10.80) 



where A* is a triangular lattice with lattice constant 



Thus, choosing 



9x 



4 V 3 9 



(10.81) 



(10.82) 



the same amplitude in the open string channel is given 
by 



Z 



1 



XX 



\ ^ w 



(10.83) 



because of eq. (|C60|I . The lattice is dual to A* , and 
is the triangular lattice with the lattice constant 



4g 



(10.84) 



This formula encodes, as usual in the BCFT approach, 
the complete list of scaling dimensions of possible bound- 
ary perturbations to the chiral fixed point. Namely, the 
scaling dimensions of all the boundary perturbations are 
given by W"^ + integer. In particular, the scaling di- 
mension of the leading perturbation corresponding to the 
minimal length basis vector of A^, is given by 



4.9 



3j,2 ^ g2g2 



(10.85) 



From this spectrum of boundary operators, and the 
ground-state degeneracy H10.82|l . we see that the bound- 
ary state with larger r and s are unstable. While any 
given chiral rotation angle ^ may be approximated with 
arbitrary precision by taking large r and s, such a choice 
gives a highly unstable fixed point. 

As a special series of chiral boundary states, let us take 
r = 1. Then, due to the parity condition H10.75|l . s must 
be odd: s = 2n — 1 for an integer n. Now the boundary 
states are characterized by the chiral rotation angle 



tan - = {2n- 1)-^. 
2 ^ ' ^ 



(10.86) 



This agrees with the chiral rotation angle for the chiral 
fixed points obtained from the mapping to the DHM. 



Indeed, the scaling dimension of the leading perturbation 
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3+ (2n- 1)V 



(10.87) 



also agrees exactly with the instanton analysis on the 
chiral fixed points in the DHM approach. In particular, 
n = 1 , corresponds to the most stable ones which we 
denoted as x±- 

We have thus succeeded in explicitly constructing the 
boundary state for the chiral fixed points. In general, 
a chiral fixed point is characterized by two non-zero in- 
tegers (r, s) with r — s (mod 2). It includes the one- 
parameter series r = 1 , s = 2n — 1 which was found in 
Sec. IVII However, except for (r, s) = (1,±1), the chiral 
fixed point is always unstable for any given value of g. 

The only chiral fixed points that can be stable are 
(r, s) = (1, ±1), that is x±. The BCFT gives us the com- 
plete spectrum of the boundary operator, and thus more 
convincing evidence of the stability. In particular, for 
1 < g < 3, x± fixed points should be stable against any 
boundary perturbation, including the asymmetry and the 
change of the fiux 0, but excluding the "resonant" type 
perturbation as discussed in Sec. lVIIl'Bl 



4. Conductance 

Let us calculate the conductance at the chiral fixed 
points x± from the above analysis on the boundary con- 
dition. Using Jj^ = T^^L^ ^^^^ 

lim I dx / dT e"^'' {J^'{v,T)r{x,Q)) 



1 

= g[(l~cosO'5'''' + sin^e^,] 
Thus, from Eq. (|D23p we have 

„2 



(10.88) 



Gjk = 5-^ [(1 - cos^)vj ■ Vk + sin^Vj x -i^] (10.89) 



(1 - cos ^){Sjk 



1 , sin £ 

3) + ^ej4lO-90) 



For the chiral rotation angle £ in Eq. I|10.76|l . this reduces 
to the Z3 symmetric form (12.511 . with 



Gs = 



Ga 



Ags^ 



h g-^r-^ + 3s 
g2 ig'^rs 



(10.91) 
(10.92) 



h g^r^ -\- 3s^ 

In particular, for the most stable chiral fixed points x±i 
where r = 1 and s = ±1, we recover the same re- 
sult (|5.92|) derived from the mapping to the DHM. 



D. Asymmetric fixed point 

In this paper, we have mostly discussed the junction 
with Z3 symmetry among the wires. However, if the 
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junction of three wires could be realized, it will inevitably 
contain some asymmetry between the wires. Therefore, 
it is important to study the effect of Z3 asymmetry. 

Here we restrict our study to the asymmetry at the 
junction, assuming the three quantum wires are identical 
in the bulk. Although there would also be some asym- 
metry in the bulk in reality, it is outside the scope of the 
present paper. 

There is a natural fixed point corresponding to the 
limit of the strong anisotropy: two of the wires are 
strongly coupled to form a single wire with perfect trans- 
mission at the junction, and the other is decoupled. 
Without losing generality, let us assume that the wires 1 
and 2 are strongly coupled and the wire 3 is decoupled. 
Let us discuss the stability of this fixed point, which we 
call Da. 



1. Heuristic analysis 

As perturbations to this asymmetric fixed point, we 
consider the following three operators: 

electron hopping between wires 3 and !& 2 
tplips + h.c. 

electron pair hopping between wires 3 and 1&: 2 

V'lV'iV'aV's + h.c. 

electron backscattering in the wire 1&: 2 



h.c. 



It is natural to expect that the leading perturbation with 
the smallest scaling dimension is given by one of these 
three. Indeed, we confirm this expectation later with the 
BCFT analysis. 

In the first two cases, the scaling dimension of the op- 
erator is given by the sum of the scaling dimensions of 
the operators on wire 3 and wire 1& 2. The discon- 
nected wire 3 is subject to the Neumann boundary con- 
dition on (Dirichlet on ^3) at the junction. Thus, 
from eq. H4.3() , the electron annihilation operator ip^ cor- 
responds to the boundary operator e^'^^l^ at the Dirich- 
let boundary condition on ^3. Its scaling dimension turns 
out to be 1/(25). 

On the other hand, for the wire 1& 2, the scaling di- 
mension of the electron creation operator V'l is given by 
the bulk scaling dimension of V'li, or V'ih- It turns out 
to be (5 + 5^^)/4. This scaling dimension determines the 
Local Density Of States for the electron tunneling into 
the TLL "sl. 

As a result, the scaling dimension of the electron hop- 
ping operator + h.c. is 



1 

2^ 



(10.93) 



This, remarkably, is exactly the inverse of the scaling 
dimension (|10.85|l of the leading perturbation at x± fixed 



points. As a consequence, it is relevant for 1 < 5 < 3 and 
irrelevant for t; < 1 or 3 < g. 

Next we consider the electron pair hopping between 
the wire 3 and 1& 2. The pair annihilation operator 
in the wire 3 ?/'3'!/'3 is given by the boundary operator 
gi\/2ip3/ ^]^g Djj-jciiigi; boundary condition on ^^3, with 
the scaling dimension Ijg. The leading operator in the 



pair creation operator "01 "01 is given by 



(10.94) 



Although this has the identical form to the above ipsips, 
its scaling dimension is different since it must be evalu- 
ated as a bulk operator. The scaling dimension is given 
by 



^ 1 2 1 
2 X - X - = -. 
4 g <7 



(10.95) 



Thus, the scaling dimension of the pair hopping operator 
is 



2 13 

g 9 9' 



(10.96) 



Finally, the scaling dimension of the backscatteringop- 
erator in the wire 1& 2 is known to be g from Ref. 



2. BCFT analysis 

Now let us confirm the boundary operator content with 
the BCFT. The Da fixed point corresponds to Dirichlet 
boundary condition on $1 = {ipi—ip2)/^/2, and Neumann 
boundary condition on $2- 

Thus, the allowed winding numbers in the oscillator 
vacua is 

A$ = ~27m-^ = 2?™—!-^ — (10.97) 

AO ^ -27rTOi?3 = 27rm(i?i +X2), (10.98) 

where n, m are integers satisfying 

n = m (mod 2), (10.99) 

from Eq. H8.76|l . The unit vectors Ki and Ri are defined 
in Eq. H4.10|l and (|6.55|l respectively. 

For each permitted oscillator vacuum | (n, n, to, to)) , we 
can construct the bosonic Ishibashi state 



(71, 71, TO, to))) = exp (^af ^ • 7^(;a^^)|(n,n,TO,m)), 

(10.100) 



k=l 



where the rotation angle Q is given by 
C TO Vs 



tan ■ 



n g 



(10.101) 



similarly to H1U.44|I . 
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The boundary state for the Da fixed point is given by 
the superposition of these Ishibashi states over all possi- 
ble n, m: 

\Da)^9d^ |(n,n,m,m))). (10.102) 

n—m (mod 2) 

Although the construction is similar to that of the chiral 
boundary state H10.78|l . \Da) is not a chiral boundary 
state. Since the chiral rotation angle C determined by 
Eq. IjlO. 10111 depends on each Ishibashi state, the bound- 
ary state as a whole cannot be simply related to a chiral 
boundary condition as in (|10.16|l . Nevertheless, it can be 
a legitimate conformally invariant boundary state. 

The diagonal cylinder amplitude in the closed string 
channel is given by 



Z 



DaDa 



9Da 



9Da^ 



E 1' 

-m (mod 2) 



(10.103) 



where is the two-dimensional Bravais lattice spanned 
by {l/^,±y/gji). Using Eq. fUHUIl . we fix 



1/4 



to obtain 



ZdaDa 



1 



(10.104) 



(10.105) 



weAj 



similar manifold corresponding to the S-matrices of the 
free dual fermion, at g = 3. For g < l, the backscattering 
in the wire 2&; 3 is relevant. It is most natural to assume 
that the system flows into the N fixed point, in which all 
the three wires are disconnected. 

For 1 < g < 3, the electron hopping is relevant. In 
the time-reversal non-invariant case 7^ 0, tt, it is most 
natural to assume that the system flows into the x± fixed 
point depending on the sign of the fiux (j>. When = 0, tt, 
the system is time-reversal invariant and thus cannot flow 
into x±- We expect the infrared fixed point to be the 
nontrivial M(ysterious) fixed point, which is stable for 
1 < (7 < 3 in the time-reversal invariant models. 

For 3 < g, the electron-pair hopping is dominant. In 
the time-reversal non-invariant case 7^ 0, tt, it is most 
natural to assume that the system flows into the Dp flxed 
point, which can be regarded as the limit of strong pair 
tunneling. In the time-reversal invariant cases cf) — 0,tt, 
the situation is subtle. As discussed in Sec. El the simple 
strong pair hopping limit in the time-reversal invariant 
case suggests that the potential minima forms a flner lat- 
tice. This flxed point would correspond to Dn boundary 
state and thus is unstable for g < 9. If this argument ap- 
plies to any time-reversal symmetric case, for 3 < g < 9 
we expect the system to flow into a nontrivial fixed point. 

In any case, our analyses imply that the system re- 
covers the Z3 symmetry, even if the microscopic model 
of the junction is not Z3 symmetric. This is remarkable 
compared to many other quantum impurity problems, in 
which a high symmetry is required to reach a nontrivial 
fixed point. This feature is favorable in observing the in- 
triguing behavior of the chiral fixed points x± in future 
experiments. 



where A^ is the Bravais lattice spanned by 
(^/2, ±v/3/(45)) and is dual of A^. 

The scaling dimensions of the boundary operators can 
be read off from the amplitude, as for W G Aa- 
These include 




9 



45 



(2^,0)^ = g. 



(0,2,& = ^, 
V 45 g 



(10.106) 
(10.107) 
(10.108) 



which exactly match the three operators considered in 
the heuristic analysis. 



3. Restoration of the Zz symmetry 

We have shown that the asymmetric fixed point Da is 
always unstable except for g = 1, 3. At g = 1, the flxed 
point belongs to the continuous manifold of flxed points 
described by free fermion S'-matrices. Thus it is natural 
to have marginal operators. Presumably it belongs to the 



XI. CONSTRAINTS ON CONDUCTANCE 
FROM ENERGY CONSERVATION 

In this section we will show how energy conserva- 
tion constraints the conductance tensor. We will also 
show that the three sets of boundary conditions on the 
bosonic flelds studied above, namely Neumann, chiral, 
and Dirichlet, correspond to a situation where all the en- 
ergy exchanged in the scattering processes off the junc- 
tion is between the zero modes. Drops or increases of 
voltages across the junction are associated with the zero 
modes, or current carrying states. The non-zero or os- 
cillator modes are neutral, and can carry energy but no 
electric current. Therefore, scattering without exciting 
the oscillator modes is the most efflcient way of trans- 
mitting current without additional dissipation. 

It is simpler to discuss the energy exchanges in the 
unfolded formalism. Let (j)^'{t,x), for i — 1,2,3, be the 
bosonic flelds in each wire, 'pf'it, x < 0) correspond to the 
incoming flelds into the junction at a: = 0, and <j>f'{t, x > 
0) are the outgoing flelds after scattering. 4)ii{x) obeys 
periodic boundary conditions. 



(I)r{x + 21) = (j)ii{x) + 2Tm. 



(11.1) 
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Away from the impurity, the dynamics is governed by 



(11.2) 



The mode expansion of the fields (l)f'{t,x), in a circle of 
length 21, reads: 



(11.3) 



1 ^ 1 



E 



2g ^ x/n 

n— 1 



-mfit-x) _|_ 



and the Hamiltonian can be written as 

TT " 1 ^ 



n=l 



(11.4) 



where the operator Qf takes integer eigenvalues Qf. 
(Notice that the Hamiltonian and the mode expansion 
for each of these three chiral fields are exactly those in 
Eqs. H8.5|l and (|8.7|l . respectively, but with /3 replaced by 
I and g ^ 1.) 

The scattering picture is as follows. If an incident pulse 
of some size w is prepared and imparted into the scat- 
terer at the origin a; = 0, it will propagate freely ac- 
cording Eq. I|11.2|l as long as the pulse does not overlap 
with the origin, i.e., its amplitude is vanishingly small at 
X = 0. In this case, the energy of this incoming pulse 
is obtained from the Hamiltonian Eq. H11.4|l . and it is 
independent of the boundary interaction at the origin, 
which is what determines the scattering at the junction. 
Past the junction, once the outgoing pulse has reached 
a vanishing overlap with the origin, its energy is again 
obtained from the Hamiltonian Eq. H11.4|l . 

Now, in the case of no tunneling between wires, the 
dynamics described in Eq. (|11.2I) is exact even at a; = 0, 
because the no tunneling regime corresponds to free prop- 
agation through the junction; this is the case of perfect 
backscattering, or N BC in the unfolded formalism [in 
which case the Hamiltonian in the circle of size 21 in 
Eq. H11.4|l should coincide with that for a finite segment 
of size I in Eq. ljC46p ]. We argue that, once tunneling 
is turned on, the Lagrangian should change due to the 
tunneling terms at x = 0, but that away from a; = 
one can still use the free Lagrangian Eq. (|11.2|) as long 
as the scattered pulse goes far enough from a; = so as 
to have vanishing overlap with the origin. Finally, we 
should consider the situation where the scattered pulse 
does not go all around the circle and scatters again; such 
scattering picture is possible in the limit of Z — > oo, so 
that scattered charges do not return to take the place of 
our prepared initial state. 

Let us prepare an incoming state |IN) consisting of only 
zero modes being occupied; their occupation is dictated 
by the voltages Vi applied to the leads. After scattering 
off the impurity at a; = 0, the state becomes |OUT), 
which in general has excitations created in both the zero 



modes and the oscillator modes (n > 0). The outgoing 
currents are solely determined by the excitations of the 
zero modes. 

Let I™ and for i = 1, 2, 3, be incoming and outgo- 
ing currents in each of the three leads, respectively. The 
incoming currents are related to the voltages on the three 
leads through 



g2 g2 



(11.5) 



whereas the outgoing currents are obtained from the 
junction conductances Gtj as in Eq. (1221) through 



E^ 



(11.6) 



where li is the net current flowing into the junction from 
lead i. 

The quantum number is the total charge of the 
state, which in the unfolded formalism is proportional to 
the total current J]" or /°"* for the |IN) or |OUT) states, 
respectively. (The currents are given by the charges mul- 
tiplied by the velocity, which we set to unity, divided by 
a length /). Define 



AEn 



E 



ttI r 



(11.7) 



This is the difference between the energy carried by the 
zero modes of the two quantum states |IN) and |OUT). 
The incoming state has no excited oscillator modes, but 
the outgoing state may have them. Energy is conserved 
in the scattering process, so it follows [see Eq. (|11.4f) ] that 



AEo < 



(11.8) 



with the equality holding when no oscillator modes are 
excited in [OUT). 

Using the expressions for the conductances 
Eqs. (|11.5I11.6I) . and substituting in Eq. (|11.7() . we 
obtain 



gAE 
Itt 



ijk 



(11.9) 



1 



2 ^ 

ijk 



^2 



^ ' GijGik - 2g^ Gij Sik ] {Vj - VkY , 



where we have used that Gij = 0, which follows from 
current conservation, in order to express the second line 
in Eq. (|11.9I) in terms of voltage differences. 

Now, the most general conductance tensor that satis- 
fies the Z3 symmetry of the three lead problem is given 
by Eq. H2.5|l . Recall that a non-zero Ga is allowed if 
time-reversal symmetry is absent. Defining dimension- 
less parameters 



, . Ig^A (11.10) 



G = Gs^ and A 
e 
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FIG. 18: Physical values of G and A, allowed from energy 
conservation, must lie inside the circle of radius 2(//3 and 
origin at {2g/3, 0) in the G — A plane. The boundary of the 
circle corresponds to the case when no energy is transfered to 
oscillator modes upon scattering off the junction {AE — 0). 
The points corresponding to boundary conditions N, Dp , and 
X± lie on this boundary, and are are marked on the figure. 



for convenience, and inserting the expression for Gij in 
Eq. (|2.5|l into the second Hne of Eq. (|11.9(l . we obtain 
after some elementary manipulations that 



gAE 



3e^ 
8h^ 



G^.I^G + A^ 
3 



jk 



(11-11) 

The last term is always positive, so in order to satisfy 
AE < one must have that 



- ^ G + < , 



G 



A^< 



(11.12) 
(11.13) 



which imposes constrains to the allowed values of G and 
A. The physical values are restricted to inside a circle of 
radius 2g/3 and origin at {2g/3, 0) in the G — A plane, as 
shown in Fig. 1181 The boundary of the circle correspond 
to the cases where AE = when no energy is transfered 
to oscillator modes. 



T-invariant case: A = 



If time-reversal is a symmetry in the problem, then A 
must vanish, so that Gij = Gji. In this case, the allowed 
physical values of G lie in the segment that spans the 
diameter of the circle in the A = axis, or 



The equality AE = occurs when G = and when 
G = 4(jr/3; these cases correspond, respectively, to N and 
Dp boundary conditions (for a summary of the conduc- 
tances for different boundary conditions, see section ITl|l . 
Hence, for these two boundary conditions, no energy is 
transfered into oscillator modes, and it is simply redis- 
tributed among the zero modes of the three boson fields. 



B. T-broken case: A / 

This is the general case, in which the physical values 
of G and A must lie within the circle, and the boundary 
corresponds to the cases AE = 0. There are two par- 
ticular cases, as shown in Fig. 1181 which correspond to 
the two chiral fixed points x± j when the angle A satisfies 
tan A = ±g/V3. In this case, the conductances are given 

by Gs = ^^ and Ga = = ±9Gs. 



XII. Y-JUNCTIONS ATTACHED TO FERMI 
LIQUID LEADS 

In this Section, we discuss the effective conductance of 
the junction in the presence of Fermi liquid (FL) leads, 
as introduced in Fig. 13 The FL leads can be thought of 
.9icad = 1 wires connected at the endpoints of the three 
Y-junction wires with g parameter. If the leads and the 
three quantum wires are connected by a point contact, 
then one can show that only one mode of the FL couples 
to the wire. This is similar to the case of the Kondo prob- 
lem, where only the s-channel couples to the point-like 
impurity. In the case of the FL-quantum wire coupling 
at a point, it is not the s-channel but instead some other 
channel that couples, [s^- 

When we derived the conductance of the junction, we 
assumed that the wires were subject to voltages Vi, i = 
1,2,3. If we have Fermi liquid leads, what happens is 
that the reservoirs are at voltages Vi, i = 1, 2, 3, and one 
must obtain the relation between the Vi and Vi. If that is 
done, all we have left to do is to re-express the currents li 
in terms of the Vi. In other words, we need to calculate 
the conductance Gij from the Gij we already computed: 



The question is how to relate the Vi to the Vi. The answer 
is that the voltage difference V — Vi across the lead/ wire 
contact is related to the current li that flows through the 
contact: 



- G, iy, ~ Vi) 



(12.2) 



< Gs < 



45 e2 
3 h 



(11.14) 



where the fixed point conductance for 17 > 1 is Gc = 
2g/{g — 1) /h. This can be derived by looking at the 
problem of a single wire connected to FL leads, or in 
general for two TLLs with gi > 32, as done in Ref. |38j . 
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in which case = (g^^ — ^)/2 e^/h. For the case 
g < 1, the fixed point conductance is zero, and there is no 
transport between the leads and the wires. Notice that 
ii g —^ 1 then Gc — *■ c», and hence Vi = Vi for any finite 
current; this makes sense, as for if both sides have g — 1, 
there is no voltage drop across the lead/wire interface. 

It is easy to get Gy now. We express Vi in terms of Vi 
and If. 

V,^V~ h/G, , (12.3) 

and then substitute in the formula for the current in 
terms of the Vi : 

fc k 

or in matrix notation 

I ^GV -Q-^G I , (12.5) 
from which wc obtain that 

I ={t + G-^G)-^ GV = GV (12.6) 

or, finally, 

G^{t + G-^G)-^ G , (12.7) 
and more simply 

G-^ ^G-^ +G-^t , (12.8) 

which can be interpreted physically as saying that the 
resistance tensors of the Y-junction must be added to 
the resistances due to the lead/wire interface resistances. 
[Notice, though, that this is a formal expression since in 
practice the matrix G is singular because of the condi- 
tion 'Yl,k^jk~^- hence one must use Eq. (|12.7(l instead 
of Eq. H12.8|l .] We would like to point out that our ex- 
pressions Eqs. (112. 7112. 8|l lead to consistent results with 
those of Refs. |^ 122] when applied to the simpler two- 
wire problem. 

As an example, let us consider the effective conduc- 
tance at the chiral fixed points x±j in the presence of 
the FL leads. All we have to do is use the expression for 
the conductance tensor G and substitute into Eq. H12.7|l . 
For the chiral fixed points x±, recall the conductance 

Eq. insu- 

Substitution of Eq. into Eq. itTTTjl leads to 

G% = ^^m,k-l)±lejk]/2, (12.9) 

which is simply what one would obtain from the con- 
ductance tensor Eq. H6.92|l if we set g = 1. This re- 
sult is in agreement with numerical results that find that 
the current vanishes in one of the leads despite interac- 
tions, if the wires are connected to long non-interacting 
segments |l9j . The scaling exponent for the correc- 
tions, however, should still depend on g, and is given 



by A = 4g/(3 -|- 5^), which is also consistent with the 
numerical findings obtained for small attractive interac- 
tions. The sensitivity to the flux (and flow to x±) should 
also be a property of the interacting junction, not the FL 
leads. 

The effective conductance Gjk at other fixed points 
can be obtained in similar manner. The results are sum- 
marized in Sec. |nj We find that the effective conduc- 
tance tensor at each known fixed points can be identified 
with a "bare" conductance tensor of the junction of three 
Fermi liquids (</ = 1) at an appropriate fixed point. This 
is rather natural, because the system with finite (inter- 
acting) wires and infinite FL leads would eventually be 
renormalized into a junction of three FL leads in the low- 
energy limit. 



XIII. OPEN PROBLEMS 

In this paper we studied the problem of a junction of 
three quantum wires connected by a ring, through which 
a magnetic flux can be applied. The electron systems in 
the three leads were described in terms of a Tomonaga- 
Luttinger liquid of spinless particles, with dimcnsionless 
interaction parameter g. One of the main difficulties in 
studying the problem with three wires is the necessity 
to include Klein factors that ensure the proper ferniionic 
statistics for electrons in different leads. 

We approached the Y-junction problem using different 
methods: mapping to the dissipative Hofstadter model, 
boundary conformal field theory, and delayed evaluation 
of boundary conditions. These methods gave consistent 
results and allowed us to calculate the low-voltage and 
low-temperature fixed point conductance tensor for the 
junction as a function of the interaction parameter g. 
(See section^]for a summary of the results in the paper.) 
In this concluding section, we would like to list a few of 
the open directions that require further investigation. 

The first one is the inclusion of electron spin in the 
problem. Even in the case of tunneling between just two 
wires, the inclusion of the spin degree of freedom leads to 
a rather rich phase diagram for the charge and spin con- 
ductances as a function of the interaction parameter Q]. 
For example, one finds situations where the spin conduc- 
tance can vanish while the charge conductance does not, 
and vice versa, or systems that conduct both charge and 
spin, or that conduct neither. In the case of the three- 
lead junction the situation could, in principle, be much 
richer. One might wonder, for instance, whether there 
can be fixed points where the junction can separate the 
spin and charge of an electron incoming from one lead 
into different outgoing leads. 

It appears that the method of delayed evaluation of 
boundary conditions may lend itself naturally to the 
study of the problem with electron spin. The reason is 
that this approach seems to be generalizable to multiple 
species. In the case of three leads with spin, one could 
try to extend the treatment laid down in this paper by 
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taking six branches altogether. 

One of the major open problems concerning the Y- 
junctions is the nature of the stable fixed point in the 
range 1 < g < 3 when time-reversal symmetry is pre- 
served {(f) — 0, tt). So far, we have not been able to iden- 
tify the correct boundary condition for this "mysterious" 
(M) fixed point. We can, however, be sure that the na- 
ture of this M fixed point is rather different from N, D 
and x±- For example, using the energy conservation ar- 
guments presented in the scattering approach of section 
IXII we can argue that energy is dissipated into oscillator 
modes upon scattering off the junctions, which was not 
the case in the other four fixed points N, D and x± ■ We 
expect that the M fixed point should also describe the 
infrared stable behavior of the resonant tunneling model 
of Nayak et al. 5] . 

It is also clearly important to test our results by nu- 
merical methods, which could also guide analytical work 
in understanding the M fixed point. One numerical 
method, used by Barnabe-Theriault et al. 19], is an 
approximation based on the functional renormalization 
group approach. Within this method, the chiral fixed 
points we predicted were observed for small attractive 
interactions, when time-reversal symmetry was broken 
by an external field. Another fixed point, with maximal 
symmetric conductance for non-interacting fermions, was 
found for the case of zero flux. It is not clear whether 
this maximal conductance value is related to the con- 
ductance of the M fixed point once the renormalization 
effect of Fermi liquid leads is taken into account (the in- 
teraction in the wires in the numerical studies was turned 
off far away from the junction), or whether the approxi- 
mation scheme can only capture conductance values that 
are equivalent to those of single particles scattering of 
a renormalized barrier. It would be interesting to push 
this approach towards larger values of the attractive po- 
tential {g > 3), when the D fixed point becomes stable, 
since the physics in this case will be that of pair tunnel- 
ing, and hence values for the conductance greater than 
the maximal one for single particle scattering should be 
attainable (seem for a discussion). 

Finally, the presence of the essential Klein factors in 
the Y-junction is related to a Fermi minus sign problem 
which would presumably make Monte Carlo methods dif- 
ficult. This is to be contrasted with the two wire case 
where the Klein factors can be trivially eliminated and 
the resulting Coulomb gas can be very effectively studied 
by classical Monte Carlo |43|. The Y-junction problem 
is perhaps the simplest system with a fermion sign prob- 
lem, which comes only from the junction. (In a lattice 
formulation, the problem corresponds to three ID sys- 
tems connected together by only one rung at one of the 
end points.) Hence, it seems natural that this system. 



together with some of the analytical results that we de- 
rive, may be used as a benchmark for approximations 
that aspire to solve the fermion sign problem. 
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APPENDIX A: FREE FERMION CASE 

In this appendix we calculate the conductivity for the 
free fermion model {g = 1). This is done by first cal- 
culating the 5'-matrix and then using the Landauer for- 
malism. We consider both lattice and continuum free 
fermion models. 



1. Tight binding model 

We first consider a tight-binding version of the con- 
tinuum model that we study using bosonization. ipnj 
annihilates a fermion on site n on wire j. Here n — 0,1, 
2, ...oo and j = 1, 2, 3, with the wire index j = identi- 
fied with j ^ 3. The Hamiltonian is (setting the hopping 
amplitude t ~ 1): 

oo 3 

H = -^^(^l_^.V„+i,j+/i.c.) 

n=0 j = l 
3 

-(f/2)^[e^*/3^J^^.^o,,-i+/z.c.]. (Al) 

Writing a single-electron state as: 

|vI/)=^$„_,V'L|0), (A2) 

where $„ j is the lattice wave-function, it is seen to obey 
the lattice Schrodinger equation: 
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-E$„j = -[<^n+l,j + ^n-l,j] [n > 1) 

E<^o., = -$i,j-(f/2)[e^'^/3$o,,-i+e-'^/'$Oj+i]. (A3) 

The scattering solutions of this equation take the form: 

<i>n.j = v4in,, e-''^" + ylout,, e^'=", (A4) 

for aU n> Q and j — 1, 2, 3. The corresponding energy, which follows from the first line of the Schrodinger equation 
Eq. imi) is: 

E = -2cosk. (A5) 

Noting that Eqs. (|A4|) and (jA5|) imply 

$1,, + E^o,, = -(An,, e'" + Aut,j e-^'=), (A6) 
we see that the second line of the Schroedinger equation Eq. (|A3p gives 

(f/2) [e^"*/3 {Ain,j-i + Aoutj-i) + e-''^/3(^.^,,,+i + Aut,,+i)] = An,, e'" + A^utj e^'^^ . (A7) 
Assembling the three components of the amplitudes, Ain.j and Aout,j into vectors. An and ^out7 this becomes: 

[e-^k _ (p/2) M] lout = - [e*'^' - (f/2) M] A, , ( A8) 

I 



where 

M = e"*''^B-^ + e-'"^'^B, (A9) 

and 

/O 1 0\ 
S = 1 . (AlO) 

VI o; 

Thus, the outgoing and incoming amplitudes are related 
by a 3-dimensional S-matrix: 

Xut = 5 An, (All) 

where: 

S = -[e-'^ - {fl2)M\-^[e'^ - (f/2) Af]. (A12) 
or, equivalently, 

5* = -e'"^^ [1 - e''=(f /2) M]-\l - e-'^{V/2) M]. (A13) 
I 



Note that we have taken $^(0) to be the average of its 
values at 0+ and 0~. This can be rewritten: 

[(r + 2i) + TM] lout - -[(r - 2i) + TM] 4n, (A18) 



2. Continuum model 



We use the continuum model with the Hamiltonian 
density: 



-5(x) {r[e''^/3^]^j-i + h.c] + rV-jV,} ■ 
j 

Here we have "unfolded" the system so that all fields are 
left-movers. We do not explicitly write the subscripts L, 
and we have set = 1- The corresponding Schrodinger 
equation is: 



(A15) 



(A16) 



(A17) 



where M is the same matrix defined in the discussion of 
the tight-binding model, Eq. HA9|I . Thus the S-matrix 



Writing the scattering states as 

$ / An,, e-f- (x>0) 

' ~ \ Aut,j e-^f- (x<0) 

with energy E — p, we obtain (for i? — > 0): 



i(An,i - Aut,i) = ^ {r[e*'^/3 (An.j-l + Aout,j-l) + e'"^/^ (An,j + 1 + ^out.j + l)] + r {Ain.j + Aut,j)} ■ 
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for the continuum model is: 



The S-matrix 



S = -[(r + 2i) + rM]-^[{r - 2i) + TM]. 



(A19) 



This will be the same S-matrix as in the tight-binding 
model [Eq. (|A12|l ] , for a particular value of fc, if we choose 
r and F so that: 



2i 



2i\ 



(A20) 



and 



An explicit expression for the S-matrix of Eq. (|A13|) 
can be found by straightforward algebra (we will drop 
the overall phase factor e*^*^ which is of no physical con- 
sequence to the conductances). This S-matrix can be 
parameterized by the most general symmetric form: 



S = So + S-B + S+B- 



(A22) 



r/2 = 



\r-2i\ 



(A21) 



It is convenient to define: 



It can be shown that the continuum limit of the tight- 
binding model, with the spectrum linearized around kp^ 
gives values of T and r satisfying this equation when 
k — kp. Note, in particular, that for the particle-hole 
symmetric case, kp = 7r/2, r = and F = F. 



z± = (f/2)e 



i(fc±0/3) 



(A23) 



The coefficients Si in the S-matrix then are given by: 



where: 



So = {z+z^ + \z+f + \z^f + z^_zl + zlz*_-l)/D 
S+ = {z*_ ~ z+ — z"^ + z^z*j_) / D 

S- = [zl- z^ - z\ + z+z*_)/D, (A24) 



D = 1- 3z+z_ - - z\. (A25) 



It can be readily seen that S —e^^^I at f ^ and S —I at f ^ 00. For f /2 = 1 a perfectly chiral S-matrix 
occurs, for some value of (/>, with 5*+ or S- = 0. In this case \z± \ = 1. Thus we see that: 

z^z+S+ = {z+ - z^z\ - z^_z+ + z'i)/D = (1 - z+z-)[z+ + z'i)/D. (A26) 

Thus 5+ = when z'iz^ = e'f''-'^) ^ -1 or 

fc-0= (2n+l)7r, (A27) 

for integer n. At this point, 

Sa = {z+z^ + l + z'izX + z\z*_)/D = {), (A28) 

and therefore only 5_ = —{z*_)'^ 7^ 0. This corresponds to an electron incident on wire j being transmitted with 
probability one to wire j — 1. Conversely, S*- = 5o = when F/2 = 1 and 

fc + 0= (2n+l)7r, (A29) 

so that an electron incident on wire j is transmitted with probability one to wire j + 1. 

In the time-reversal invariant cases, = or tt, S+ = S^; an electron incident on wire j has equal amplitude to 
be transmitted to wire j — 1 or j ' + 1. The maximum possible value of |S'+| = jS*-! consistent with unitarity of the 
S-matrix is \S+\ = 2/3, j^ol — 1/3. For = 0, one obtains 

- (l4("l-2z) ' ^^''^ 
with z = (r/2)e''^. For any value of fc, \S+\ reaches 2/3 for some value of F of order one. 
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4. Conductance 



In the free fermion case, we obtain the conductance directly from the S-matrix by the Landauer formahsm. We 
assume a thermal distribution of electrons heading towards the junction from distant reservoirs with different chemical 
potentials. The total current on wire j is this incoming current minus the reflected current minus the current 
transmitted from the other 2 wires. Thus: 




1A311 



where np{e) is the Fermi distribution function. Using 
V — (de/dk)/h and expanding to linear order in the V^'s 
gives at T = 0: 

= ^{[1 - is^m - \s,,,+i\%+i - 

(A32) 

where the S-matrix is evaluated at k = kp- Thus the 
conductance tensor is: 



(35,fc - 1) 



1^4 



15- 



2 2 

(A33) 

where S± are given by Eq. (|A24p . In particular, for the 
perfectly chiral S-matrices discussed above, this reduces 
to: 



Gjk 



2h 



[{36,k - 1) ± e,k] 



(A34) 



APPENDIX B: REVIEW OF BOUNDARY 
CONFORMAL FIELD THEORY TECHNIQUES 

Here we briefly review the boundary conformal field 
theory (BCFT) techniques, developed largely by J. 
Cardy '35], which have been applied to various quantum 
impurity problems. 

We consider a general conformal field theory, such as a 
collection of free bosons or fermions, defined on the 1/2- 
line, a; > 0, with a conformally invariant boundary condi- 
tion at a; = 0. In general two-point correlation functions 
will be affected by the boundary. While the correlation 
function of some operator, 0{x,t) may behave as: 



1 



[{x - x')^ + (r - r') 



A21A ■ 



(Bl) 



in the bulk (i.e. in the absence of a boundary or far from 
the boundary), in the limit x, x' <^ \t — t'\ it behaves as: 



{0{x,t)0{x' ,t')) cx 



r'|2Af 



(B2) 



Here A is the bulk scaling dimension of the operator, 
O and As is its boundary scaling dimension. In gen- 
eral. As will depend on the boundary conditions. The 
boundary conditions are always assumed to imply: 



P(t,0) = 0, 



(B3) 



where: 



V{t, x) = nnit -x)~HL{t + x). 



(B4) 



V is the momentum density and H-l/r are the left and 
right-moving parts of the Hamiltonian density. 

We can determine all boundary scaling dimensions for 
an arbitrary conformally invariant boundary condition, 
A, denoted A^, from the finite size spectrum of the 
Hamiltonian on a strip of length I with boundary con- 
ditions A at both ends. This is done by the conformal 
transformation: 



Here 



(B5) 



(B6) 



covers the infinite half-plane and w ^ u -\- iv the strip 
with Q < V < I. The correlation function for two points 
at the edge of the strip can be obtained by this conformal 
transformation: 



{0{UI)0{U2)) = 



— smh— (ui - U2) 

TT ZL 



-2A/ 



(B7) 



As U2 — Ui 

{0{UI)0{U2)) 



00, this approaches 
7rN2A 

7 



exp 



-ttA, 



(m2 - Ml) 



I 



(B8) 

On the other hand, since u is imaginary time, we may 
evaluate the correlation function on the strip by inserting 
a complete set of states: 



E 



|(0|O|n)AA|'exp[-^„(u2-ui)]. 



(B9) 

Here |0) is the groundstate and \n) an arbitrary excited 
state for the strip Hamiltonian with b.c.'s A at both ends. 
(The groundstate occurs here since the strip has infinite 
length in the imaginary time, u, direction.) En is the 
energy of the n*'' state (measured from the groundstate 
energy). As U2 — U1 — > 00, the lowest energy excited state 
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that can be created from the groundstate by the operator 
O dominates so we conclude that this state has energy: 

El = n^. (BIO) 

There is a one-to-one correspondance between boundary 
operators and the finite size spectrum with the boundary 
scahng dimensions and finite size energies related by Eq. 

(IbToIi . 

Another very useful quantity to consider is the par- 
tition function with conformally invariant b.c.'s A 
and B at the 2 ends of a finite system of length I, at 
inverse temperature, f3. This may be expressed in terms 
of the finite size spectrum on the strip as: 

Zab = tr expi-Z^i/^-^S] ^ exp[-/3S,f (Bll) 

n 

where the E:^^ are the energies on the strip with b.c.'s 
A and B at the two ends. This picture is sometimes 
referred to as "open string channel" . On the other hand, 
we may interchange space and imaginary time and write 
instead: 

Zab = {A\cM-IH^]\B), (B12) 

where \A) and \B) are boundary states that correspond 
to the b.c.'s A and B respectively. Here denotes the 
Hamiltonian with periodic b.c.'s on a ring of length f3. 
This point of view is sometimes referred to as "closed 
string channel" . All boundary states obey the condition: 

V{x)\A)=0. (B13) 



The eigenstates of are of the form: 

< - ^A„, (B14) 

where the A„ are bulk scaling dimensions. This follows, 
similarly to the correspondance between boundary scal- 
ing dimensions and boundary energies, from a conformal 
transformation from the infinite plane to the cylinder of 
circumference f3. 



APPENDIX C: REVIEW OF BOUNDARY 
CONFORMAL FIELD THEORY FOR STANDARD 
FREE BOSONS 

1. One component free boson 

A simple example of boundary dimensions and energies 
is given by the periodic boson. Thus we consider the bulk 
Lagrangian density: 

'C = £(M', (CI) 

and assume that ip{x, t) is a periodic variable so that 
is identified as 

(p^(p + 2Trn (neZ). (C2) 



I 

These boundary conditions determine the bulk mode expansion for ip{x, t) on a circle of circumference P, < x < (3: 



ip{t,x) =<^n^ — 
P 



P'x + -Pt 
9 



, exp 



2n 

'J 



exp 



27r 



h.c. 



Here: P is the momentum operator conjugate to the constant term, ipg: 

[00, P] = h 



(C3) 



(C4) 



and P' is another momentum operator. They both have integer eigenvalues. The integer eigenvalues of P' follow from 
the periodic b.c.'s and the angular nature of ip: 



ip{l3) = (^(0) + 27rn. 



(C5) 



The integer eigenvalues of P follow from the fact that it is conjugate to an angular variable, ipo. Explicitly, the 
wave- function, exp[— iP</Jo], with P an eigenvalue of P, must be single valued. 



x± = tzL X, 



(C6) 



L /H 

and a„ are boson creation operators for the left and right moving finite momentum modes. This mode expansion 
is consistent with the equal-time commutation relations: 



fix), 



difix') 



dt 



— — 10(X — X ), 
9 



(C7) 
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which follow from the normalization of the Lagrangian. We may decompose ip into left and right moving modes: 



and then write the dual field: 

This has the mode expansion: 
2tt 



9{t,x) 



13 



[Px + gP't] 



an exp 



2n 



'inxj. 



a-n exp 



27r' 



h.c. 



P' is the momentmn conjugate to Oq: 



[Oo,P']^i, 

and we see from this mode expansion that 9(t, x) is also an angular variable: 

0{t,x) ^ e{t,x) + 2Tm (n e Z). 



(C8) 
(C9) 

(CIO) 

(Cll) 
(C12) 



r 



We may equally well write the Lagrangian density in 
terms of 0: 



(C13) 



The bulk primary operators are 

exp{i[rnp{t,x) + m9{t,x)]} with scaling dimension: 



gm 



25 



(n, m G Z) 



(C14) 



(All descendent operators are simply primaries multiplied 
by products of multiple derivatives of ip. These have 
dimensions equal to that of the corresponding primary 
plus a positive integer. ) Inserting the mode expansion 
of if into the Hamiltonian, gives the finite size spectrum 
with periodic b.c.'s on a circle of circumference /3: 



determines (pR{t, x) (for a; > 0) as the analytic continua- 
tion of ipLit,x) to the negative x-axis: 

ipR{t,x) = (po - ipLit,-x), (a; > 0). (C19) 

The correlation function for ipL(t, x): 

(ipL(t,x)(pL(0,0)) = ---^ln.'E+ + constant, (C20) 
2.9 

is unaffected by the boundary. Thus we see that the 
boundary operators can be rewritten by the replace- 
ments: 



0(t,O) 



H 



27r 
T 



p2 

25 



- 2 

P' 9 



(C15) 

We see that the relationship, Eq. (jB14|) between the bulk 
scaling dimensions and finite size energies with periodic 
b.c.'s is obeyed. 

Now consider the Dirichlet (D) b.c. on (p: 

Lp(t,0) — (fio, (a constant, independent of t), (C16) 

implying: 



(C17) 



We may calculate the dimensions of all boundary opera- 
tors directly by recognizing that the D b.c: 



ipR(t,0) = ipo - (/?L(t, 0), 



(C21) 

Thus the non-trivial operators are: 

exp[zm6'(t, 0)] — > exp[2zm(/52,(i, 0) — irrupo], (C22) 
of scaling dimension: 

Ad = gm^. (C23) 

Note that this is twice the hulk scaling dimension of 
exp[zm6']. This factor of 2 arises from the b.c. The 
bulk scaling dimension of exp[im0] consists of equal 
contributions of gm?' /A from the left and right factors: 
exp[±z(y9i/fl]. Upon imposing the D b.c. the right fac- 
tor vanishes and the left factor has an extra factor of 2 
in the exponent which quadruples the scaling dimension: 



(C18) gm74 grr? 
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We may check the general BCFT results on this simple example. Consider the finite size spectrum with the same 
D b.c. on both ends of the strip of length I (with the same value of i^o)- The D b.c. essentially sets P — and 



-a^. The mode expansion of Eq. (|C3p becomes: 



27r - 



1 ^ 1 o 



sm 



m:x\ 



cxp 



innt 



(C24) 



P' must have integer eigenvalues due to the D b.c. and the periodic nature of ip. Thus the finite size spectrum with 
D b.c.'s is: 



~ 2 \ ^ 

9P' +E 



(C25) 



We see that the general relation, Eq. ljB10|l between the dimensions of boundary operators and finite size spectrum 
on a strip with the corresponding boundary conditions is obeyed. 

We may also find the corresponding boundary state. This must obey the operator equation: 



Note that this implies: 



ip{0,x)\D{ipt))) = (po|£'(<y3o)), (independent of x). 
9(^(0, x) 



dx 



\D{^^))=Q. 



(C26) 



(C27) 



Note that, compared to the operator boundary condition of Eq. (jC17|) . t and x have been interchanged. This is a 
consequence of the interchange of space and time involved in going between the two interpretations of Zab ■ In the 
boundary state representation, the boundary corresponds to the circle t — Q and Eq. ljC27|) is the condition that tp 
(when acting on the boundary state) be constant along the boundary. The D boundary state is: 



exp 



exph*P^o]|(0,P)). 



(C28) 



J P=-oo 



Here |(0, P)) is the eigenstate of P' with eigenvalue 0, the eigenstate of P with (integer) eigenvalue P, and the 
groundstate of all the harmonic oscillators. The condition Eq. HC26|) follows using the explicit form of the wave- 
fimction |(0,P)): 



(V3ol(0,P)) oc exp[iPv3o] 



(C29) 



[We apologize for the confusing, but unfortunately standard, notation here. In Eq. (jC29|) . ifQ is a co-ordinate whereas 
in Eq. ljC28|l . (^0 is a fixed number, corresponding to an eigenvalue of the co-ordinate.] It is straightforward to 
calculate Zjjd in the boundary state representation: 



exp 



21 



.27r P 



(C30) 



Here we have introduced the Dedekind //-function: 

00 

l^(q)^q''^^Yl{l-it), (C31) 



11=1 



and the convenient notation: 



q = e . 



(C32) 



The factor of l/r?(g) in Eq. IjCSOjl comes from the oscilla- 
tor mode factor in the boundary state of Eq. ljC28|) and 



we have included the universal groundstate energy for 
free bosons (from the zero point motion/ Casimir effect): 



Eo ^ (C33) 
The sum in Eq. (|C30p can also be written in terms of q: 



(C34) 
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In order to check that Zjjd indeed gives the correct finite 
size spectrum with D b.c.'s we need to re-express it in 
terms of: 



= p. I 



q = e 



(C35) 



i.e. perform a modular transformation. The modular 
transformation of the Dedekind r^-function is: 



(C36) 



v{q) = \I Yi 



The modular transformation of the sum in Eq. ljC30l) can 
be computed using the Poisson summation formula, i.e. 
the Fourier transform of the periodic (5- function, dp{x): 



inserting Eqs. fC36|) and fC37t into Eq. fC30^ gives: 



Z 



DD 



(CSS) 



Using the definition of q in Eq. I|C35|I and the repre- 
sentation Eq. (|B11|I for Zab, we extract the finite size 
energies with D b.c.'s at both ends of the strip: 



^ ,-^V(4.) ^ Y: exp 
Pez Pez 



Trip' 



21 



5/3 
— 

P'ez 
P'ez 



I 



The Neumann (N) b.c. is: 



^^"^ = t[5P''+ integers], 



(C37) in agreement with Eq. lfC25|l 
I 



This is equivalent to 89 /dt = or equivalently: 
This implies: 

'PRi't,^) ^ -9a/g + (pL{t,-x), (x > 0). 
Now the non-trivial primary boundary operators are: 

exp[m(^(i, 0)] cx exp[2zn(pi(t, 0) — in9Q/g], 

of dimension: 

An = n'^/g, 

again twice the dimension of the corresponding bulk operator. The mode expansion is: 



. , 27r - 1 ■c-^ 1 „ /nTrx\ 



a„ exp 



imTt\ 



I J 



h.c. 



(C39) 



(C40) 

(C41) 
(C42) 
(C43) 
(C44) 

(C45) 



with P having integer eigenvalues. The corresponding spectrum can be read from after rewriting the Hamiltonian as: 



ttNN 

=7 



p2 



(C46) 



The corresponding boundary state, obeying: 



9{0,x)\N{9o)) ^9q\N{9o)), 



(C47) 
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l^(^o)> - (f 



1/4 



exp 



E 



J2 e^p[^zP'eo]\iP',0)). 



(C48) 



J P' = -oo 



r 



2. Multi-component free boson 

In some applications, including the subject of the 
present paper, we need to consider a multi-component 
free boson field theory. While Ref. [s^ also contains 
a similar review, here we summarize the basics of the 
multi-component free boson field theory conforming to 
the normalizations of the present paper. Let 



(C49) 



be a c-component free boson field, with the Lagrangian 
density 



(C50) 



At this point, the theory is just a collection of c free 
bosons which are independent of each other. However, 
we will be interested in various possible boundary in- 
teractions which couples different components. In fact, 
we often introduce a multi-dimensional generalization of 
the periodicity (|C2|) . which is often called compactifica- 



tion. Sometimes the compactification ties different com- 
ponents together, so that they cannot be regarded as 
completely independent, even before we introduce an in- 
teraction at the boundary. 

The general form of the compactification would be 
given as the identification 



(f ^ ip + 2t:R, 



(C51) 



where R G A for a Bravais lattice A. The different com- 
ponents would be completely independent only if A is 
rectangular. For a given Bravais lattice A, we can define 
a reciprocal lattice A* so that 



K -ReZ 



(C52) 



for any vectors K G A* and i? S A. [Note that the 
lattices, X]i=i "•i^i ^^'^ Si=i for arbitrary integers 

Ui and rrii, where the Ki are defined in Eq. (|4.1(J|I and 
the Ri in Eq. (|6.55() . are reciprocal.] 



J 



Imposing periodic boundary conditions in the space direction, we can generalize the mode expansion ljC3|l to the 
present case as 



N ^ 27r 



Rx + -Kt 
9 . 



1 °° 1 f 



exp 



2tt 

'T 



exp 



2n 



+ h.c. 



(C53) 



Here we have replaced the zero-mode "momentum" operators with their eigenvalues R E A and K E A* . The dual 
field is defined similarly to Eq. ljC9|l . It has a similar mode expansion to the above, generalizing Eq. (jC10|l . As a 
result, 6 can be regarded as compactified as 



e + 2ttK, 



(C54) 



where K e A* . 

The vacua of the oscillator modes 



\iR,K)) 



(C55) 



are then labelled by K and R which are the eigenvalues of P and P' . 



Generalizing ljC28(l . the Dirichlet boundary state cor- where Vb(A) is the c-dimensional volume of the unit cell 
responding to — ipQ at the boundary is given as of the Bravais lattice A, and 



(C56) 



1(0, X))) = exp 



1(0, X)), (C57) 
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is the bosonic Ishibashi state. The partition function for 
the cyHnder with the same Dirichlet boundary condition 
at both ends (diagonal cyhnder amphtude) is given, gen- 
erahzing (|C30|) . as 

{D{0o)\exp[-lH^]\D{0o)) 



^DD 



In order to find the spectrum of the boundary operators, 
we must modular transform Eq. I|C58(I . For this purpose, 
the multi-dimensional generalization of Eq. ljC37|) is use- 
ful: 



1 



E 



(C59) 

Sometimes it is convenient to renormalize the lattices as 
~^ ^* ^^'^ \/9^ ~^ to obtain the equivalent 
expression 



V£A' 



weA 



Using eq. (|C59|I . we find 



B.eA 



(C60) 



(C61) 



In fact, the prefactor (2g) 



-c/2_ 



Vo(A) 



in the bound- 



ary state (|C56|I . which represents the generally non- 
integer "ground-state degeneracy" was chosen so 
that Eq. ljC61l) satisfies Cardy's consistency condition. 
Namely, the coefficient of Eq. (|C61|I must be unity (or 
integer) to allow the interpretation as in Eq. (|Blip . From 
Eq. (|C61|I we can read off the scaling dimensions of the 
boundary operators with Dirichlet boundary condition as 



A^i = gR^ + integer, 



(C62) 



where i? e A as usual. These are the dimensions of the 
boundary operators exp[iQ ■ R], consistent with the com- 
pactification of Eq. (|C54p . 

Likewise, the Neumann boundary state corresponding 
to ^ = 00 at the boundary is given as 



where 



V^o(A) J2 exp[-ii? • eo]\{R, 0))), 

(C63) 



rga 



|(i?,0))) =exp 



+ E«' 

n=l 



|(i?,0)), (C64) 



is the bosonic Ishibashi state. The cylinder partition 
function for the Neumann boundary condition at both 
ends reads 



ZNNiq) = (.9/2)^/2^0 (A) 



1 



ReA 



(C65) 



in the closed string channel. Modular transforming using 
Eq. HC59() gives the same partition function in the open 
string channel as 



ZNN{q) 



( 1 



Y^q-- 



/a 



(C66) 



A'eA* 



The scaling dimensions of the boundary operators with 
Neumann boundary conditions are now given by 



Aat = h mteger, 

9 



(C67) 



where X G A* . These are the dimensions of the boundary 
operators exp[i$ • K], consistent with the compactifica- 
tionofEq. (jUHJ. 

The discussions in this Appendix naturally form a ba- 
sis for the problem of the junction of three quantum 
wires, which is mapped to the boundary problem of two- 
component free boson. However, there is an important 
modification due to the fcrmionic nature of the electrons. 



This is discussed in Sec. lVIlH before applying the present 
formalism to the junction problem in Sec. Ixl 



APPENDIX D: BOUNDARY CONDITIONS AND 
THE CONDUCTANCE 



Physical properties of the multiple wire junctions can 
be related to the boundary conditions imposed on the 
bosonic fields describing each wire. In preparation 
for discussing conformally invariant boundary conditions 
that correspond to RG fixed points for the 3-wirc junc- 
tion problem, let us show below how the conductivity 
tensor is extracted from the boundary conditions. 

Let us start with a simple example first, that of a single 
quantum wire, described by fields 'p(x, t) and 9[x, r) with 



(Dl) 



or dual action 



S 



Airg 



dr dx {d^ef 



(D2) 



As in Ref. Jj], we calculate the conductance within the 
linear response theory. Applying an AC electric field in 
an interval Q < x < L and taking the DC limit, we obtain 
the Kubo formula for the conductance of the wire as 
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G 



1™ -r~~T 



dx{TrJ{y,T)Jix,0)) 



1 

G= lim ^ — 



dx 



{TrMy,T)jR{x,0)) + {TrJLiy,T)JL{x,0)) 

-{TrJR{y,T)Mx,0)) - {TrJLiy,T)jR{x,0)) 
I 



(D3) 



where the current J{x, r) = —idrOix, t). The current J can be expressed in terms of chiral components a,s J = Jr — Jl, 
with Jr = 89 = [dx — idr)d/'2. and Jj^ ~ 80 = [8^ + i8r)6/2. In terms of the chiral currents, one has 



(D4) 
(D5) 



1. Conductance of an infinite wire 

Let us consider first an infinite wire. In this case, the 
correlation function between the different chiral compo- 
nents of the currents (JrJl) vanishes. Thus we need only 
the "diagonal" contributions. 

{jR{y,T)jR{x,0)) - 8^9{z)9{0)), (D6) 

where z = ir + {y ^ x). Because 

(0(z)0(O)>^-£ln|z|^ (D7) 

we obtain 

{jR{y,r)jR{x,0))=-^^. (D8) 
47r 

To calculate its contribution to the conductance, we first 
perform the Fourier transformation. A simple contour 
integral gives (for lu > 0) 




dre"^^- -2 = 27rw e-"", (D9) 



where u — y — x and H{t) is the Heaviside step function. 
Combining with a similar calculation on the (JlJl) part, 
the DC conductance for the infinite wire reads 

G = 9^^l^ dx[H{x^y) + H{y~x)]=g^. (DIO) 

In this calculation, the length L of the section on which 
the voltage is applied does not affect the final result of 
the conductance. Moreover, the location y where one 
measures the current does not matter, and it can even be 
outside of the region where the voltage is applied. This 
is natural because in the DC limit, the current would be 
uniform throughout the wire even if the voltage is applied 
to a particular section of the wire. 

However, the result (|D10I) should be treated with cau- 
tion. As it has been discussed by several authors [s^, |^ 
Is^ . the observed DC conductance per channel in a quan- 
tum wire is generally the free electron value G = e^//i in- 
dependent of the TL parameter g, rather than the renor- 
malized value eq. ljD10|) . This is related to the fact that 



the voltage is applied in a physical setting corresponds to 
the potential drop between the two reservoirs to which 
the wire is connected. The assumption of the electric field 
applied on the finite interval of the wire does not reflect 
this. Nevertheless, the calculation leading to eq. (|D10|) is 
useful as a starting point for various applications. In fact, 
the physical result G = /h can be recovered within the 
present approach, by modeling the reservoirs by Fermi 
liquid leads (TL liquid with g = 1). 

In this paper, we first ignore the issue of the reservoirs 
and calculate the conductance simply using the present 
approach. Later, in Sec. IXIll we will discuss the physical 
DC conductance by including the Fermi liquid leads. 



2. Conductance of a semi infinite wire 

As the simplest example involving a boundary, let us 
consider a half-infinite wire with an "open" boundary at 
X = Q (through which no current can flow.) The wire 
is assumed to extend for a; > 0. It corresponds to the 
Dirichlet boundary condition on the field 9 or equiva- 
lently the Neumann boundary condition on the field ip. 
The boundary condition 6'(0, r) = is translated to 

Jfl(O) - Jl(0) = 0. (Dll) 

A convenient trick to respect this condition is to ana- 
lytically continue the right-mover current to x < and 
identify 

Jl{x,t) = Jr{-x,t) (Df2) 

for a; > 0. 

As a general rule, correlation functions of the chiral 
operators are not affected by the boundary condition. 
Therefore the contributions from (JrJr) and (JlJl) re- 
mains the same as in eq. (jDIO|) . 

Cross correlations between the different chiral compo- 
nents are generally non-vanishing in the presence of the 
boundary, and are determined by the boundary condi- 
tion. In the present case, all the left-mover current may 
be replaced by the right-mover current (in fictitious, neg- 
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ative X region) according to eq. (|D12p . As a result, 

-2 (D13) 



{jRiy,T)JL{x,0)) = {jR{y,T)jR{-X,0)) 

9 1 



{iT+ (x + y))' 



Upon Fourier transform, this term contributes to the con- 
ductance a term 



Grl — —g- 



hL 



dx H{x + y) = —g- 



(D14) 



because a;, y > 0. On the other hand, the corresponding 
one from {JlJr) vanishes because it involves the step 
function H{—x — y) for x, y > 0. 

Summing up all the contributions, we obtain the DC 
conductance 

G = (D15) 
for the open wire. This is an intuitively obvious result. 



3. Conductance of a half-infinite wire attached to a 
superconductor 

As another simple example, let us consider the half- 
infinite wire attached to a superconductor aX x — 0. 
Again the wire is assumed to extend to a; > 0. Now 
the appropriate boundary condition on the boson field is 
given by the Dirichlet boundary condition on (/>, since it 
fixes the phase of a Cooper pair. In terms of current, it 
is now given as Jr{Q) + Jl{0) = and it is solved by 
identifying 



Jl(x,t) = -Jr{-x,t) 



(D16) 



for X > 0. The calculation exactly follows the previous 
case in Sec. ID 21 except that the contribution from the 
cross term (JrJl) flips the sign. As a result, we find the 
DC conductance 



G = 2g- 



(D17) 



which is doubled from the bulk one. This may be inter- 
preted as a consequence of the (perfect) Andreev reflec- 
tion. 



J 



4. Conductance of the Y-junction 

Now let us turn to our problem of the Y-junction. We 
define three half-infinite wires 1,2 and 3 which extends 
for X > 0, and a; = is represents the junction. Each 
wire is bosonized separately and described by the boson 
fields (fij and 9j . Currents are also defined respectively. 

We discuss the conductance matrix introduced in 
Eq. H2.2|l . In the present approach, the conductance ma- 
trix is given as 



lim 



1 



w^o+ h nujL 
where J, 



dx 



dT e'"^(J,(y,T)Jfe(x,0)), 



J — -idrOj is the current operator for wire j. 

In analyzing the Y-junction, it is convenient to work 
on the rotated basis defined in Eqs. (|4.5I4.6I) . We can 
express the fields ipj and 9j in this basis as 



1 



eo 



(D19) 
(D20) 



where the vectors Kj are defined in Eq. 1)4. 9|l . and recall 
that i = («>i,$2) and 9 = (61,62). 

It also defines the transformation of the currents as 



(D21) 

where summation over ^ is implicitly assumed, 

= -idrO^, (D22) 

and the coefficients can be read off from 

Eqs. (11)1911)2011 : 



1/V3 



/i = 
M= 1,2 



From Eq. (|D18|I . the conductance in this basis is given by the new basis, 



G 



jk 



lim 

^^0+ h TTLOL 



dx 



dT e^'^^ v.^vk. {J^{y,T)r{x,0)), 



(D23) 



where the summation over fi, v is implicitly assumed. 

Conservation of total charge implies that the "center of mass" boson is always subject to the Neumann bound- 
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ary condition. As a consequence, the cross-terms {J'^ J^) 
vanish for ^ ^ 0. Moreover, the contribution from {J'~'J'~') 
to the conductance also vanishes, because the calculation 
is identical to that in the open wire. Therefore, the con- 
ductance of the Y-junction is given by Eq. (|D23|I with 
the /i, sum restricted to over just 1,2. 

The boundary conditions will determine the {Jf^J") 
correlations. We will use the expression Eq. (|D23() ex- 
plicitly in the calculations carried out for the different 
boundary conditions we analyze in Sec. Ixl 



APPENDIX E: JUNCTIONS OF THREE 
FRACTIONAL QUANTUM HALL EDGES 

For the case g = 1 and for flux cj) — ±7r/2, we showed 
that we could understand the scattering in terms of an 
S'-matrix that took incoming/outgoing states cyclically 
from lead 1 — > 2, 2 3, and 3^1. We can gain a 
great deal of intuition beyond the 5=1 case by studying 
the case of tunneling in junctions of three edge states 
of fractional quantum Hall (FQH) liquids, in particular 
those belonging to the Laughlin sequence (filling fraction 
ly = l/(2m -t- 1), TO odd). In this case, the dynamics of 
the edge modes is described by chiral bosons, and the 
fractionally charged edge quasiparticles are constructed 
from these chiral bosons 37] . We will show here that one 
can still understand the scattering in junctions of three 
FQH edges as cycling between the three leads. 

One of the results below is that we can construct the 
Klein factors for the multi-lead problem using the zero- 
modes of the chiral bosons, instead of the Pauli matrices 
used by Nayak et al. 0. 

The geometry we consider is shown in Fig. 1191 The 
shaded areas correspond to Hall liquids, and the white 
regions to vacuum. Electrons tunnel through the vac- 
uum or the FQH liquid, while quasiparticles tunnel only 
through the FQH hquid. 

In order to understand the quasiparticle tunneling pic- 
ture for the three leads, we have to ensure that the quasi- 
particle operators for the three leads respect the appro- 
priate fractional statistics relations: this is the general- 
ization of the Klein factors needed in the fermionic prob- 
lem of Nayak et al.. 



1. Klein factors from the boson zero-modes 



Consider three right moving boson fields 



i)i(t — x) ^ Xi-\ — Y Pi {t — x) + Oscillator modes 

(El) 



•4 

4 





(a) electron tunneling (b) quasiparticle tunneling 

FIG. 19: Two tunneling geometries in a three lead FQH junc- 
tion, (a) Electron tunneling geometry, where the Hall liquid 
is separated in three disjoint pieces (dark regions); tunneling 
occurs through the connected vacuum (white region), so only 
electrons can tunnel between the tips of the three FQH re- 
gions, (b) Quasiparticle tunneling geometry, where there is a 
single connected FQH liquid (dark region) and so fractionally 
charged quasiparticles can tunnel between the three closely 
spaced tips of the liquid. 



for i — 1,2,3 labeling the three leads. The zero- modes 
satisfy [Xi, Pj] = i Jy. The Lagrangian is given by 



>C = - ^ ^ dx<t>i {dt + dx) > 

i=l 



(E2) 



(In this appendix we adopt the standard normalization 
for edge modes, see Ref. [s^l- This normalization is 
related to our standard normalization throughout the 
rest of the paper by taking Att 2tt in the Lagrangian 
Eq.El) 

The quasiparticle/electron operators are written as 



iXXi iX 

e e 



'-Pi{t-x) 



(E3) 



where ipi contains only oscillator modes. (A — 
for electrons and A = ^/v for quasiparticles within the 
Laughlin sequence ly = 2m+i •) These ipi have the correct 
statistics with respect to themselves, but not among each 
other (they simply commute). To fix this, define the 
Klein factors 

= e'^^^"'^-^' , (E4) 
which satisfy [ry^, rjj] — 0, and 

= i)j r]i X e*"^"'^' . (E5) 

If we define new edge operators 

= V^A , (E6) 

we do not spoil the statistics between same edge opera- 
tors as long as we choose au = for i = 1,2,3, so that 
[r]i, tjji] = 0. Now consider the statistical phase when two 
different edge vertex operators '^i and for i ^ j (so 
that [ipi,4'j] = 0) are interchanged: 



^'i^i = A Vj 



(E7) 



61 






FIG. 20: The six possible tunneling between two edges are 
depicted on the left, in terms of vectors that span the trian- 
gular lattice. On the right, the phases accumulated due to the 
Klein factors are shown for elementary up and down triangles, 
according to Eqs. lEQIElOt . 



It sufRces to consider only anti-symmetric matrices so 



that Ui 



-a 



ji, with miit entries a^- 



±1. The right 



statistical angle is picked up if we choose 9X 



ttA or 
and for 

J when z 7^ J. 

Notice that there is still a sign ambiguity in choosing 



9 = TrX: for electrons "iif-^f - 
quasiparticles = ^f^'f 



X e 



the three 



-a-j 



±1 for i 7^ j. To consistently 



fix the correct signs, one must look, for example, at the 
relative phases as compared to those obtained by clo sing 
the system into a single edge, as studied in Ref. 
Doing so, we fix ai2 — a^i = 1 and 023 = — 1. 



2. Tunneling terms and the connection with the 
Callan-Freed model 

Now we turn to the tunneling between the three edges: 
Tij — ^'l^'j. It is useful to also define 



-iX( 



(E8) 



which does not contain the Klein factors. In terms of the 
<I> basis, defined in Eqs. (|4.7() and (|4.8|l . we can write six 
possible operators: 



ti2 = e 



-iXK3<S> 



I23 — e , 



^31 = e 



-iXK2<S> 



with tji — t*j where the Ki are defined in Eqs. H4.10|l . 

The six operators correspond to the six choices of ±Ka, 
a = 1,2,3. 

Now, the Klein factors are responsible for extra phases 
coming when considering time-ordered products of the 
Tij operators instead of the t^'s. 

To illustrate this point, let us now calculate the phases 
obtained when taking products of operators along the 
elementary up and down triangles in the lattice spanned 
by the vectors Ki and K2- 

Using the commutation relations of the rj and ^, one 
can show that 



— ^12*23*31 



(E9) 



so there is no extra phase due to the Klein factors when 
traversing the elementary up triangles (see Fig. I20|l . On 



the other hand, 

= ^13^21*32 X e'^^' 



(ElO) 



when traversing down triangles. 

This is equivalent to a flux flux ttj/"^ in the elec- 
tron tunneling problem (jriy in the quasiparticle tunneling 
problem) for down triangles, and a flux for up triangles 
(See fig.EOll- 

Once we have identified the statistical phases for the 
two elementary triangles in the lattice spanned by the 
vectors Ki and K2, we can match order by order a 
Coulomb gas expansion of electron/quasiparticle tun- 
neling operators with a Coulomb gas expansion in the 
Callan-Freed model. This procedure is explicitly carried 
out in Ref. 22] and in full detail in section IVll 



3. Callan-Freed applied to the FQH junction 

In the case of Fig. lTW a). electron tunneling is irrelevant 
at low energies, so the configuration is stable. Now, for 
Fig. liyr b). quasiparticle tunneling is a relevant perturba- 
tion, so the question is what is the strong coupling fixed 
point that is reached. Here we focus on the Z3 symmet- 
ric case. Since we will concentrate on the quasiparticle 
tunneling problem, we choose A — ^/V. 

Following Ref. ^2] and section IVll we can match the 
Coulomb gas expansion of the FQH tunneling problem 
to that of the Callan-Freed problem model by choosing 
a, (3 in the dissipative Hofstadter as follows: 



a2 +/J2 



and 



47r A 



2Tm 



(Ell) 



(E12) 



where A — \/3/4 is the area of the elementary trian- 
gle of unit length, and n is an integer. The first condi- 
tion Eq. IjElip is fixed by the scaling dimension v of the 
quasiparticle tunneling operator. The second condition 
Eq. l|E12p is fixed by the statistics originating from the 
Klein factors explained above. 

We can then solve for a, (3, and using that the lattice 
dual to the triangular lattice of unit length has lattice 
constant 2 / \/3 (this dual lattice contains the minima of 
the periodic potential in the problem, and the distance 
between the nearby minima are the "instantons"), we ob- 
tain the scaling dimension {2/^/3)^a (of the strong cou- 
pling tunneling operators: 



An = 



4t/ 



3iy2 + {2n-iyy 



(E13) 



As in Ref. '22'|, there appears to be a whole family of 
minima and their associated instanton operators labelled 
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by different integer n. Notice tfiat the larger \2n—v\, the 
more relevant is the operator. 

In what follows, let us consider the cases n = ±1,0. 
One can easily show that, for < 1, the choices ±1 lead 
to relevant operators A„ < 1, so the associated strong 
tunneling fixed points are unstable. Clearly, this will also 
be the case if \n\ > 1, since A„ decreases as \n\ increases. 
So the only choice for a stable fixed point for < 1 is 
n = 0. Notice that when i/ = 1, there are two possible 
choices n = 0, 1 (see and section IVIl for a discussion). 
In this case, we obtain a simple result, 

Ao = - . (E14) 



This result has a very simple physical interpretation: that 
the strong quasiparticle tunneling regime of Fig. I19r bl 
leads to a dual description similar to that in Fig. I19f a') 
where the interior of the original single FQH region is 
pinched, leading to three disconnected FQH leads, be- 
tween which electrons can tunnel. Indeed, Aq = is the 
scaling dimension of electron tunneling operators. 
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